GDL/GMF- 90-01 


Structure  of  Monopropellant 
Spray  Flames  at  Elevated 
Pressures 

by 

G.M.  Faeth,  T.-W.  Lee,  M.E.  Kounalakis 
and  L.-K.  Tseng 


FINAL  REPORT 


15  January  1990 


U.S.  Army  Research  Office 
Contract  No.  DAAL03-86-K-0154 


Department  of  Aerospace  Engineering 
The  University  of  Michigan 


Approved  for  Public  Release;  Distribution  Unlimited 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
authors  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


n  it»aismHa3fcgi3MMg£a 


REPORT  DOCUMENTATION  PAGE 


utrifort  sicurity  classification 

unclaoaiiiod 


O.  SECURITY  CLASSIFICATION  AUTHORITY 


20.  OICLASSIFlCATION/ DOWNGRADING  SCHEDULE 


4.  FIRFORMING  ORGANIZATION  at, OUT  NUMItR(S) 


1  OISTRUIJTION/AVAIUUIUTY  OF  REFORT 

Approved  for  public  release; 
distribution  is  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMSER(S) 


NAM*  OF  RIRFORMING  ORGANIZATION 

Dept,  of  Aerospace  Engineering 
The  University  of  Michigan 


AOORISS  (aty,  J(MI,  A nd  Abeam) 

217  Aerospace  Engineering  Building 
Ann  Arbor,  MI  48109-2140 


?*.  NAM  I  op  MONITORING  ORGANIZATION 

U.S.  Army  Research  Office 


7b.  aooriss  (Cty.  sun.  end  zip  cot*) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


k  nami  of  funding /sponsoring 

ORGANIZATION 

U.S.  Army  Research  Office 


it  aooriss  (CHr.  sun.  on*  zjrcmW 


lb.  OFFICt  SYMIOL 
f It 


9.  PROCURIMINT  INSTRUMINT  IOINTIFICATION  NUMIIR 

DAAL03-86-K-0154 


IUMINT  NO. 


P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


1 1 .  ntu  (InehM »  Soeufrty  CUmfctt 

Structure  of  Monopropellant  Spray  Flames  at  Elevated  Pressures  (U) 


1 2.  PfRSONAL  AUTWORCS) 

G.M.  Faeth,  T.-W.  Lee,  M.E.  Kounalakis  and  L.-K.  Tseng 


13b  TIM*  COVIRiO  II*.  OAT*  OF  RIFORT  (Yoor,  Month,  Ooy) 

FROM  9/1/86  T0  12/3 l/8p  15  January  1990 


II.  SUFFtlMINTARr  NOTATION  The  view,  opinions  and/or  findings  contained  in  .this  report  are  those 
of  the  author (s)  and  should  not  be  construed  as  an  official  Department  od  the  Army  position. 


C0SAT1  COOIS 


GROUP  SUBGROUP 


It.  SUUICT  TIRMS  (Conbnuo  on  royom  H  memory  and  Monttty  by  Mac*  rmmbot) 

Spray  combustion 
Liquid  propellants 
Monopropellant  combustion 


(Cononuo  on  mono  it  memory  ord  Mormty  by  Mock  nttmbor) 


The  combustion  properties  of  HAN-based  monopropelluiu  (LGP184S  and  1846)  were  studied,  considering;  liquid  surfsce 
properties  and  conditions  required  for  the  thermodynamic  critical  point  to  be  reached,  drop  combustion  properties  in  combustion  gis 
environments  at  pressures  of  0.2-7, 0  MPa,  and  pressure-atomized  spray  combustion  properties  in  combustion  gas  environments  at 
pressures  of  3-9  MPa.  Computations  of  liquid  surface  properties  in  flame  environments  indicated  unusually  high  liquid  surface 
temperatures,  ca.  1000K,  and  critical  combustion  pressures,  ca.  230  MPa,  since  these  propellants  have  thermodynamic  properties  similar 
to  molten  salts.  Drop  burning  was  dominated  by  subsurface  reaction  leading  to  bubble  formation  and  microexplosions  which 
mechanically  removed  liquid,  at  low  pressures  (less  than  2.1  MPa);  and  conventional  heteieogeneous  combustion  from  the  surface  at  high 
pressures  (greater  than  2.1  MPa).  The  resulting  apparant  burning  rates  were  relatively  independent  of  pressure  and  were  large,  ca.  10 
mm/s.  which  is  consistent  with  earlier  strand  burning  rate  measurements  of  jelled  propellants  at  pressures  greater  than  10  MPa  by 
McBrauiey  (1980,  1981),  Measurements  of  drop  size  distributions,  liquid  flow  rates,  and  liquid  mass  fluxes  m  combusting  sprays  were 
consistent  with  the  individual  drop  burning  rate  measurements  and  exhibited  strong  effects  of  separated  flow.  Deterministic  and 
stochastic  separated  flow  models  were  developed  which  yielded  predictions  that  were  similar  to  each  other  and  were  in  reasonably  good 
agreement  with  the  measurements.  A  locally-homogeneouj  flow  model  was  also  developed  to  provide  estimates  of  spray  behavior  as 
thermodynamic  critical  conditions  are  approached;  this  approach  indicates  strong  effects  of  the  degree  of  flow  development  at  the 
injector  exit  on  combustion  properties,  with  fully-developed  turbulent  flow  yielding  the  smallest  combustion  volumes. 


20.  DISTRIBUTION /AVAIIAI1UTY  OF  AISTRACT 
3  unclassifuoajnumitiq  3  same  At  RFT 


224.  NAMI  OF  RISPONSIIU  INOIVIOUAL 

David  M.  Mann 


00  Form  1473,  JUN  M 


□  OTiG  USER* 


21  Alt  TRACT  SICURITY  CLASSIFICATION 


L*liT-#fn  M 


220  TEKPHONtf/mftd)  MoCadOi  2*e  .OFFKI  .SYMIOL 

(919)  549-0641 


Frvvious  adit! o/il  try  odnsNfd. 


MION  OF  THIS  aaGE 


ii 


Unclaiaif led 


Table  of  Contents 


Page 

Abstract . . _____ . . . . ii 

List  of  Tables _ _____ _ _ _ ___ _ _ ___________ _ , _ _  iv 

List  of  Figures _ _ _ _ _  v 

Nomenclature _ _  vi 

List  of  Appendices  . . _______ . .  vii 

1 .  Statement  of  Problem _ _ _ _ ______ _  1 

2.  Summary  of  Most- Important  Results  _ _ 3 

2. 1  Supercritical  Combustion _ 3 

2.2  Drop  Combustion _ ___ _  5 

2.3  Spray  Combustion . _______________ . . 11 

2.4  Conclusions _ _ _  21 

3.  List  of  Publications _ 21 

4 .  List  of  Participating  Scientific  Personnel _ 22 

References  . 22 

Appendix  A:  Archival  Publications . . . .  25 


List  of  Tables 


Table  Title  Page 

1  Combustion  Properties  of  HAN-Based  Monopropellants.  From 

Kounalakis  and  Faeth  (1988). . . . 4 


List  of  Figures 

Figure  Caption  Page 

1  Sketch  of  a  combusting  pressure-atomized  monopropellant  spray. _  2 

2  Liquid  and  gas  phase  properties  at  the  liquid  surface  for  combusting 

LGP1845.  From  Kounalakis  and  Faeth  (1988).  _ 6 

3  Effect  of  the  TEAN-water  binary  interaction  parameter  on  the  predicted 

critical  combustion  pressure  of  LGP1845.  From  Kounalakis  and  Faeth 
(1988). _ _____ _ 7 

4  Drop  diameter  histories  at  0.51  and  2.1  MPa.  From  Lee  et  al.  (1989).  9 

5  Apparent  drop  and  strand  burning  rates  of  LGP1845.  From  Lee  et  al. 

(1988).  _ _ _  10 

6  Measured  (Birk  and  Reeves,  1987)  and  predicted  (LUF  method)  liquid 

volume  fractions  along  axis  of  combusting  LGP1846  sprays. _  13 

7  Measured  (present  study  and  Birk  and  Reeves  (1987))  and  predicted  (LHF 

method)  liquid  volume  fractions  along  the  axis  of  combusting  LGP1845 
sprays.  From  Lee  et  al.  (1989). _ _ _  14 

8  Measured  (present  study)  and  predicted  (DSF  and  SSF  methods)  SMD 
variation  with  streamwise  distance  for  combusting  LGP1845  sprays. 

From  Lee  and  Faeth  ( 1 990) .  . . . .  1 5 

9  Measured  (present  study)  and  predicted  (DSF  and  SSF  methods)  drop 

size  distributions  for  combusting  LGP1845  sprays.  From  Lee  and  Faeth 
(1990).  _  _ _ _ _ _  17 

10  Predicted  (DSF  method)  phase  velocity  variation  along  axis  of  combusting 

LGP1845  sprays  at  10  MPa.  From  Lee  and  Faeth  (1990). _ 18 

1 1  Predicted  (LHF  method)  mean  scalar  properties  along  axis  of  combusting 
HAN-based  monopropellant  sprays  for  fully-developed  flow  at  injector 

exit.  From  Lee  et  al.  (1988). _ 19 


12  Predicted  (LHF  method)  mean  scalar  properties  along  axis  of  combusting 
HAN-based  monopropellant  sprays  for  slug  flow  at  injector  exit.  From 
Lee  et  al.  (1988). _ _ _ _  20 


v 


Nomenclature 


Symbol 

c 

d 

L 

P 

SMD 

u 

x 


Description 

reaction  progress  variable 
injector  diameter 
injector  passage  length 
pressure 

Sauter  mean  diameter 
streamwise  velocity 
streamwise  distance 


Subscripts 

c 

P 

0 


liquid  volume  fraction 


centerline  value 
drop  properties 
initial  condition 


Superscripts 


(  )  time-averaged  mean  quantity 

O  Favre-averaged  mean  quantity 


VI 


List  of  Appendices 


A.  Archival  Publications 

A.l  "Combustion  of  HAN-Based  Liquid  Monopropellants  near  the 
Thermodynamics  Critical  Point,"  by  M.E.  Kounalakis  and  G.M.  Faeth. 

A. 2  "Analysis  of  Combusting  High-Pressure  Monopropellant  Sprays"  by  T.-W. 

Lee,  J.P.  Gore,  G.M.  Faeth  and  A.  Birk. 

A. 3  "Separated-Flow  Considerations  for  Pressure-Atomized  Combusting 

Monopropellant  Sprays,"  by  T.-W.  Lee,  L.-K.  Tseng  and  G.M.  Faeth. 

A. 4  "Structure  and  Mixing  Properties  of  Combusting  Monopropellant  Sprays," 

by  T.-W.  Lee  and  G.M.  Faeth. 


vii 


1 .  Statement  of  Problem 


Combusting  monopropellant  sprays  are  an  important  fundamental  combusting  flow, 
since  they  are  the  premixed  counterpart  of  more  widely  studied  spray  diffusion  flames. 
Additionally,  combusting  monopropellant  sprays  have  applications  to  regenerative  liquid- 
propellant  guns,  throttable  thrustors,  and  underwater  propulsion  systems.  Motivated  by 
these  observations,  the  present  investigation  considered  monopropellant  combustion 
properties  both  as  individual  drops  and  pressure-atomized  sprays.  The  investigation  was 
limited  to  hydroxyl-ammonium-nitrate  (HAN)-based  monopropellants  (LGP1845  and 
1846)  that  are  of  interest  for  several  high-pressure  monopropellant  combustion  systems. 

The  general  nature  of  a  combusting  pressure-atomized  monopropellant  spray  is 
illustrated  in  Fig.  1.  This  involves  injection  of  liquid  from  a  passage  into  a  combustion 
chamber  where  gas  velocides  are  modest  near  the  injector  exit.  Heat  losses  from  practical 
monopropellant  combustion  systems  are  small;  therefore,  gas  temperatures  are  near  the 
adiabatic  combustion  temperature  of  the  propellant.  The  atomization  breakup  regime  is 
illustrated,  where  drops  begin  to  form  from  the  surface  of  the  liquid  right  at  the  injector 
exit,  since  this  condition  is  most  important  in  practice.  Atomization  breakup  yields  a  liquid 
core,  much  like  the  potential  core  of  a  single-phase  jet,  surrounded  by  a  multiphase  mixing 
layer  containing  drops  and  other  irregular  liquid  elements  (Ruff  et  al.,  1989,  1990). 
Beyond  the  end  of  the  liquid  core,  the  flow  becomes  a  dilute  spray  which  has  a  finite  length 
due  to  the  combustion  of  the  drops.  Since  monopropellant  burning  rates  are  relatively  high 
(with  liquid  regression  rates  never  much  lower  than  lOmm/s  for  HAN-based 

monopropellants)  flame  thicknesses  over  liquid  surfaces  are  small  (ca.  l|im)  so  that  the 
actual  combustion  process  occurs  in  thin  heterogeneous  reaction  zones  which  cover  all 
liquid  surfaces  —  separating  liquid  reactants  and  combustion  product  gases  (Lee  et  al. 
1988,  1989;  Lee  and  Faeth,  1990).  Major  issues  with  respect  to  combusting 
monopropellant  sprays  involves  the  extent  of  the  liquid-containing  region  and  the  rate  of 
energy  release,  so  that  adequate  combustion  volumes  can  be  provided;  the  presence  and 
nature  of  combustion  instabilities;  and  requirements  for  ignition  and  extinction  of  the 
spray.  None  of  these  issues  can  be  addressed,  however,  without  some  understanding  of 
drop  combustion  properties  and  spray  structure;  therefore,  these  simpler  problems  were 
the  main  focus  of  the  present  investigation. 

Several  questions  concerning  monopropellant  drops  and  spray  combustion  must  be 
resolved  for  an  adequate  understanding  of  the  process.  First  of  all,  monopropellant 
combustors  operate  over  a  wide  range  of  relatively  high  pressures,  which  implies  that 
liquid  surfaces  can  approach  or  exceed  their  thermodynamic  critical  point  (Faeth,  1972, 
1987;  Kounalakis  and  Faeth,  1988).  Combustor  pressures  when  the  critical  point  is 
reached  must  be  identified,  since  below  this  pressure  a  multiphase  flow  like  Fig.  1  is 
present,  while  above  it  the  flow  becomes  a  single-phase  dense  fluid  flow  with  very 
different  properties.  Secondly,  the  combustion  properties  of  individual  drops  must  be 
known,  for  the  pressure  range  where  drops  are  present,  since  drop  combustion  mainly 
controls  combustion  rates  of  the  spray.  Next,  information  concerning  the  structure  of 
combusting  sprays  is  needed;  particularly  the  extent  of  the  liquid-containing  region,  drop 
sizes  that  are  present,  and  the  importance  of  separated-flow  phenomena  like  the  relative 
velocities  between  the  drops  and  the  gas  and  the  combustion  rates  of  drops  within  a 
combusting  spray.  Finally,  the  wide  range  of  operating  conditions  for  combusting 
monopropellant  sprays  precludes  exhaustive  testing  to  find  spray  structure  for  all 
conditions  —  effective  models  of  the  process  are  needed. 

The  present  investigation  sought  to  make  a  contribution  to  each  of  these  issues, 
with  the  following  specific  objectives: 
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COMBUSTION  PRODUCTS 


Figure  1.  Sketch  of  a  combusting  pressure-atomized  monopropellant  spray0 


1 .  Liquid  surface  properties  and  pressures  required  for  the  liquid  surface  to  exceed  the 
thermodynamic  critical  point  (the  critical  combustion  pressure)  were  investigated 
theoretically. 

2.  The  combustion  of  individual  drops  in  combustion  gas  environments  was  observed 
at  pressures  of  0.2-7  MPa,  yielding  drop  burning  rates  over  this  range  of 
conditions. 

3.  The  combustion  of  pressure-atomized  sprays  in  combustion  gas  environments  was 
observed  at  pressures  of  3-9  MPa,  yielding  measurements  of  liquid  volume 
fractions,  drop  size  distributions,  and  liquid  mass  flow  rate  and  mass  flux 
distributions.  Additionally,  models  of  the  process  were  developed,  both  ignoring 
and  considering  effects  of  separated  flow,  and  evaluated  using  the  new 
measurements. 

Supercritical  combustion,  drop  combustion  and  spray  combustion  are  discussed  in  turn  in 
the  next  chapter,  followed  by  some  concluding  remarks.  Subsequent  chapters  summarize 
publications  and  personnel  of  the  investigation.  The  present  discussion  is  brief;  more 
information  can  be  found  in  Kounalakis  and  Faeth  (1988),  Lee  et  al.  (1988, 1989)  and  Lee 
and  Faeth  (1990),  which  appears  as  Appendices  A.1-A.4.  Finally,  Lee  (1990)  reports 
additional  details  of  both  measurements  and  predictions  as  well  as  a  complete  tabulation  of 
data. 

2.  Summary  of  Most-Important  Results 
2. 1  Supercritical  Combustion 

The  objective  of  this  portion  of  the  investigation  was  to  find  gas  and  liquid  phase 
properties  at  the  liquid  surface  for  combusting  HAN-based  monopropellants  at  various 
pressures.  The  pressures  required  for  the  surface  to  reach  its  thermodynamic  critical 
pressure,  the  critical  combustion  pressure,  was  also  of  interest  since  this  defines  the 
pressure  regime  where  conventional  spray  processes  must  be  considered.  The  details  of 
this  work  are  described  by  Kounalakis  and  Faeth  (1988),  which  appears  in  Appendix  A.l. 

Hetereogeneous  monopropellant  flames  at  high  pressures  are  very  thin,  ca.  lpm, 
and  will  probably  never  be  resolved  by  measurements.  Therefore,  liquid  surface  properties 
were  found  theoretically,  using  methods  developed  earlier  in  this  laboratory  (Faeth,  1972, 
1987).  This  involved  considering  transport  in  the  gas  phase,  at  the  asymptotic  limit  of 
infinitely  high  activation  energies  of  all  reactions  so  that  reaction  occurs  within  a  thin  flame 
sheet  located  at  a  distance  from  the  liquid  surface.  Transport  in  the  gas  phase  allowed  for 
real  gas  effects  (but  only  considering  concentration  diffusion  with  equal  binary  diffusivities 
of  all  species)  and  nonunity  Lewis  numbers.  Phase  equilibrium  was  assumed  at  the  liquid 
surface  with  real-gas  effects  primarily  treated  using  the  Soave  equation  of  state. 

Combustion  product  properties  were  needed  for  the  computations.  They  were 
found  assuming  thermodynamic  equilibrium  at  the  downstream  end  of  the  flame,  and 
adiabatic  combustion,  using  the  Gordon  and  McBride  (1971)  algorithm.  The  resulting 
properties  of  the  reactants  and  the  combustion  products  are  summarized  in  Table  1  for 
LGP1845  and  1846  at  various  pressures.  Adiabatic  flame  temperatures  are  relatively  low, 
ca.  2000-2200  K;  therefore,  effects  of  dissociation  are  small  and  the  composition  and 
temperature  of  the  combustion  products  varies  very  little  over  the  range  1-100  MPa.  The 
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Table  1.  Combustion  Properties  of  HAN-based  Monopropellants3 
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aAdiabatic  constant  pressure  combustion  with  the  reactant  at  298.15  K. 
bRe actant  composition  (%  by  mass):  (HAN),  63.2;  (TEAN),  20;  and  H2O,  16.8. 
cReactant  composition  (%  by  mass):  (HAN),  60.8;  (TEAN),  19.2;  and  H2O,  20.0. 
dMajor  species  only.  Minor  species  include  CO,  H2,  NO,  OH  and  O2. 


combustion  products  involve  unusually  large  concentrations  of  water  vapor,  with  the 
remainder  being  split  about  equally  between  carbon  dioxide  and  nitrogen. 

Uncertainties  in  predicdons  of  liquid  surface  properties  are  large  due  to  the 
uncertainties  of  the  thermochemical  and  transport  properties  of  the  HAN-based 
monopropellants.  Best-estimate  predictions  of  liquid  surface  properties  for  LGP1845  are 
plotted  as  a  function  of  pressure  in  Fig.  2  (results  for  LGP1846  were  similar  within 
property  uncertainties).  Gas-phase  compositions  are  relatively  independent  of  pressure  but 
liquid-phase  concentrations  vary  substantially,  due  to  increased  liquid-phase  solubilities  for 
combustion  products  as  the  thermodynamic  critical  point  is  approached.  An  interesting 
feature  of  these  results  is  that  the  concentration  of  TEAN  at  the  liquid  surface  is  unusually 
high  at  low  pressures,  due  to  its  low  volatility,  even  though  its  initial  concentration  in  the 
propellant  is  only  20%  (see  Table  1).  This  causes  liquid  surface  temperatures  to  be  rather 
high  at  low  pressures,  ca.  1000  K,  before  decreasing  to  ca.  800  K  as  the  critical  point  is 
approached.  Such  high  temperatures  are  rarely  encountered  at  the  liquid  surface  for  liquid 
monopropellants  and  bipropellants;  they  occur  for  HAN-based  monopropellants  since 
these  materials  have  thermochemical  properties  like  molten  salts. 

Approach  to  the  thermodynamic  critical  point  is  indicated  by  gas  and  liquid 
concentrations  approaching  one  another  at  high  pressures  in  Fig.2.  Predictions  ultimately 
yield  a  critical  combustion  pressure  of  2500  atm.  (250  MPa),  with  an  uncertainty  of  50 
percent  due  to  property  uncertainties.  Similar  to  liquid  surface  temperatures,  the  critical 
combustion  pressures  of  the  HAN-based  propellants  are  unusually  high  in  comparison  to 
other  monopropellants  and  bipropellants,  which  have  critical  combustion  pressures  of  100- 
200  atm.  (Faeth,  1971,  1987).  Thus,  HAN-based  monopropellants  involve  effects  of 
drops  and  sprays  over  much  of  their  range  of  application,  except  for  perhaps  the  highest- 
pressure  portions  of  liquid  gun  cycles.  Again,  this  unusual  behavior  follows  from  the 
molten  salt  character  of  the  monopropellant. 

Examination  of  the  sensitivity  of  critical  pressure  predictions  to  property 
uncertainties  indicated  that  uncertainties  were  dominated  by  the  binary  interaction  parameter 
of  TEAN  and  water  in  the  Soave  equation  of  state.  The  effect  of  this  parameter  of  critical 
combustion  pressure  predictions,  as  well  as  the  effect  of  using  the  Redlich-Kwong 
equation  of  state,  is  illustrated  in  Fig.  3.  Use  of  the  Redlich-Kwong  equation  of  state  is  felt 
to  be  less  reliable  than  the  Soave  equation  of  state  and  this  approach  will  not  be  considered 
any  further  (Kounalakis  and  Faeth,  1988).  However,  estimates  of  critical  combustion 
pressures  are  clearly  quite  sensitive  to  binary  interaction  parameter  TEAN  and  water  vapor. 

In  summary,  the  present  study  of  liquid  surface  properties  during  the  combustion  of 
HAN-based  monopropellants  has  shown  the  following:  results  for  LGP1845  and  1846  are 
similar;  critical  combustion  pressures  are  unusually  high,  250  MPa  with  an  uncertainty  of 
50  percent;  liquid  surface  temperatures  are  unusually  high,  800-1000  K;  liquid  surface 
properties  tend  to  be  dominated  by  TEAN  at  low  pressures,  due  to  its  low  volatility;  and 
the  uncertainties  of  the  predictions  are  dominated  by  the  TEAN-water  vapor  binary 
interaction  parameter.  The  major  conclusion  with  respect  to  spray  combustion  is  that 
combustion  occurs  by  conventional  monopropellant  spray  processes,  involving  liquid 
surfaces  and  drops,  at  all  but  the  very  highest  pressures  of  liquid  gun  cycles. 

2.2  Drop  Combustion 

With  the  importance  of  drop  and  spray  combustion  established  for  the  HAN-based 
monopropellants,  the  next  phase  of  the  investigation  considered  the  combustion  properties 
of  individual  drops.  Drop  combustion  has  been  considered  earlier  by  Beyer  (1986,  1988), 
Beyer  and  Teague  (1986)  and  Zhu  and  Law  (1987),  who  observed  significant  subsurface 
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Figure  2.  Liquid  and  gas  phase  properties  at  the  liquid 

surface  for  combusting  LGP1845.  From  Kounalakis 
and  Faeth  (1988) . 


CRITICAL  COMBUSTION  PRESSURE  (ATM.) 


reaction  leading  to  bubble  formation  and  growth,  and  shattering  of  drops  by  bursting 
bubbles.  However,  these  results  were  either  found  in  low  temperature  gas  environments, 
Beyer  and  coworkers,  or  at  atmospheric  pressure,  Zhu  and  Law,  so  that  the  relevance  of 
the  findings  to  practical  combusting  monopropellant  sprays  was  uncertain.  Additionally, 
McBratney  (1980,  1981)  and  Vosen  (1988)  reported  measurements  of  strand  burning  rates 
of  HAN-based  monopropellants,  which  can  be  related  to  drop  burning  rates.  However, 
McBratney  (1980,  1981)  used  jelled  propellants  so  that  the  jelling  agent  could  have 
influenced  his  measurements;  while  Vosen  (1988)  encountered  liquid  surface  instabilities 
with  unjelled  propellants  which  precludes  using  these  data  to  estimate  drop  burning  rates. 
Thus,  the  present  investigation  considered  unjelled  drops  in  high  temperature  and  pressure 
environments.  The  details  of  this  work  are  described  by  Lee  et  al.  (1989),  which  can  be 
found  in  Appendix  A. 3. 

Drop  combustion  was  observed  with  drops  supported  in  the  post-flame  region  of  a 
premixed  flame  stabilized  within  a  pressure  vessel  at  pressures  of  0.2-7  MPa.  The 
combustion  temperature  of  the  premixed  flame  was  generally  in  the  range  2200-2300  K, 
which  is  slightly  higher  than  the  adiabatic  constant  pressure  flame  temperature  of  the 
propellants  (see  Table  1).  The  drops  were  supported  on  quartz  fibers  (with  initial  drop 

diameters  in  the  range  300-1200  p.m)  within  a  retractable  shield  to  protect  them  from 
disturbances  when  the  premixed  flame  was  ignited.  Once  the  premixed  flame  had 
stabilized,  the  shield  was  retracted  which  submerged  the  drops  in  the  combustion  gas.  The 
drops  were  backlighted  and  photographed  with  a  high-speed  motion  picture  camera. 

Some  typical  photos  of  drop  diameter  as  a  function  of  time  are  illustrated  in  Fig.  4. 
The  origin  of  these  photos  is  arbitrary  due  to  uncertainties  concerning  the  time  that  the 
drops  were  actually  submerged  in  the  flame,  due  to  disturbances  of  the  flame  by  the 
retractable  shield.  Results  at  0.51  MPa  are  representative  of  low  pressure  behavior  where 
subsurface  bubble  formation  and  bursting,  and  mechanical  removal  of  liquid  by 
microexplosions,  caused  erratic  diameter  variations.  This  behavior  is  consistent  with 
observations  of  Beyer  (1986,  1988)  and  Beyer  and  Teague  (1986)  for  lower  ambient 
temperatures  and  Zhu  and  Law  (1987)  for  lower  ambient  temperatures  and  pressures.  This 
behavior  follows  from  the  large  liquid  surface  temperatures  discussed  in  connection  with 
Fig.  2,  and  the  fact  that  thermal  diffusivities  are  higher  than  mass  diffusivities  in  liquids. 
This  allows  subsurface  liquid  having  high  concentrations  of  HAN  to  become  heated, 
leading  to  bubble  formation  by  either  vaporization  or  reaction.  Results  at  higher  pressures, 
e.g.  2.1  MPa  in  Fig.  4,  exhibit  more  conventional  monopropellant  behavior  with  little 
bubble  formation  in  the  liquid  and  gasification  occurring  primarily  at  the  drop  surface  in 
conjunction  with  a  hetereogeneous  flame.  This  occurs  since  surface  regression  rates  are 
high  at  high  pressures,  so  that  the  thermal  wave  cannot  propagate  very  fax  into  the  liquid, 
and  liquid  residence  times  at  high  temperatures  become  small. 

The  variation  of  drop  diameter,  excluding  cases  where  drops  burst  completely  or 
fell  off  the  quartz  probe  between  consecutive  pictures,  was  averaged  over  several  tests  to 
give  effective  drop  burning  rates.  These  results  are  illustrated  in  Fig.  5,  along  with  the 
strand  burning  rate  measurements  of  McBratney  (1980,  1981)  and  Vosen  (1988).  The 
combination  of  drop  bursting  at  low  pressures,  and  conventional  hetereogeneous 
combustion  at  high  pressures,  causes  present  drop  burning  rates  to  be  relatively 
independent  of  pressure.  They  were  relatively  independent  of  diameter  as  well,  which  is 
consistent  with  a  classical  linear  monopropellant  drop  burning  rate  law,  e.g.  -dTp/dt  =  K  = 
const.  The  present  results  are  a  crude  extension  of  McBratney's  (1980,  1981) 
measurements  with  jelled  propellants  suggesting  that  use  of  a  jelling  agent  did  not  have  a 
strong  effect.  Vosen's  (1988)  results  are  much  higher,  due  to  the  effects  of  liquid  surface 
disturbances,  noted  earlier. 
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and  strand  burning  rates  of  LGP1845. 


In  summary,  present  measurements  indicate  drop  burning  rates  of  ca.  10  mm/s  for 
drop  diameters  of  300-1200  |im  and  pressures  of  0.2-7  MPa.  Present  results  are  consistent 
with  earlier  strand  burning  rate  measurements  of  McBratney  (1980,  1981)  at  pressures 
greater  than  10  MPa.  Taken  together,  these  findings  suggest  a  relatively  weak  effect  of 

pressure  on  the  burning  rate,  ca.  p1/3,  except  at  the  lowest  pressures  where  subsurface 
gasification  causes  bubble  formation  and  microexplosions  that  mechanically  remove  liquid 
from  the  drop  at  a  relatively  high  rate.  However,  these  low  pressure  effects  are  not  very 
important  in  practice,  since  it  is  difficult  to  reliably  ignite  HAN-based  monopropellant 
sprays  at  pressures  below  3  MPa  (Lee  et  al.  1989). 

2.3  Sprav  Combustion 

The  final  phase  of  the  investigation  involved  consideration  of  combusting  HAN- 
based  pressure-atomized  sprays.  Birk  and  Reeves  (1987)  reported  earlier  flow 
visualization  results  for  flows  of  type.  During  the  present  investigation,  this  information 
was  processed  to  yield  liquid  volume  fraction  distributions.  In  addition,  a  spray 
combustion  apparatus  was  developed  and  used  for  measurements  of  liquid  volume 
fractions,  drop  size  distributions  and  liquid  flow  rate  and  liquid  flux  distributions.  Finally, 
several  models  of  combusting  monopropellant  sprays  were  developed  and  evaluated  using 
the  new  measurements,  as  follows:  a  locally-homogeneous  flow  (LHF)  model,  where 
interphase  transport  rates  were  assumed  to  be  infinitely  fast;  a  deterministic  separated  flow 
(DSF)  model,  where  finite  interphase  transport  rates  were  considered  but  drop/turbulence 
interactions  were  ignored;  and  a  stochastic  separated  flow  (SSF)  model  where  both  finite 
interphase  transport  rates  and  drop/turbulence  interactions  were  considered.  The  details  of 
this  work  are  described  by  Lee  et  al.  (1988,  1989)  and  Lee  and  Faeth  (1990),  which  appear 
in  Appendices  A.2-A.4. 

The  experiments  involved  observations  of  combusting  sprays  in  a  combustion  gas 
environment  at  pressures  of  3-9  MPa.  The  test  arrangement  was  similar  to  Birk  and 
Reeves  (1987)  but  the  combustion  chamber  had  a  larger  internal  volume.  The  combustion 
gas  was  produced  by  filling  the  chamber  with  a  combustible  mixture  and  then  igniting  it 
with  a  spark  to  reach  the  desired  final  pressure  and  temperature.  The  propellant  was 
pressure  fed  through  various  injectors  into  the  hot  combustion  products,  where  it  ignited 
and  burned  reliably  at  pressures  greater  than  2.7  MPa.  Combustion  gas  temperatures  were 
somewhat  greater  than  the  adiabatic  flame  temperature  for  constant  pressure  combustion  of 
the  monopropellant  (see  Table  1).  Measurements  involved  flash  shadowgraph  motion 
pictures  which  were  reduced  to  find  liquid  volume  fractions;  and  slide  impaction 
measurements  of  drop  size  distributions,  liquid  flow  rates  and  liquid  fluxes. 

Under  the  LHF  approximation,  relative  velocities  between  the  phases  are  assumed 
to  be  small  and  the  flow  is  treated  like  an  equivalent  single-phase  fluid  so  that  properties 
that  are  difficult  to  estimate,  like  initial  drop  size  distributions,  play  no  role  in  the 
predictions.  This  type  of  mixing-limited  behavior  is  expected  to  work  best  when  ambient 
pressures  and  liquid  injection  velocities  are  high,  which  yields  relatively  small  drops  after 
breakup  (Ruff  and  Faeth,  1990).  Similarly,  a  mixing-limited  treatment  of  premixed 
combustion  was  used,  assuming  a  thin  flame  sheet  under  the  laminar  flamelet  concept  of 
Bray  (1980).  Turbulent  mixing  was  treated  using  a  Favre-averaged  turbulence  model  that 
had  been  calibrated  using  measurements  from  a  variety  of  constant  and  variable  density 
turbulent  single  and  multiphase  jets  (Bilger,  1976;  Faeth,  1987;  Jeng  and  Faeth,  1984;  Ruff 
et  al.,  1989,  1990). 
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Initial  evaluation  of  LHF  predictions,  using  measurements  of  liquid  volume 
fractions  obtained  from  the  experiments  of  Birk  and  Reeves  (1987),  was  quite  promising 
(Lee  et  al.,  1988).  The  results  are  illustrated  in  Fig.  6.  Measured  liquid  volume  fractions 
have  a  line  of  sight  bias,  however,  this  effect  becomes  small  near  the  downstream  end  of 
the  visible  liquid-containing  region  so  the  method  provides  some  indication  of  the 
penetration  of  liquid  into  the  flow.  Since  the  degree  of  flow  development  at  the  injector  exit 
was  not  well  known,  predictions  were  carried  out  for  a  variety  of  injector  passage  lengths 

ranging  from  slug  flow,  L/d  =  0,  to  fully-developed  turbulent  pipe  flow,  L/d  =  ». 
Predictions  and  measurements  are  in  fair  agreement  for  L/d  in  the  range  2-10,  which  was 
reasonable  for  the  test  conditions.  Predictions  and  measurements  were  also  relatively 
independent  pressure.  Finally,  the  predictions  are  clearly  significantly  influenced  by  flow 
development  lengths  (L/d)  —  this  behavior  was  confirmed  by  independent  measurements 
for  nonevaporating  sprays  at  atmospheric  pressure  (Ruff  et  al.,  1989,  1990). 

Subsequendy,  the  LHF  predictions  were  evaluated  using  the  present  measurements 
(Lee  et  al.,  1989).  Results  analagous  to  Fig.  6  are  illustrated  in  Fig.  7.  The  measurements 
include  the  results  obtained  from  Birk  and  Reeves  (1987),  dark  symbols,  as  well  as  the 
present  measurement,  open  symbols,  which  largely  involved  nearly  fully-developed 
turbulent  flow  at  the  injector  exit.  Predictions  are  for  fully-developed  flow,  either  ignoring 
line-of- sight  bias  effects,  dashed  line,  considering  them  using  a  stochastic  approach 
developed  by  Faeth  et  al.  (1988),  solid  line.  As  noted  earlier,  however,  estimates  of  the 
downstream  end  of  the  liquid-containing  region  are  not  strongly  effected  by  line-of-sight 
bias.  Clearly,  present  measurements  yielded  a  much  longer  liquid-containing  region  than 
Birk  and  Reeves  (1987).  This  was  attributed  to  potential  effects  of  unstable  injector  flow 
or  cavitation  during  the  earlier  measurements.  Comparing  LHF  predictions  with  the  new 
measurements  yielded  more  conventional  behavior  for  the  predictions  based  on  this 
approximation;  namely,  that  the  LHF  method  generally  overestimates  the  rate  of 
development  of  sprays  (Faeth,  1987).  Single -drop  trajectory  calculations,  using  the  LHF 
predictions  to  define  the  velocity  field,  confirmed  this  behavior  and  also  indicated  the 
presence  of  significant  amounts  of  liquid  beyond  the  visible  end  of  the  liquid-containing 
region. 


The  separated  flow  models  were  developed  motivated  by  these  findings.  They  also 
involved  use  of  the  thin-flame  approximation  which  implies  constant  density  gas  and  liquid 
phases.  The  approach  involved  Eulerian  calculations  for  gas-phase  properties,  with 
turbulent  mixing  treated  similar  to  the  LHF  approach  but  allowing  for  drop  source  terms  in 
the  governing  equations  for  conservation  of  mass  and  momentum.  Lagrangian  calculations 
of  drop  trajectories  were  used  for  the  liquid  phase,  specifying  a  constant  drop  burning  rate 
of  lOmm/s  based  on  the  single  drop  measurements  (this  burning  rate  is  sufficiently  high 
that  forced  convection  and  curvature  have  negligible  effects  on  the  burning  rate).  Initial 
drop  size  and  velocity  distributions  were  specified  near  the  injector  exit,  based  on  velocity 
and  liquid  flux  estimates  from  the  LHF  predictions,  and  drop  size  distributions  extrapolated 

from  the  measurements.  Notably,  initial  drop  sizes  (SMD  of  240  pm)  were  consistent  with 
estimates  of  aerodynamic  breakup  (Reitz  and  Bracco,  1982)  or  secondary  breakup  theories, 
after  allowing  for  effects  of  turbulent  flow  at  the  injector  exit  (Lee  and  Faeth,  1990).  The 
DSF  approach  involved  computing  drop  trajectories  based  on  mean  gas  velocities.  The 
SSF  approach  involved  random  walk  calculations  with  the  drops  interacting  with 
instantaneous  velocities  found  from  random  selections  satisfying  the  local  moments  of 
velocity  fluctuations. 

Predicted  and  measured  SMD  as  a  function  of  distance  from  the  injector  are 
illustrated  in  Fig.  8  for  injector  diameters  of  0.3  and  0.6  mm.  In  terms  of  x/d,  the  small 
injector  has  a  much  larger  liquid-containing  region.  This  is  indicative  of  strong  separated 
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Figure  6.  Measured  (Birk  and  Reeves,  1987)  and  predicted  (LHF  method)  liquid  volume  fractions 
along  axis  of  combusting  LGP1846  sprays. 
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Figure  7.  Measured  (present  study  and  Birk  and  Reeves  (1987))  and  predicted  ( LHF  method)  liquid 
volume  fractions  along  the  axis  of  combusting  LGP1845  sprays.  From  Lee  et  al.  (1989) 
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Figure  8.  Measured  {present  study)  and  predicted  (DSF  and  SSF 
methods)  SMD  variation  with  streamwise  distance  for 
combusting  LGP1845  sprays.  From  Lee  and  Faeth  (1990) 


flow  effects,  e.g.,  drop  sizes  and  velocities  were  similar  for  the  two  injectors  and  the  drops 
survive  for  similar  distances,  implying  a  larger  x/d  for  the  smaller  injector.  The  DSF  and 
SSF  predictions  are  not  very  different,  in  view  of  uncertainties  in  specifying  initial 
conditions;  the  DSF  method  is  recommended  due  to  its  reduced  computational 
requirements.  Finally,  predictions  are  in  reasonably  good  agreement  with  measurements, 
implying  that  the  present  single  drop  burning  rate  measurements  are  consistent  with 
burning  rates  in  the  combusting  sprays,  since  the  predictions  are  very  sensitive  to  this 
parameter. 

The  evolution  of  predicted  and  measured  drop  size  distributions  for  the  0.6  mm 
diameter  injector  are  illustrated  in  Fig.  9  (results  for  d  =  0.3  mm  are  similar). 
Measurements  indicate  a  rather  broad  distribution  of  drop  sizes  at  x/d  =  320,  with  the  drops 
becoming  progressively  concentrated  at  small  sizes  as  combustion  proceeds.  Predictions 
were  obtained  by  extrapolating  results  at  x/d  =  320  to  the  injector  exit  to  specify  initial 
conditions;  therefore,  agreement  at  this  point  is  not  surprising.  However,  predictions  still 
treat  the  subsequent  evolution  of  the  drop  size  distribution  reasonably  well. 

The  reason  for  the  importance  of  separated  flow  effects  can  be  seen  from  the 
predictions  illustrated  in  Fig.  10.  Favre-averaged  drop  velocities  and  time-averaged  gas 
velocities  (which  are  the  same  as  the  Favre  averages  when  the  density  is  constant)  are 
plotted  as  a  function  of  distance  from  the  injector  for  d  =  0.3  and  0.6mm.  Velocity 
differences  between  the  phases  are  large,  except  near  the  injector  exit  where  the  LHF 
approximation  is  reasonable  (Ruff  et  al.,  1989,  1990)  and  near  the  end  of  the  liquid- 
containing  region  where  drops  are  small  and  responsive.  Relative  velocities  are  larger  for 
the  small  injector  diameter  since  it  has  faster  deceleration  rates  due  to  its  smaller  scale.  In 
view  of  these  results,  it  is  not  surprising  that  use  of  the  LHF  approximation  was  not  very 
successful  for  present  test  conditions. 

Increasing  pressures  do  not  improve  the  response  of  a  particular  drop  size  to  gas 
motion,  however,  higher  ambient  pressures  reduce  drop  sizes  upon  breakup  so  that  the 
LHF  approach  becomes  more  effective  (Ruff  and  Faeth,  1990).  Reduced  surface  tension 
as  the  thermodynamic  critical  point  is  approached  also  acts  to  reduce  drop  sizes  and 
improve  LHF  predictions.  Thus,  it  is  worthwhile  to  consider  some  LHF  results  relating  to 
high  pressure  combusting  HAN-based  monopropellants. 

Some  typical  LHF  predictions  are  illustrated  in  Figs.  11  and  12.  This  involves 
Favre-averaged  reaction  progress  variable  (which  is  0  in  the  reactants  and  1  in  the  products) 
and  time-averaged  liquid  volume  fractions  along  the  spray  axis  for  slug  and  fully-developed 
turbulent  flow  at  the  injector  exit.  At  the  LHF  limit,  effects  of  Reynolds  number  are  small; 
thus,  pressure  is  the  only  parameter  on  the  plots.  The  results  indicate  that  fully-developed 
flow  yields  the  fastest  mixing  rates,  but  this  effect  is  relatively  small  at  the  highest 
pressures  (100  MPa  and  above)  where  the  LHF  approximation  is  likely  to  be  most  correct. 
For  both  flows,  however,  the  effect  of  pressure  is  quite  significant  with  the  liquid 
containing  region  becoming  smaller  at  high  pressures.  This  is  an  effect  of  increased 
entrainment  due  to  higher  ambient  gas  densities  at  high  pressures  (Faeth,  1987).  Since 
LHF  predictions  clearly  overestimate  the  length  of  the  liquid-containing  region  at  low 
pressures  (near  10  MPa)  the  reduction  of  the  length  of  the  liquid  containing  region  with 
increasing  pressure  would  be  even  larger,  due  to  combined  effects  of  better  atomization  and 
increased  entrainment.  However,  quantitative  estimates  of  spray  properties  at  high 
pressures  using  the  present  LHF  model  should  be  approached  with  caution,  pending 
evaluation  of  predictions  with  measurements. 


16 


RELATIVE  NUMBER  FREQUENCY  (%) 


Figure  9.  Measured  (present  study)  and  predicted  (DSF  and  SSF  methods) 
drop  size  contributions  for  combusting  LGP1845  sprays.  From 
Lee  and  Faeth  (1990)  . 
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Figure  11.  Predicted  (LHF  method)  mean  scalar  properties  along  axis  of  combusting  HAN-based 

monopropellant  sprays  for  fully-developed  flow  at  injector  exit.  From  Lee  et  al.  (1988) 
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Figure  12.  Predicted  (LHF  method)  mean  scalar  properties  along  axis  of  combusting  HAN-based 
mcnopropellant  sprays  for  slug  flow  at  injector  exit.  From  Lee  et  al.  (1988) . 


2.4  Conclusions 

Major  conclusions  of  the  present  investigation  can  be  summarized  as  follows: 

1 .  Critical  combustion  pressures  of  the  HAN-based  monopropellants  axe  unusually 
high,  250  MPa  with  an  uncertainty  of  50  percent.  This  implies  that  subcritical 
combustion  and  spray  processed  are  relevant  for  combusting  monopropellant 
sprays,  except  for  the  highest-pressure  portions  of  gun  cycles. 

2.  Estimates  of  liquid  surface  properties  and  critical  combustion  pressures  are 
unusually  sensitive  to  the  binary  interaction  parameter  of  TEAN  and  water. 
Measurements  of  this  and  other  properties  concerning  phase  equilibrium  are  needed 
to  reduce  uncertainties  of  critical  combustion  pressures. 

3.  Measurements  yielded  effective  drop  burning  rates  of  ca.  lOmm/s  for  drop 
diameters  of  300-1200  (im  and  pressures  of  0.2-7  MPa.  The  effective  drop 
burning  rate  involved  subsurface  gas  evolution  and  mechanical  removal  of  liquid  by 
microexplosions,  which  dominated  the  process  at  low  pressures  (below  2.1  MPa), 
and  conventional  hetereogeneous  combustion  from  the  liquid  surface,  which 
dominated  the  process  at  high  pressures. 

4.  Present  measurements  of  drop  burning  rates  at  pressures  of  0.2-7  MPa  were 
consistent  with  McBratney's  (1980,  1981)  strand  burning  rate  measurement  for 
jelled  propellants  at  10-100  MPa,  and  present  drop  combustion  rates  within 
combusting  sprays. 

5 .  Estimated  mean  drop  sizes  near  the  injector  exit  were  consistent  with  estimates  of 
aerodynamic  or  secondary  drop  breakup  theory  (Reitz  and  Bracco,  1982;  Ruff  and 
Faeth,  1990),  however,  this  approach  is  only  provisional  pending  study  of  effects 
of  liquid  turbulence  and  approach  of  the  flow  to  the  LHF  limit. 

6.  The  combusting  monopropellant  sprays  exhibited  strong  separated  flow  effects  for 
present  test  conditions,  with  the  length  of  the  liquid-containing  region  being  nearly 
the  same  for  injector  diameters  of  0.3  and  0.6mm  when  initial  SMD  and  injector 
velocities  were  nearly  the  same. 

7.  Both  separated  flow  models  were  in  reasonably  good  agreement  with 
measurements,  however,  the  DSF  approach  is  favored  since  it  requires  less 
extensive  computations. 

8.  The  LHF  approach  overestimated  the  rate  of  development  of  the  flow  for  present 
test  conditions  but  still  could  be  of  value  at  high  pressure  conditions  approaching 
and  exceeding  the  critical  combustion  pressure.  The  present  approach  is 
provisional,  however,  pending  proper  calibration  of  its  treatment  of  turbulent 
premixed  combustion. 
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Combustion  of  Han*Based  Liquid  Monopropellants  Near  the 
Thermodynamic  Critical  Point 

M.  E.  KOUNALAKIS  and  G.  M.  FAETH 
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The  high-preuure  combustion  properties  of  liquid  monopropellenu  involving  mixtures  of  hydroxyl  immonium 
nitrite  (HAN),  tnethanol  ammonium  nitrate  (TEAN).  and  water  are  considered  theoretically.  Liquid  surface 
properties  and  the  critical  combustion  pressure  (the  pressure  required  for  the  propellant  surface  to  exceed  its 
thermodynamic  critical  point)  were  found  allowing  for  real-gas  phenomena  and  the  presence  of  dissolved  combustion 
product  gases  in  the  liquid.  Critical  combustion  pressures  for  the  HAN-based  monopropellants  were  found  to  be 
unusually  high.  ca.  2500  atm.  with  an  estimated  uncertainty  of  50%.  Predictions  were  unusually  sensitive  to  the 
critical  temperature  of  TEAN  and  the  binary  interaction  parameter  between  TEAN  and  water:  both  must  be  known 
more  accurately  for  definitive  estimates  of  the  liquid-surface  properties  of  HAN-based  monopropellants  at  high 
pressures. 

NOMENCLATURE  Subscripts 


c  molar  concentration 

d  droplet  diameter 

D  effective  binary  diffusivity 

hj  partial  enthalpy  of  species  i 

k(j  binary  interaction  parameter 

h"  total  molar  flux 

p  pressure 

Pr  Prandtl  number 

R  gas  constant 

Re  Reynolds  number 

Sc  Schmidt  number 

T  temperature 

u  streamwise  velocity 

v  specific  volume 

x  distance 

Yj  mole  fraction  of  species  i 

a  thermal  diffusivity 

6  characteristic  convection  thickness 

€j  mole  flux  fraction  of  species  / 

X  mixture  thermal  conductivity 

w  acentric  factor 
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c  thermodynamic  critical  point 

f  flame  condition 

R  reactant 

o  far-upstream  condition 

INTRODUCTION 

The  combustion  of  liquid  monopropellant  sprays 
involves  pressure-atomized  injection  of  the  liquid 
into  a  high-pressure  and  high-temperature  envi¬ 
ronment  for  typical  applications,  e.g.,  gas  genera¬ 
tors  and  regenerative  liquid-propellant  guns.  After 
injection,  the  surfaces  of  ligaments,  drops,  etc., 
become  heated,  the  propellant  vaporizes  and, 
Anally,  reacts  to  combustion  products  in  the  gas 
phase.  At  sufficiently  high  pressures,  however, 
liquid  surfaces  approach  the  thermodynamic  criti¬ 
cal  point,  causing  transition  to  the  supercritical 
combustion  regime  at  the  critical  combustion 
pressure.  The  process  is  equivalent  to  the  single¬ 
phase  premixed  combustion  of  a  dense  gas  within 
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the  supercritical  combustion  regime,  and  liquid 
surfaces  are  no  longer  observed.  Clearly,  it  is 
important  to  know  critical  combustion  pressures 
and  liquid  surface  properties  to  define  the  pressure 
range  where  sprays  are  observed  and  to  find  spray 
properties  within  this  regime.  Analysis  was  under¬ 
taken  to  find  these  properties  for  hydroxyl  ammo¬ 
nium  nitrate-based  monopropellants  during  the 
present  investigation.  HAN-based  monopropet- 
lanis  are  of  interest  because  they  are  being 
considered  for  several  high-pressure  monopropel¬ 
lant  combustion  systems  [1], 

Near-critical  phenomena  of  high-pressure  liquid 
combustion  has  been  studied  earlier  in  this  Labora¬ 
tory  [2-7],  Various  liquid  monopropcllams,  such 
as  nitrate  esters  [2],  ethylene  oxide  [3],  and 
hydrazine  [4]  were  studied  in  strand  combustion 
configurations,  by  analyzing  gas  phase  transport 
using  the  thin  flame  approximation.  High-pressure 
phenomena,  such  as  real-gas  effects  and  the 
presence  of  gaseous  combustion  products  dis¬ 
solved  in  the  liquid  phase,  were  treated  using  the 
Redlich-Kwong  equation  of  state  [8].  Predictions 
generally  agreed  with  measured  liquid  surface 
temperatures,  within  uncertainties  anticipated  due 
to  limitations  concerning  the  high-pressure  ther¬ 
mophysical  and  transport  properties  of  combusting 
monopropellant  systems  [2-4],  Critical  combus¬ 
tion  pressures  were  in  the  range  100-200  atm  for 
the  nitrate  ester  monopropellants  (normal  propyl 
nitrate,  ethyl  nitrate,  and  propylene  glycol  dini¬ 
trate)  [2]  and  ethylene  oxide  (3),  which  is  typical 
of  other  liquid  combustion  processes  that  have 
been  studied  [5-7].  These  predictions  were  rela¬ 
tively  sensitive  to  values  of  transport  properties 
and  had  estimated  uncertainties  of  30%.  Unfortu¬ 
nately,  experimental  evaluation  of  predicted  criti¬ 
cal  combustion  pressures  was  not  possible,  due  to 
the  appearance  of  unstable  combustion  waves 
associated  with  the  loss  of  the  stabilizing  effect  of 
surface  tension  near  the  thermodynamic  critical 
point  [2,  3]  as  the  critical  combustion  pressure  was 
approached. 

Similar  studies  of  the  liquid  surface  properties 
of  HAN-based  monopropellants  have  not  been 
reported,  although  McBratney  [9,  10)  has  mea¬ 
sured  the  strand  burning  rates  of  some  typical 
blends.  A  liquid  surface  was  clearly  observed  at 
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290  atm,  while  jelled  samples  maintained  stable 
combustion  surfaces  up  to  600  atm.  These  findings 
suggest  significantly  higher  critical  combustion 
pressures  for  HAN-based  monopropellants  than 
for  other  monopropcllams  studied  s  far;  how¬ 
ever.  the  presence  of  a  jell  complicates  the 
interpretation  of  these  results.  The  results  also 
indicate  that  measurements  of  the  liquid  surface 
properties  of  HAN-based  monopropeilants  would 
be  very  difficult.  First,  the  burning  rates  of  HAN- 
based  monopropeilants  are  on  the  order  of  20  mm 
s'1  [9,  10),  which  is  an  order  of  magnitude  faster 
than  other  propellants  that  have  been  studied  [2, 
3],  Such  burning  rates  imply  characteristic  flame 
thicknesses,  a/u  for  each  phase,  and  characteristic 
flame  residence  times,  a/u1  for  each  phase,  on  the 
order  of  l  jim  and  l  fts  based  on  typical  transport 
properties.  Thus,  adequate  resolution  and  re¬ 
sponse  for  surface  temperature  measurements,  or 
even  for  distinguishing  a  liquid  surface  in  such  a 
thin  flame  zone,  is  unlikely.  Finally,  past  attempts 
to  measure  critical  combustion  pressures,  or  liquid 
surface  properties  near  the  thermodynamic  critical 
point,  have  not  been  successful  [2,  3]. 

In  view  of  the  experimental  problems,  theory 
was  used  to  find  liquid  surface  properties  and 
critical  combustion  pressures  during  the  present 
investigation.  Earlier  theoretical  methods  [2-4] 
were  extended  to  treat  the  specific  features  of 
some  representative  HAN-based  monopropeilants. 
Predictions  require  thermophy steal  and  transport 
properties  that  must  be  estimated  for  monopropel¬ 
lants;  therefore,  the  effects  of  these  uncertainties 
were  evaluated  using  sensitivity  analysis.  The 
study  was  limited  to  two  HAN-based  monopropei¬ 
lants,  LGP  1845  and  LPG  1846,  which  are 
mixtures  of  HAN,  triethanol  ammonium  nitrate, 
and  water.  These  choices  were  made  because 
McBratney  studied  LGP  1845  [10],  because  recent 
spray  studies  have  used  LGP  1846  [l],  and 
because  similar  formulations  are  candidates  for 
practical  applications. 

THEORY 

General  Description 

The  theory  is  an  extension  of  the  approach  used  in 
[3].  Phase  equilibrium  considerations  were  im- 


27 


LIQUID-MONOPROPELLANT  COMBUSTION 


181 


proved  by  considering  the  more  highly  developed 
modified  Redlich-Kwong  equation  of  state  due  to 
Soave  (11,  12  J,  as  well  as  the  basic  Redlich- 
Kwong  equation  of  state  [7]  used  by  this  labora¬ 
tory  in  earlier  work  [2,  3].  The  earlier  methods 
were  developed  for  the  analysis  of  monopropellant 
strand  combustion  (2,  3(;  the  relevance  of  this 
approach  for  estimating  liquid  surface  properties 
and  critical  combustion  pressures  in  sprays  was 
also  considered. 

As  in  past  work  [2-4],  the  combustion  process 
was  assumed  to  be  a  steady,  laminar  one-dimen¬ 
sional  deflagration  wave.  Typical  combusting 
spray  conditions  involve  pressures  greater  than  10 
MPa  and  drops  or  other  liquid  elements  having 
diameters  or  radii  of  curvature  greater  than  10  p.m, 
because  smaller  liquid  elements  approximate  the 
dynamics  of  the  gas  phase  in  any  event  (13). 
Under  these  conditions,  pressure  drops  across  the 
wave  are  less  than  100  Pa  even  for  the  high 
burning  rates  of  the  HAN-based  propellants; 
therefore  the  approximation  of  a  constant-pressure 
deflagration  wave  can  be  adopted  with  little  error. 
Based  on  the  burning  rate  measurements  of 
McBratney  [9,  10],  combustion  lifetimes  of  the 
liquid  elements  of  interest  are  greater  than  1  ms, 
which  is  much  greater  than  the  estimated  charac¬ 
teristic  residence  time  (1  ns)  in  the  wave;  there¬ 
fore,  the  steady  (quasisteady)  wave  approximation 
is  appropriate  as  well. 

Forced  convection  and  turbulence  in  sprays  can 
cause  departures  from  one-dimensional  waves 
near  liquid  surfaces.  Multidimensional  effects  due 
to  forced  convection  are  small  when  the  character¬ 
istic  convection  thickness  is  greater  than  the 
characteristic  flame  thickness  (14].  The  convec¬ 
tion  thickness  for  drops  can  be  estimated  with 
sufficient  accuracy  for  present  purposes  as  follows 
[13] 

5  =  4/(2  +  0.556  ReT(Pr  or  Sc)7).  (1) 

Maximum  drop  Reynolds  numbers  in  sprays  are 
on  the  order  of  10 2  (13);  therefore,  Eq.  1  implies 
that  characteristic  convection  thicknesses  are 
greater  than  the  1-^m  characteristic  flame  thick¬ 
ness  of  present  propellants  for  drop  diameters 
greater  than  10  /im.  Existing  information  on  the 


length  microscales  of  turbulence  in  spray  flames  is 
very  limited;  however,  one  estimate  for  the  fully 
developed  portion  of  sprays  having  a  typical 
Reynolds  number  of  104  suggests  that  the  micro¬ 
scales  are  on  the  order  of  10  times  the  distance 
from  the  injector  [13],  This  implies  that  micro¬ 
scales  are  greater  than  the  characteristic  flame 
thickness  for  distances  greater  than  1  mm  from  the 
injector  exit;  therefore,  turbulence  should  have 
little  effect  on  the  one-dimensionality  of  flame 
properties  for  most  portions  of  monopropellant 
spray  flames.  Based  on  these  considerations,  the 
assumption  of  one-dimensional  laminar  flow  ap¬ 
pears  to  be  relevant  for  spray  combustion  of  the 
present  monopropellants. 

Combustion  was  assumed  to  be  adiabatic  with 
thermodynamic  equilibrium  reached  at  the  down¬ 
stream  edge  of  the  reaction  zone,  which  are 
reasonable  assumptions  for  high- intensity  spray 
combustion  at  high  pressures.  The  properties  of 
the  combustion  products  were  computed  under 
these  assumptions  using  the  Gordon  and  McBride 
code  [15];  results  for  LGP  1845  and  LPG  1846  are 
summarized  in  Table  1  for  pressures  of  1  MPa,  10 
MPa,  and  100  MPa.  Both  reactant  mixtures  are 
stoichiometric  and  have  virtually  identical  com¬ 
bustion  product  compositions;  the  main  difference 
between  the  two  is  that  the  flame  temperature  of 
LGP  1845  is  roughly  100  K  higher  than  that  of 
LGP  1846.  The  combustion  products  are  nearly 
70%  water  vapor.  Due  to  the  relatively  low 
combustion  temperatures  and  high  pressures,  dis¬ 
sociation  effects  are  small,  and  combustion  prod¬ 
uct  concentrations  arc  nearly  independent  of  pres¬ 
sure. 

The  propellants  have  relatively  low  flame  tem¬ 
peratures;  therefore,  effects  of  radiation  can  be 
neglected.  For  example,  if  we  conservatively 
assume  that  the  emissivity  of  the  combustion 
products  is  unity  and  that  reabsorption  between  the 
flame  and  the  liquid  surface  is  negligible,  the 
burning  rate  measurements  of  McBratney  [9,  10] 
imply  that  radiation  contributes  less  than  4%  to  the 
enthalpy  rise  of  gasification  of  the  liquid. 

Precise  treatment  of  combustion  chemistry  is 
not  possible,  because  the  mechanisms  and  rate 
constants  are  not  known  for  the  high-pressure 
decomposition  of  HAN-based  monopropellants; 
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TABLE  I 

Combustion  Properties  of  HAN- Based  Mono  propel  lams* 


Propellant 

LGP  1843* 

LGP  t846r 

Pressure  (MPa) 

l 

10 

100 

1 

10 

100 

Reactant  density  (kg  m‘J) 

1432 

1454 

1476 

1430 

1432 

1454 

Combustion  product  properties 

Density  (kg  m-!) 

1.30 

12.9 

129 

1.36 

13.5 

135 

Temperature  (K) 

2133 

2152 

2163 

2027 

2039 

2045 

Composition  (%  by  volume)4 

Water  vapor 

68.8 

69.2 

69.3 

70.6 

70.8 

70.9 

Carbon  dioxide 

12.7 

12.9 

13.0 

12.2 

12.3 

12.3 

Nitrogen 

17.3 

17.4 

17.4 

16.5 

16.5 

16.6 

*  Adiabatic,  constant  pressure  combustion,  with  the  liquid  reactant  at  298.13  K. 
4  Reactant  composition  (%  by  mass):  HAN.  63.2;  TEAN,  20;  and  HA  16.8. 
f  Reactant  composition  <*  by  mass):  HAN.  60.8:  TEAN,  19.2;  and  HA  20.0. 
4  Major  species  only.  Minor  species  include:  CO.  H;,  NO.  OH.  and  <V 


therefore,  the  limiting  approximation  of  an  infi¬ 
nitely  thin  gas  phase  reaction  zone  located  at  some 
distance  from  the  liquid  surface  was  made,  as  in 
past  work  [2 -4).  This  large  activation  energy  limit 
is  favored  for  the  present  monopropeilants  due  to 
their  relatively  low  flame  temperatures. 

Consistent  with  the  thin  flame  approximation, 
liquid  phase  reactions  were  also  ignored.  This  is 
favored  by  the  very  short  residence  times  in  the 
high-temperature  portions  of  the  liquid  phase,  ca. 
100  ns.  Typical  unimolecular  decomposition  reac¬ 
tions,  which  are  at  least  representative  of  the  first 
stage  of  reaction  in  the  liquid,  have  characteristic 
reaction  times  two  or  more  orders  of  magnitude 
longer  than  the  liquid  residence  time  for  present 
conditions,  providing  some  justification  for  this 
assumption. 

Phase  equilibrium  was  assumed  at  the  liquid 
surface,  as  in  past  work  [2-4].  This  is  a  good 
approximation  for  drop  diffusion  flames  at  pres¬ 
sures  greater  than  0.1  MPa  [12];  however,  the 
assumption  was  reexamined  due  to  the  high 
burning  rates  of  the  HAN  propellants.  Following 
Dunn  and  Reay  [16],  the  burning  rate  measure¬ 
ments  of  McBratney  [9,  10]  imply  differences 
between  actual  and  saturated  vapor  pressures  at 
the  liquid  surface  less  than  0.3  MPa  due  to  finite 


rates  of  evaporation.  This  is  clearly  small  in 
comparison  to  the  combustion  pressures  of  inter¬ 
est,  which  were  greater  than  10  MPa.  The 
corresponding  temperature  jumps  at  the  surface, 
due  to  finite  evaporation  rates,  were  found  using 
the  present  phase  equilibrium  analysis.  The  result¬ 
ing  temperature  jumps  were  less  than  0.3  K,  which 
is  also  small  in  comparison  to  the  temperature  rise 
in  the  liquid  phase.  Thus,  phase  equilibrium  at  the 
surface  can  be  adopted  with  little  error. 

Gas  phase  transport  was  treated  similar  to 
earlier  analyses  of  liquid  combustion  processes  at 
high  pressures  [2-7],  Only  concentration  diffusion 
was  considered,  taking  the  binary  diffusivities  of 
all  species  to  be  equal,  and  neglecting  the  Dufour 
energy  flux.  The  assumption  of  an  effective  binary 
diffusion  coefficient  is  more  questionable  for 
blends  like  the  HAN-based  propellants  than  for 
blends  studied  in  earlier  work  where  only  a  single 
species  was  diffusing  and  the  approximation  was 
exact.  However,  the  present  approach  seems 
reasonable  in  view  of  the  other  approximations  of 
the  analysis;  the  effect  of  the  approximation  was 
examined  by  sensitivity  calculations.  The  specific 
treatment  of  thermophysical  and  transport  proper¬ 
ties  will  be  discussed  later.  Effects  of  compress¬ 
ibility  and  varying  Lewis  numbers  were  consid- 
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ered,  particularly  because  variations  in  the  Lewis 
number  are  known  to  influence  predictions  of 
critical  combustion  pressures  [2-7]. 

Monopropellants  typically  have  negligible  dis¬ 
solved  gas  concentrations  prior  to  injection;  there¬ 
fore,  carbon  dioxide  and  nitrogen  are  not  present 
in  the  bulk  liquid  and  have  negligible  mass  fluxes 
upstream  of  the  flame.  However,  the  propellant 
components  HAN,  TEAN,  and  water,  all  have 
finite  mass  fluxes  upstream  of  the  flame,  the  last 
being  both  a  propellant  component  and  a  combus¬ 
tion  product. 

Transport  Analysis 

Because  only  liquid  surface  properties  are  of 
interest  (not  burning  rates),  the  equations  to  be 
solved  are  identical  in  all  one-dimensional  coordi¬ 
nate  systems:  Cartesian  for  large  liquid  elements, 
cylindrical  for  ligaments,  and  spherical  for  drops. 
A  Cartesian  formulation,  sketched  in  Fig.  1,  is 
used  to  simplify  the  notation.  The  analysis  is 
confined  to  the  region  between  the  flame  and 
liquid  surface,  corresponding  to  the  outer  region 
of  asymptotic  analysis  of  monopropellant  flames  at 
the  high  activation  energy  limit.  Effects  of  multi- 
component  diffusion  and  the  presence  of  dissolved 
gases  cause  reactant  (HAN  and  TEAN)  concentra¬ 
tions  to  vary  in  the  liquid  phase.  There  is  a 
concentration  jump  at  the  liquid  surface,  dictated 
by  phase  equilibrium  requirements.  The  corres¬ 
ponding  temperature  jump  is  small,  however,  as 
noted  earlier.  Properties  at  the  flame  surface  are 


Fig.  1.  Sketch  of  the  liquid  monopropellant  combustion 
process. 


given  by  the  adiabatic  flame  calculations  discussed 
earlier,  see  Table  1 . 

Gas  phase  transport  was  formulated  in  terms  of 
molal  quantities  to  simplify  matching  with  the 
phase  equilibrium  analysts.  Under  present  as¬ 
sumptions,  the  equations  governing  conservation 
of  energy  and  species  between  the  liquid  surface 
and  the  flame  sheet  are  as  follows  [2-4] 

"  dr 

d  *  Y  A,o)  =  X  —  (2) 

,T!  dx 

d  Yi 

d ’  ( Yi-ti) ~ cD  — ■  ,i*l,  N.  (3) 

dx 

Initial  conditions  for  Eqs.  2  and  3  can  be  specified 
in  terms  of  known  flame  properties  as  follows 

r-r„  yj-r*/-!,  (4) 


Denote  the  species  present  in  the  original  mono¬ 
propellant  as  /  =  1 ,  •  *  • ,  R  <  N.  Then, 
eliminating  spatial  derivatives  among  Eq.  3,  inte¬ 
grating,  and  applying  the  concentration  boundary 
conditions  of  Eq.  4  yields 


(*>-*/> 


(5) 


where  i  *  j,  and  j  can  be  any  species  y  =  1, 

R. 

Eliminating  spatial  derivatives  between  Eq.  2 
and  Eq.  3  for  /  =  j  yields 


dYj  ( \/cD)(tj-Yj ) 

—  =  ~~  (o) 

dr 

i- 1 

where  j  is  selected  in  the  same  manner  as  Eq.  5. 
The  initial  condition  for  Eq.  6  is 

T=Th  Yj=Yj(.  (7) 

As  noted  earlier,  the  effects  of  pressure  on 
flame  properties  were  small;  therefore,  a  single 
flame  condition  was  used  for  each  propellant. 
These  properties,  as  well  as  the  e,  that  are  fixed  by 
the  propellant  composition,  are  summarized  in 
Table  2  for  each  propellant. 
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TABLE  2 


Mole  Flux  Fraction*  and  Flame  Properties  for  Calculations 


Fuel 

LGP  1843 

LGP  1846 

Mole  Flux  Fractions  (%) 

HAN 

39  1 

34.5 

TEAN 

5.6 

4.9 

H,Q 

55.3 

60,6 

Flame  Properties: 

Temperature  (K) 

2146 

2040 

Composition  (%  by  volume): 

HjO 

69 

71 

CO, 

13 

12 

N, 

IS 

17 

Given  Eq.  5  and  mixture  properties  as  a 
function  of  temperature  and  pressure,  Eq.  6  can  be 
integrated  from  the  flame  toward  the  liquid  sur¬ 
face,  providing  the  transport  locus  of  the  variation 
of  species  mole  fractions  with  temperature.  Liquid 
surface  conditions  are  reached  when  the  properties 
along  the  transport  locus  reach  a  dew  point  of  the 
phase  equilibrium  analysis.  Dew-point  computa¬ 
tions  are  considered  next. 

Phase  Equilibrium  Analysis 

Two  methods  were  used  to  compute  phase  equilib¬ 
rium  at  the  liquid  surface:  the  Redlich-Kwong 
equation  of  state  with  mixing  rules  developed  by 
Prausnitz  and  Chueh  [8],  which  was  used  during 
earlier  studies  of  high-pressure  liquid  combustion 
[2-7];  and  the  more  highly  developed  modified 
Redlich-Kwong  equation  of  state  due  to  Soave 
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[11,  12].  The  use  of  these  equations  of  state  for 
polar  compounds  like  HAN,  TEAN,  and  waier  is 
tentative,  for  lack  of  a  viable  alternative.  Prausnitz 
et  al.  [17]  found,  however,  that  polar  corrections 
of  the  equation  of  state  are  small  when  reduced 
temperatures  are  greater  than  0.95;  this  condition 
was  generally  satisfied  for  the  most  polar  com¬ 
pounds  during  present  computations. 

Phase  equilibrium  at  the  liquid  surface  requires 
that  the  temperature,  pressure,  and  fugacity  of 
each  species  is  the  same  in  both  phases.  Given  the 
pressure  and  the  gas  phase  composition  along  the 
transport  locus,  the  corresponding  liquid  phase 
composition  and  equilibrium  temperature  can  be 
computed  from  the  equation  of  state.  This  corres¬ 
ponds  to  a  dew-point  calculation,  using  the  com¬ 
puter  codes  found  in  [  1 1]  and  [7]  for  the  Soave  and 
Relich-Kwong  equations  of  state.  Liquid  surface 
conditions  are  reached  when  the  transport  and 
phase  equilibrium  temperatures  are  the  same. 


Thermophysical  Properties 

The  equations  of  state  require  critical  properties 
and  acentric  factors,  while  computations  of  flame 
properties  require  enthalpies  of  formation.  The 
values  of  these  properties  used  during  the  calcula¬ 
tions  are  summarized  in  Table  3.  The  properties  of 
H^O,  COj,  and  Nj  listed  in  this  table  are  drawn 
from  standard  references  [18.  19],  and  present  no 
problems. 

HAN  and  TEAN  decompose  at  high  tempera¬ 
tures;  therefore,  their  thermophysical  properties 
must  be  estimated.  The  Lyderson  method  [20,  21] 


TABLE  3 

Thermophysical  Properties 


Species 

HAN 

TEAN 

H,0 

CO, 

N, 

Critical  Properties 

Temperature  (K) 

763 

1294 

647.3 

304.2 

126.2 

Pressure  (atm.) 

76 

36 

217.6 

72.8 

33.5 

Volume  (craVgmol) 

196 

528 

56 

94 

89.5 

Acentric  factor 

0.68 

1.41 

0.344 

0.225 

0.040 

Enthalpy  of  formation 

-95.3 

-  185.3 

-57.8 

-94.0 

0.0 

at  298.15  K(kcal  gmol-') 

31 


LIQUID-MONOPROPELLANT  combustion 

was  used  to  estimate  the  critical  pressure,  the 
critical  volume,  and  the  ratio  of  the  normal  boiling 
and  critical  temperatures.  Acentric  factors  were 
computed  from  their  definition,  using  the  Lee  and 
Kesier  generalized  vapor-pressure  relationship 
[22].  The  critical  temperature  was  then  computed 
from  the  Lee  and  Kesier  expression  [22]  for  the 
compressibility  factor,  applied  at  the  critical  point, 
as  follows 

Tt  =  (pcuc//?)/(  0,291  -  O.OSOu)  (8) 

These  methods  are  routinely  used  for  hydrocar¬ 
bons,  but  they  have  not  been  tested  for  compounds 
like  HAN  and  TEAN.  Thus,  the  approach  was 
evaluated  by  applying  it  to  a  variety  of  organic  and 
inorganic  compounds  whose  critical  properties 
and  acentric  factors  were  known:  NHj,  C02, 
H20,  SOj,  N02,  n-hexanol,  benzyl  alcohol,  n, 
/t-dimethylaniline,  n-hexadecane,  n-eicosane,  and 
n-butanol.  The  average  errors  were  as  follows: 
critical  temperature,  -5%;  critical  pressure, 
-2.5%;  critical  volume,  -9.6%;  and  the  acen¬ 
tric  factor,  0.6%.  Based  on  these  findings,  the 
approach  appears  to  be  reasonable  for  present 
purposes. 

The  enthalpies  of  formation  of  HAN  and  TEAN 
listed  in  Table  3  were  found  using  the  Vertna  and 
Doraiswamy  group  contribution  method  [23], 
Uncertainties  in  these  properties  are  not  very 
important,  however,  because  flame  properties  are 
dominated  by  the  enthalpies  of  formation  of  H20, 
C02,  and  N2,  which  are  well  known  [18,  19]. 

(deal-gas  enthalphy  changes  with  temperature 
were  found  from  standard  sources  for  H20,  CC^, 
and  N2  [18,  19].  The  group  contribution  specific 
heat  correlation  of  Rihani  and  Doraiswamy  [24] 
was  used  to  compute  ideal-gas  enthalpy  changes 
for  HAN  and  TEAN.  Enthalpy  deviations  and  c 
were  found  directly  from  the  equations  of  state. 

Binary  interaction  parameters  k,,  are  used  in  the 
equations  of  state.  Present  values  of  the  kfJ  are 
summarized  in  Table  4.  These  estimates  were 
found  from  (l  l,  12]  for  the  Soave  equation  of 
state,  and  from  [5,  8]  for  the  Redlich-Kwong 
equation  of  state.  Measurements  of  kij  are  not 
available  for  binary  pairs  involving  HAN  and 
TEAN;  therefore,  the  values  for  the  hydrocarbon 
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TABLE  4 


Binary  Interaction  Parameters  (*,,)“ 


Substance 

H.O 

CO, 

N, 

Soave  equation  of  state 

HAN 

0.381 

0  136 

0.140 

TEAN 

0.368 

0.138 

0.164 

H.O 

0.000 

0.102 

0.140 

CO; 

0.102 

0.000 

-0.022 

Redlich-Kwong  equation  of  state 

HAN 

0.383 

0.283 

0.333 

TEAN 

0.790 

0.650 

0,743 

*  k„  =  0:  k,,  ■  k,,\  other  interaction  parameters  not  shown 
were  taken  to  be  rero.  The  system*  were  modeled  by 
homo  morphs. 


homomorphs  of  these  species  were  used,  as  in  past 
work  [2-7], 

Transport  Properties 

Mixture  thermal  conductivities  and  effective  bi¬ 
nary  diffusivities  are  needed  to  integrate  Eq.  6. 
Correlations  of  mixture  thermal  conductivities 
require  the  viscosities  of  all  species  as  well.  Ideal- 
gas  thermal  conductivities  and  viscosities  were 
taken  from  Svehla  [25]  for  H20,  CO*,  and  N2, 
while  the  Eucken  model  [25]  and  the  Reichenberg 
method  [26]  were  used  to  estimate  these  properties 
for  HAN  and  TEAN.  Ideal  gas  mixture  thermal 
conductivities  were  then  found  from  the  Wassil- 
jewa  equation  [27],  following  Lindsay  and  Brom¬ 
ley  [26].  The  effect  of  compressibility  on  the 
mixture's  thermal  conductivity  was  determined  by 
the  Stiel  and  Thodos  method  [29],  using  the 
Prausnitz  and  Gunn  rules  for  mixture  pseudocriti- 
cal  properties  [18]. 

Low-pressure  binary  diffusivities  were  found 
for  ail  binary  pairs  in  the  system  using  the  method 
of  Fuller,  Schettler  andGiddings  [30, 31]  based  on 
the  average  of  these  values  at  each  temperature. 
The  effect  of  compressibility  on  the  effective 
binary  diffusivity  was  computed  following  Daw¬ 
son  et  al.  [32],  using  the  Prausnitz  and  Gunn 
modified  rules  for  mixture  pseudocritical  proper¬ 
ties  [18],  as  before. 
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RESULTS  and  discussion 
Critical  Combustion  Conditions 

Results  for  LGP  1845  and  LPG  1846  were  very 
similar,  in  view  of  property  uncertainties;  there¬ 
fore,  findings  for  LGP  1845  will  be  talcen  as 
representative  in  the  following. 

The  nature  of  the  solution  for  liquid  surface 
properties,  in  the  vicinity  of  the  thermodynamic 
critical  point,  can  be  seen  by  reference  to  Fig.  2. 
Transport  and  phase  equilibrium  loci  are  plotted  as 
functions  of  the  gas  phase  mole  fraction  of  HAN 
and  the  temperature  at  various  pressures.  The 
phase  equilibrium  locus  was  found  using  the  Soave 
equation  of  state  with  all  ky  -  0;  however,  the 
properties  of  the  plots  are  similar  to  findings  for 
the  range  of  ky  considered  during  the  investiga¬ 
tion.  Both  loci  satisfy  Eq.  5,  which  relates  the 
mole  fractions  of  TEAN,  HjO,  COi ,  and  to  the 
mole  fraction  of  HAN  in  the  gas  phase.  In 
addition,  the  transport  locus  satisfies  the  conserva¬ 
tion  of  energy  requirements  of  Eq.  6,  while  the 
phase  equilibrium  locus  satisfies  the  requirements 
of  thermodynamic  equilibrium  at  the  liquid  sur¬ 
face.  The  intersection  of  these  two  loci,  at  a  given 
pressure,  is  a  solution  for  the  gas  phase  properties 
at  the  liquid  surface,  where  Eqs.  5  and  6  and  the 
requirements  for  phase  equilibrium  are  all  satis¬ 
fied. 

At  the  lowest  pressure  illustrated  in  Fig,  2,  500 
atm,  the  transport  and  phase  equilibrium  loci 
intersect,  yielding  a  liquid-surface  solution  point. 
This  behavior  persists  for  all  pressures  below  640 
atm.  However,  the  transport  locus  becomes  tan¬ 
gent  to  the  phase  equilibrium  locus  at  640  atm;  this 
is  the  highest  pressure  where  a  liquid  surface  is 
present  for  the  conditions  illustrated  in  Fig.  2,  and 
it  corresponds  to  the  critical  combustion  pressure. 
Pressures  greater  than  640  atm  are  within  the 
supercritical  combustion  regime  where  the  trans¬ 
port  and  phase  equilibrium  loci  no  longer  inter¬ 
sect,  see  the  results  illustrated  in  Fig.  2  for  650 
atm.  In  this  regime,  flow  properties  simply  vary 
along  the  transport  locus  until  the  initial  conditions 
in  the  bulk  liquid  are  reached,  and  a  liquid  surface 
is  never  observed. 


Liquid  Surface  Properties 

The  properties  of  the  liquid  surface  for  the 
subcritical  regime,  using  the  same  solution  param¬ 
eters  as  Fig.  2  (the  Soave  equation  of  state  with  all 
kij  =  0),  are  illustrated  in  Fig.  3.  The  surface 
temperature  and  species  mole  fractions  in  the  gas 
and  liquid  phases  are  plotted  as  a  function  of 
pressure  up  to  the  critical  combustion  pressure. 

Gas  phase  concentrations  are  relatively  indepen¬ 
dent  of  pressure  in  Fig.  3;  however,  liquid  phase 
concentrations  vary  substantially  due  to  the  pres¬ 
ence  of  dissolved  gases  in  the  liquid  at  higher 
pressures.  Other  liquid  combustion  systems  be¬ 
have  in  a  similar  manner  [2-7],  At  low  pressures, 
concentrations  of  dissolved  combustion  product 
gases  (COj  and  Nj)  are  quite  small,  and  the  more 


Fig.  2.  Transport-  and  phase -equilibrium  loci  for  subcritical, 
critical,  and  supercritical  conditions;  Soave  equation  of  state, 
k„  =■  0. 
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Fig.  3.  Properties  of  the  liquid  surface  as  a  function  of 
pressure:  Soave  equation  of  state,  k,,  »  0. 

volatile  propellant  components  (HAN  and  HjO) 
are  depleted  in  the  liquid  phase,  leaving  a  high 
concentration  (a  mole  fraction  near  0.9)  of  the 
relatively  nonvolatile  TEAN.  As  the  pressure  is 
increased,  however,  concentrations  of  the  more 
volatile  components  increase  in  the  liquid,  with 
water  having  the  highest  concentrations  in  the 
liquid  for  pressures  greater  than  400  atm.  Concen¬ 
trations  of  the  gas  and  liquid  are  not  the  same  at  the 
critical  combustion  condition,  and  neither  are  the 
enthalpies;  therefore,  the  critical  combustion  con¬ 
dition  does  not  correspond  to  a  thermodynamic 
critical  point  in  this  case.  Thus,  while  the  critical 
combustion  condition  often  is  a  thermodynamic 
critical  point  [2-7],  such  behavior  is  not  always 
observed. 

At  low  pressures,  the  surface  temperatures 
plotted  in  Fig.  3  increase  with  increasing  pressure, 
which  is  typical  of  the  behavior  of  other  monopro- 
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pellants  that  have  been  studied  [2-4].  However, 
the  surface  temperature  reaches  a  maximum,  and 
then  declines  as  the  critical  combustion  pressure  is 
approached.  Similar  behavior  has  been  observed 
for  ethylene  oxide  [3],  The  effect  is  due  to  the 
greater  concentrations  of  more  volatile  species, 
which  have  relatively  low  critical  temperatures,  in 
the  liquid  phase  at  high  pressures.  Due  to  the 
relatively  low  volatility  of  HAN  and  TEAN, 
however,  liquid  surface  temperatures  are  unusu¬ 
ally  high;  e.g.,  peak  temperatures  approach  1000 
K  in  Fig.  3,  as  opposed  to  500  K  for  the  nitrate 
esters  and  ethylene  oxide  [2,  3]. 

The  best  estimate  of  liquid  surface  properties 
for  LGP  1845,  using  the  Soave  equation  of  state 
and  the  standard  ky  values,  is  illustrated  in  Fig.  4. 


Fig.  4.  Properties  of  the  liquid  surface  as  a  function  of 
pressure:  Soave  equation  of  state,  standard  k,r 
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The  change  in  the  ky  causes  a  substantial  increase 
'  in  the  predicted  critical  combustion  pressure,  from 
640  atm  for  ky  =  0,  to  2600  atm  for  the  standard 
ky.  Such  strong  effects  of  the  ky  on  liquid  surface 
properties  were  not  observed  during  earlier  studies 
of  liquid  monopropellants  (1-3J.  For  the  results 
illustrated  in  Fig.  4,  the  critical  combustion 
condition  coincides  with  a  thermodynamic  critical 
point;  this  is  indicated  by  the  equality  of  gas  and 
liquid  compositions  at  the  critical  combustion 
pressure.  Other  differences  when  the  standard  ku 
are  used  included  a  greater  variation  of  gas  phase 
compositions  and  higher  concentrations  of  HAN  in 
the  liquid  phase,  cf..  Figs.  3  and  4.  However,  the 
range  of  liquid  surface  temperatures  and  the 
presence  of  a  peak  temperature  at  low  pressures 
are  similar  in  Figs.  3  and  4. 

The  Soave  equation  of  state  was  felt  to  be  the 
most  reliable,  because  the  ky  values  for  the 
Redlich-Kwong  equation  of  state  were  unusually 
high,  and  solutions  with  this  approach  were 
numerically  very  stiff.  Nevertheless,  results  were 
obtained  using  the  Redlich-Kwong  equation  of 
state,  because  it  has  been  used  for  monopropel¬ 
lants  in  the  past  [2,  3].  Liquid  surface  properties 
computed  in  this  manner  for  LGP  1845  are 
illustrated  in  Fig.  5  using  standard  values  of  the 
ky.  The  main  difference  between  the  results  using 
the  Soave  and  Redlich-Kwong  equations  of  state, 
cf.  Figs.  3  and  4,  is  that  the  latter  yields  a  much 
higher  critical  combustion  pressure,  higher  than 
the4000-atm  pressure  range  considered  in  Fig.  5. 
Thus,  the  results  in  Fig.  5  are  similar  to  the  low- 
pressure  region  of  Fig.  4:  the  gas  phase  composi¬ 
tions  are  relatively  independent  of  pressure;  the 
liquid  phase  concentrations  of  the  more  volatile 
components  (C02 ,  N2,  and  H20)  increase  with 
pressure;  and  the  surface  temperature  increases 
with  pressure.  Predicted  liquid  surface  tempera¬ 
tures  become  quite  high  as  the  pressure  is  in¬ 
creased,  which,  if  true,  raises  concerns  about 
potentially  significant  effects  of  reaction  near  the 
liquid  surface. 

Sensitivity  Analysis 

The  equation  of  state  and  the  values  of  the  ky 
clearly  have  a  significant  effect  on  predictions  of 
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Fig.  5.  Properties  of  the  liquid  surface  is  a  function  of 
pressure;  Redlich-Kwong  equation  of  state,  standard  k,,. 

liquid  surface  conditions  and  the  critical  combus¬ 
tion  pressure.  Sensitivity  analysis  was  undertaken 
to  more  precisely  identify  parameters  that  have  a 
strong  influence  on  predictions.  The  baseline  for 
the  sensitivity  calculations  was  the  Soave  equation 
of  state  with  all  ky  =  0.  An  exhaustive  search 
procedure  was  used,  where  the  effect  of  each 
parameter  was  evaluated  independently. 

Table  5  is  a  summary  of  the  results  of  the 
sensitivity  analysis,  considering  only  the  most 
sensitive  parameters.  The  percent  change  in  the 
critical  combustion  pressure,  the  critical  combus¬ 
tion  temperature,  the  pressure  where  the  surface 
temperature  is  maximum,  and  the  maximum  sur¬ 
face  temperature  are  given  as  a  function  of  the 
percent  change  of  the  parameters  (or  the  change  in 
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TABLES 


Sensitivity  Analysis  (LOP  1843) 


Output 

variable 

Critical 

combustion 

pressure 

(atm) 

Critical 

combustion 

temperature 

(K) 

Pressure  at 
maximum 
temperature 
(atm) 

Maximum 

temperature 

(K) 

Standard  Value4 

640 

840 

200 

1030 

Parameter* 

Variation 

Percent  change 

Thermal 

+  23* 

-  18 

10 

0 

-0 

conductivity 

-23* 

6 

2 

0 

-0 

DiflusioQ 

+  23* 

6 

2 

0 

-0 

coefficient 

-23* 

-  18 

10 

-14 

-0 

Critical  temperature 

+  10* 

44 

5 

0 

11 

of  TEAN 

-10* 

-34 

-4 

-23 

-18 

Acentric  factor 

+  20* 

13 

2 

0 

2 

of  TEAN 

-20* 

-  18 

2 

0 

2 

V 

0.23 

80 

-11 

23 

0 

TEAN-HjO 

4  Sands ni  property  values  except  all  ky  =*  0.  Soave  equation  of  state. 
9  Results  for  the  most  sensitive  parameters  are  listed. 

‘  All  other  kv  »  0. 


value  of  the  ky).  The  most  sensitive  parameters 
include  the  mixture  thermal  conductivity,  the 
effective  diffusion  coefficient,  the  critical  temper¬ 
ature  and  acentric  factor  of  TEAN,  and  the  binary 
interaction  parameter  between  TEAN  and  water. 
For  continuous  variables,  such  as  the  mixture 
thermal  conductivity  or  the  effective  diffusion 
coefficient,  the  parameter  variations  involved 
changing  the  parameter  from  the  standard  value, 
by  the  percentage  shown,  at  every  point  in  the 
calculation.  Except  for  the  ky,  which  will  be 
considered  subsequently,  the  variations  are  repre¬ 
sentative  of  anticipated  uncertainties  [18,  25], 
similar  to  past  work  [2,  3]. 

Uncertainties  in  the  critical  combustion  pressure 
due  to  uncertainties  in  the  thermal  conductivity, 
the  diffusion  coefficient,  and  the  acentric  factor  of 
TEAN,  are  comparable  to  the  uncertainties  in  the 
parameter  itself,  which  is  similar  to  findings  for 
other  monopropellants  (2,  3],  The  effect  of 
changes  in  these  parameters  on  other  output 


variables,  however,  is  relatively  small.  On  the 
other  hand,  except  for  the  critical  combustion 
temperature,  all  output  variables  exhibit  unusually 
strong  sensitivity  to  uncertainties  in  the  critical 
temperature  of  TEAN.  This  effect  appears  to  be 
related  to  the  vapor  pressure  of  TEAN,  with  the 
reduced  volatilities  of  TEAN  at  higher  critical 
temperatures  driving  the  system  toward  higher 
critical  combustion  pressures.  The  importance  of 
TEAN  for  critical  combustion  phenomena,  in  spite 
of  its  relatively  low  concentration  in  the  propel¬ 
lant,  is  probably  due  to  its  low  volatility,  which 
causes  it  to  have  relatively  high  concentrations  at 
the  liquid  surface,  see  Fig.  3. 

The  high  concentrations  of  both  TEAN  and 
water  at  the  liquid  surface  are  also  responsible  for 
the  strong  sensitivity  of  the  critical  combustion 
pressure  to  variations  in  the  value  of  the  binary 
interaction  parameter  between  TEAN  and  water. 
Higher  values  of  this  binary  interaction  parameter 
tend  to  increase  concentrations  of  TEAN  in  the 
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liquid  at  the  surface,  acting  in  a  manner  similar  to 
the  critical  temperature  of  TEAN.  Comparable 
changes  of  ail  other  had  a  small  effect, 
influencing  output  parameters  less  than  10%, 
similar  to  past  findings  concerning  effects  of 
uncertainties  in  the  ktj  [2,  3].  Thus,  the  large 
increase  in  critical  combustion  pressure  between 
Fig.  3  for  all  k^  =  0,  and  Fig.  4  for  the  standard 
kjj,  is  largely  due  to  change  of  the  binary 
interaction  parameter  between  TEAN  and  water. 

Findings  were  similar  using  the  Redlich-Kwong 
equation  of  state.  Predictions  were  most  sensitive 
to  the  binary  interaction  parameter  between  TEAN 
and  water,  and  to  a  lesser  degree  to  the  critical 
temperature  of  TEAN. 

The  direct  effect  of  the  binary  interaction 
parameter  between  TEAN  and  water  on  critical 
combustion  pressure  predictions  is  illustrated  in 
Fig.  6.  Results  are  illustrated  for  both  the  Soave 
and  Redlich-Kwong  equations  of  state.  All  param¬ 
eters  are  at  their  standard  values,  except  for  the 
TEAN-water  binary  interaction  parameter  which 
varies  as  shown.  The  solution  at  the  standard  value 
of  the  interaction  parameter  is  marked  for  the 
Soave  equation  of  state:  the  standard  critical 
combustion  condition  for  the  Redlich-Kwong 
equation  of  state  is  outside  the  pressure  range 
illustrated  in  Fig.  6.  The  critical  combustion 
pressure  is  relatively  insensitive  to  the  binary 
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Fig.  6.  Effect  of  the  TEAN-w»ter  binary  interaction  perime¬ 
ter  on  the  predicted  critic*!  combustion  pressure. 


interaction  parameter  between  TEAN  and  water  at 
low  values  of  this  parameter,  but  becomes  very 
sensitive  in  the  range  of  the  standard  values.  For 
the  Soave  equation  of  state,  a  25  %  variation  of  the 
TEAN-water  interaction  parameter  about  the 
standard  value  would  result  in  a  50%  variation  in 
the  critical  combustion  pressure.  This  sensitivity  is 
smaller  than  the  sensitivity  to  the  critical  tempera¬ 
ture  of  TEAN  shown  in  Table  5;  however, 
uncertainties  in  the  binary  interaction  parameter 
between  TEAN  and  water  are  larger  than  those  for 
the  critical  temperature  of  TEAN.  Thus,  this 
interaction  parameter  will  have  to  be  known  more 
accurately  for  definitive  estimates  of  liquid  surface 
temperatures  and  critical  combustion  pressures  of 
LGP  1845  and  LPG  1846. 

Based  on  these  considerations,  the  best  estimate 
of  the  critical  combustion  pressure  for  LGP  1845 
and  LPG  1846  is  2500  atm,  with  an  uncertainty  of 
50%.  This  follows  from  predictions  using  the 
Soave  equation  of  state  with  an  uncertainty  of  10% 
for  the  critical  temperature  of  TEAN,  which  is 
reasonable  in  view  of  the  assessment  of  estimates 
of  this  property  discussed  earlier,  and  uncertain¬ 
ties  of  all  other  properties  of  25  % ,  which  is  typical 
of  past  practice  [2-7].  A  low  estimate  of  the 
critical  combustion  pressure  would  be  600  atm, 
based  on  the  asymptotic  of  the  critical  combustion 
pressure  plot  for  the  Soave  equation  of  state 
illustrated  in  Fig.  6.  The  low  estimate  is  in  the 
region  where  McBratney  [9,  10]  observed  the 
onset  of  unsteady  combustion  waves,  which  often 
been  associated  with  near-critical  phenomena; 
however,  the  use  of  jelled  propellants  during  these 
tests,  and  the  absence  of  direct  evidence  concern¬ 
ing  the  onset  of  supercritical  combustion,  raises 
questions  concerning  the  interpretation  of  these 
findings.  The  Redlich-Kwong  equation  of  state 
yields  critical  combustion  pressures  greater  than 
4000  atm;  however,  estimated  binary  interaction 
parameters  for  this  approach  are  unusually  high. 
Therefore,  this  finding  is  felt  to  be  less  reliable 
than  results  using  the  Soave  equation  of  state. 

In  spite  of  the  uncertainties,  it  seems  clear  that 
the  critical  combustion  pressures  of  the  HAN- 
based  monopropellants  LGP  1845  and  LPG  1846 
are  unusually  high,  roughly  an  order  of  magnitude 
higher  than  other  liquid  monopropellant  and  bipro- 


37 


UQUID-MONOPROPELLANT  COMBUSTION 

pellant  combustion  processes  that  have  been  stud¬ 
ied  (2-7).  Present  uncertainties  in  estimates  of 
liquid  surface  properties  and  critical  combustion 
pressures  could  be  reduced  if  liquid  surface 
temperatures  were  measured;  however,  this  is  a 
very  challenging  experimental  problem  due  to  the 
unusually  high  burning  rates  of  these  monopropel¬ 
lants. 

CONCLUSIONS 

The  main  conclusions  of  the  present  investigation 
are  as  follows; 

1 .  Present  findings  suggest  unusually  high  critical 
combustion  pressures  for  typical  HAN-based 
monopropellants  (LGP  1845  and  LPG  1846), 
on  the  order  of  2500  atm  with  an  uncertainty  of 
50%.  This  suggests  that  subcritical  combustion 
and  spray  processes  are  relevant  to  the  combus¬ 
tion  of  these  monopropellants  for  most  applica¬ 
tions. 

2.  Estimates  of  the  liquid  surface  properties  of  the 
present  HAN-based  monopropellants  are  un¬ 
usually  sensitive  to  the  binary  interaction  pa¬ 
rameter  between  TEAN  and  water,  and  to  a 
lesser  degree,  to  the  critical  temperature  of 
TEAN. 

3.  Depending  on  the  thermophysical  and  transport 
properties  of  the  system,  conditions  at  the 
liquid  surface  at  the  highest  pressure  where  a 
liquid  surface  is  observed  (the  critical  combus¬ 
tion  condition)  may  or  may  not  correspond  to  a 
thermodynamic  critical  point. 

4.  The  present  HAN-based  monopropellants  ex¬ 
hibited  relatively  high  liquid  surface  tempera¬ 
tures,  800-1000  K,  in  the  pressure  range  of 
interest  (pressures  greater  than  100  atm).  This 
reduces  polar  liquid  effects,  which  complicate 
the  thermodynamics  at  lower  temperatures,  but 
increases  the  potential  for  significant  liquid 
phase  reaction  near  the  surface,  although  the 
short  residence  times  found  under  these  condi¬ 
tions  (ca.  100  ns)  mitigate  the  potential  reaction 
effects  to  some  degree). 

It  should  be  noted  that  McBratney  [9,  10]  has 
observed  unstable  combustion,  which  is  often 
associated  with  critical  combustion  conditions  [2, 
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3],  at  pressures  of  600-1000  atm;  therefore, 
present  estimates  may  be  high.  New  measure¬ 
ments  are  needed  to  reduce  the  uncertainties  of 
present  findings  and  to  definitively  establish  the 
range  of  pressures  where  spray  phenomena  must 
be  considered  for  these  monopropellants. 
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( Recent  d  May  l.  1987;  in  final  form  August  J.  1987) 

Abstract— A  simplified  analysis  of  monopropeilam  spray  combuslion  was  developed,  based  on  the 
locally-homogeneous-flow  approximation  of  multiphase  (low  theory  and  the  thin-fiame  approximation 
of  turbulent  premixed  flame  theory.  The  performance  of  the  analysis  was  evaluated  using  shadowgraphs 
of  spray  flames  lor  a  hydroxyl  ammonium  nitrate  iHANbbased  liquid  monopropeilam  al  ambient 
pressures  of  b-tj  MPa.  Predictions  showed  that  these  spray  flames  are  very  sensitive  to  the  degree  of 
flow  development  at  the  injector  exit,  with  fully-developed  turbulent  flows  requiring  significantly 
smaller  combustion  volumes  than  slug  flows  having  low  iniiiaj  turbulence  levels.  There  was  encouraging 
agreement  between  predictions  and  measurements;  however,  uncertainties  concerning  injector  exit 
conditions  for  the  experiments  precluded  definitive  assessment  of  the  analysis 


INTRODUCTION 

Combusting  monopropellant  sprays  have  applications  for  throttlable  thrustors, 
underwater  propulsion  systems,  and  regenerative  liquid-propellant  guns. 
Monopropellant  spray  flames  are  also  an  important  fundamental  problem  of 
combustion  science,  as  the  premixed  counterpart  of  the  spray  diffusion  flame. 
Motivated  by  these  considerations,  the  structure  of  combusting  monopropeilam 
sprays  was  examined  during  the  present  investigation.  The  main  objective  was  to 
develop  a  simplified  analysis  of  the  process,  and  to  evaluate  predictions  by 
comparison  with  measurements. 

Theories  of  combusting  monopropellant  sprays  are  not  highly  developed; 
therefore,  the  present  analysis  was  simplified.  Multiphase  effects  were  treated 
using  the  locally-homogeneous-fiow  (LHF)  approximation  of  multiphase  flow 
theory.  The  LHF  approximation  implies  that  relative  velocities  between  the  phases 
are  negligible  (negligible  slip),  yielding  a  single-fluid  formulation  which  is 
independent  of  the  specifics  of  spray  breakup  and  initial  drop-size  distributions. 
The  LHF  approach  generally  provides  useful  qualitative  information  concerning 
spray  structure  (Faeth,  1983,  1987;  Mao  ei  al.,  1980,  1981;  Shearer  et  al,  1979). 
Additionally,  \Mi  et  al  (1983,  1984)  report  encouraging  quantitative  performance 
of  LHF  analysis  for  pressure-atomized  injection  at  elevated  pressures— conditions 
which  are  most  often  encountered  for  combusting  monopropellant  sprays.  Finally, 
recent  work  in  this  laboratory  (Ruff  et  al,  1988),  suggests  that  LHF  analysis  is 
effective  for  near-injector  conditions  in  the  atomization  breakup  regime,  based  on 
predictions  and  measurements  of  the  properties  of  pressure-atomized  water  sprays 
in  still  air  at  normal  temperature  and  pressure. 


Address  correspondence  to;  G.  M.  Faelh,  217  Aerospace  Engineering  Building,  The  University  of 
Michigan.  Ann  Arbor,  Ml  4810V-2I40. Telephone: 1 3 1 3 >  764-7202. 

¥5 


41 


96 


T.  \V.  LEE  et  al. 


A  second  major  assumption  of  the  analysis  involved  use  of  the  laminar  flamelet 
concept  of  premixed  turbulent  single-phase  combustion,  proposed  by  Bray  ( 1978, 
1980).  Liquid  monopropellants  generally  have  relatively  high  burning  rates  at 
elevated  pressure  (McBratney,  1980.  1981),  which  implies  relatively  thin  reaction 
zones;  therefore,  the  limit  of  an  infinitely-thin  flame  was  used  to  provide  a  simple 
formulation  with  minimal  empiricism  (Bray,  1978, 1980). 

The  performance  of  the  analysis  was  evaluated  using  measurements  of  pressure 
atomized  sprays  from  a  round  injector  passage,  obtained  with  an  apparatus 
developed  by  Birk  and  Reeves  (1987).  The  measurements  consisted  of 
shadowgraphs  of  spray  flames  produced  by  a  hydroxl-ammonium-nitrate  (HAN)- 
based  liquid  monopropellant  at  ambient  pressures  of  6-8  MPa.  The  data  was 
reduced'  to  yield  mean  and  fluctuating  liquid  volume,  fractions  as  a  function  of 
position,  for  comparison  with  predictions. 

Experimental  and  theoretical  methods  are  described  in  the  next  two  sections  of 
the  paper.  Predictions  of  spray  structure  are  then  compared  with  measurements. 
The  paper  concludes  with  discussion  of  additional  theoretical  findings  for 
conditions  where  test  results  are  not  available. 


FIGURE  I  Injector  assembly. 
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EXPERIMENTAL  METHODS 

Experimental  methods  will  be  only  briefly  described,  more  details  are  provided  by 
Birk  and  Reeves  (1987).  The  experiments  were  carried  out  within  a  windowed 
cylindrical  test  chamber  having  an  inside  diameter  of  57  mm  and  an  interior  length 
of  406  mm,  yielding  a  volume  of  roughiy  1050  cc.  The  axis  of  the  chamber  was 
vertical,  with  the  pressure-atomized  spray  injected  vertically  upward  along  the  axis. 
The  propeilant  was  injected  into  a  hot  precombusted  gas  mixture,  in  order  to 
simulate  injection  into  an  adiabatic  combustion  environment. 

The  arrangement  of  the  injector  assembly  is  illustrated  in  Figure  1.  The 
propellant  was  pressurized  and  fed  through  the  injector  passage  using  a 
hydraulically  driven  piston  arrangement,  with  the  hydraulic  system  ultimately 
actuated  by  a  pneumatic  system  (Birk  and  Reeves.  1987).  The  injector  passages 
were  round  (1  and  2  mm  diameter),  and  were  2.4  mm  long,  yielding  passage 
length-to-diameter  ratios  of  1.2  and  2.4.  Due  to  steps  in  the  propellant  feed  system 
and  effects  of  the  vena-contracta  at  the  passage  inlet,  however,  it  is  likely  that  the 
effective  length-to-diameter  ratios,  as  a  measure  of  flow  development  at  the 
passage  exit,  are  not  very  well  defined  and  are  probably  somewhat  larger:  they  are 
taken  to  be  in  the  range  2-10  in  the  following.  The  outputs  of  pressure 
transducers  (Kistler  601B)  were  recorded  to  provide  the  test  chamber  pressure 
and  the  pressure  differential  across  the  injector  passage. 

The  experiments  were  limited  to  the  HAN-based  monopropeilant,  LGP  1846, 
see  Table  I  for  propellant  properties.  The  hot  gas  within  the  chamber,  used  to 
simulate  an  adiabatic  combustion  environment,  was  produced  by  spark  igniting  a 
combustible  gas  mixture,  initially  at  1  MPa.  having  the  following  composition  (by 
volume):  hydrogen.  20  percent;  oxygen,  10  percent;  and  argon.  70  percent.  Argon 
served  as  a  diluent  to  prevent  detonation  of  the  gas  mixture.  Due  to  the  relatively 
large  surface-to-volume  ratio  of  the  test  chamber,  heat  losses  from  the  burning  gas 
mixture  were  significant,  yielding  gas  temperatures  and  pressures  of  roughly 
1800  K  and  5.5  MPa  after  combustion.  The  temperature  level  was  roughly  10 
percent  below  estimated  adiabatic  flame  temperatures  for  the  propellant  (see  Table 
I),  providing  a  resonable  approximation  adiabatic  combustion  conditions.  The 


TABLE  I 

Properties  of  LGP  1846* 


Pressure  iMPai 

1 

10 

100 

Reactunt  Properties 

Densitv  <kgy  m'i 

1430 

1432 

1454 

Product  Properties 

Densitv  ikg/m'i 

1.36 

13.5 

134  7 

Temperature  !Kl 

Mass  Fractions  t%f' 

2027 

2039 

2045 

Water  Vapor 

53.8 

53.9 

54.0 

Carbon  Dioxide 

24.3 

24,6 

24.7 

Nitrogen 

21.0 

21.0 

21.1 

'Adiabatic  constant  pressure  combustion  with  the  reaciant  at  298.15  K.  Reactant  composition  (%  by 
massi:  hydroxyl-ammonium  nitrate  tHANl.  60.8;  triethanol-ammonium  nitrate  (TEANL  19.2:  and 
water.  20,0. 

"Maior  species  only.  Minor  species  include  CO.  H,.  NO.  OH  and  O.. 
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chamber  pressure  increased  as  the  propellant  burned:  present  measurements  were 
obtained  at  chamber  pressures  of  6-8  MPa.  The  density  of  the  burned  gas  mixture 
is  greater  than  the  density  of  the  adiabatic  combustion  products  of  the  propellant, 
due  to  the  lower  temperature  and  the  presence  of  argon;  therefore,  the  ambient 
density  of  the  experiments  at  6-8  MPa.  approximates  ambient  densities  for 
adiabatic  propellant  combustion  at  pressures  of  10-13  MPa. 

The  combusting  spray  was  observed  using  motion  picture  shadowgraphs.  A 
Photec  camera  was  used,  with  Kodak  VNF  7240  film.  Typical  film  speeds  were 
5000  frames  per  second.  Backlighting  was  provided  by  flash  sources  synchronized 
with  the  camera.  Maximum  flash  durations  were  less  than  2  ,us;  therefore,  the 
image  of  the  spray  was  effectively  stopped  on  the  films. 

A  typical  spray  shadowgraph  is  illustrated  in  Figure  2,  The  spray  or  liquid- 
containing  region  appears  as  the  dark  irregular  zone  near  the  center  of  the 
photograph.  The  boundaries  of  the  spray  are  reasonably  well-defined.  The 
shadowgraphs  were  unusually  clear  in  comparison  to  combusting  hydrocarbon 
sprays  at  comparable  pressures  (Mao  et  aL,  1981).  Factors  responsible  for  this 
behavior  include  the  absence  of  particulates,  like  soot,  in  the  monopropellant 
combustion  products;  and  the  relatively  uniform  temperature  of  the  combustion 
products,  which  minimizes  sharp  density  gradients  that  cause  irregular 
shadowgraph  patterns.  The  dark  zones  are  wispy,  suggesting  that  the  liquid 
behaves  somewhat  like  tracer  particles  in  the  outer  region  of  the  liquid-containing 
region  that  can  be  seen  on  the  photographs.  This  appearance  tends  to  support  use 
of  the  LHF  approximation  as  a  point  of  departure  for  analysis. 


FIGURE  2  Photograph  of  combusting  monopropellani  spray  ( 1  mm  diameter  injector.  6-8  MPa). 


The  motion  picture  films  were  analyzed  to  yield  time-averaged  mean  and 
fluctuating  liquid  volume  fractions.  Maintaining  the  thin-fiame  concept,  each 
picture  frame  was  analyzed,  assigning  dark  zones  to  unbumed  liquid  reactant  and 
light  zones  to  gaseous  combustion  products.  For  each  test,  15-25  frames  were 
available  for  analysis  during  the  steady  flow  portion  of  the  spray  combustion 
process.  Separating  dark  and  light  zones  on  the  film  was  somewhat  subjective;  and 
since  the  measurements  correspond  to  line-of-sight  projections,  they  are  biased 
downstream  and  radially  outward  from  correct  point  measurements  of  mean  and 
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fluctuating  liquid  volume  fractions.  These  effects  were  not  quantified,  but  are  not 
felt  to  be  large  in  comparison  to  other  experimental  uncertainties,  as  will  be  seen 
subsequently. 

Measurements  to  be  considered  were  obtained  from  three  tests,  denoted  tests  2, 
4  and  7  in  the  following.  A  1.0  mm  diameter  injector  was  used  during  test  2,  while 
a  2.0  mm  diameter  injector  was  used  during  tests  4  and  7.  In  the  portions  of  the 
films  that  were  analyzed,  liquid  velocities  were  in  the  range  45-55  m/s  for  both 
injectors.  Based  on  liquid  properties  summarized  by  Birk  and  Reeves  (1987),  this 
yields  injector  Reynolds  numbers  in  the  range  1-2  x  104  and  Ohnesorge  numbers 
of  0.023  (test  2)  and  0.016  (tests  4  and  7).  The  Reynolds  numbers  are  sufficiently 
high  to  insure  reasonably  turbulent  flow,  while  the  Reynolds  number/Ohnesorge 
number  combinations  suggest  that  operation  for  all  tests  was  well  within  the 
atomization  breakup  regime  (Miesse,  1955;  Ranz,  1958;  Ruff  etaL,  1988). 


THEORETICAL  METHODS 

The  analysis  involves  the  LHF  and  thin-flame  approximations,  with  turbulent 
mixing  treated  using  a  Favre-average  turbulence  model  along  the  lines  of  Bray 
(1978,  1980)  and  Bilger  (1976).  This  approach  provides  a  useful  limit  and  is 
consistent  in  the  sense  that  both  multiphase  and  chemical  reaction  phenomena  are 
assumed  to  be  controlled  by  turbulent  mixing,  vastly  reducing  the  empiricism  and 
amount  of  input  data,  e.g.,  initial  drop  size  and  velocity  distributions,  chemical 
kinetic  properties,,  etc.,  needed  to  define  the  problem  (Bray,  1978,  1980;  Faeth, 
1983,  1987).  As  noted  earlier,  the  LHF  approach  has  proven  to  be  effective  for 
providing  at  least  qualitative  information  on  the  structure  of  pressure-atomized 
sprays  at  high  pressures  (Faeth,  1983,  1987;  Mao  etaL,  1980,  1981;  Shearer  et  al, 
1979;  and  Wu  et  al,  1983,  1984);  while  more  recent  work  (Ruff  et  al,  1988) 
suggests  that  LHF  analysis,  similar  to  the  approach  used  here,  has  potential  for 
useful  quantitative  predictions  for  the  near-injector  region  of  turbulent  pressure- 
atomized  sprays  in  the  atomization  breakup  regime— conditions  which  are  met  for 
present  tests.  Furthermore,  the  approach  does  not  seem  unreasonable  in  view  of 
the  wispy,  gas-like,  appearance  of  the  spray-containing  region  seen  on  the  spray 
shadowgraphs,  see  Figure  2. 

High  pressure  conditions,  which  are  generally  of  interest  for  monopropellant 
sprays,  are  favorable  for  application  of  the  thin-flame  approximation.  Strand 
burning  rate  measurements  reported  by  McBratney  (1980,  1981)  for  HAN-based 
monopropellants  suggest  liquid  regression  rates  of  20  mm/s  at  pressures  on  the 
order  of  10  MPa.  Using  typical  transport  properties,  this  implies  a  characteristic 
flame  thickness  on  the  order  of  1  /tm  (representative  of  both  the  liquid-phase 
preheat  zone  and  the  distance  from  the  liquid  surface  to  the  reaction  zone  in  the 
gas  phase).  The  liquid-containing  regions  of  the  sprays  had  characteristic 
dimensions  on  the  order  of  10-100  mm.  see  Figure  2;  therefore,  the  probability  of 
a  typical  spatial  position  being  within  the  reaction  zone  is  quite  small,  justifying  the 
thin-flame  approximation.  As  noted  earlier,  the  rather  sharp  distinction  between 
liquid  and  gas,  and  the  absence  of  refractive  index  variations  in  the  gas-containing 
regions,  based  on  photographs  like  Figure  2,  are  also  supportive  of  this  hypothesis. 

It  seems  probable  that  the  main  structural  features  of  turbulent  pressure- 
atomized  sprays  in  the  atomization  breakup  regime  are  generally  similar  for  both 
noncombusting  and  combusting  monopropellant  sprays.  Based  on  what  is  known 
about  noncombusting  sprays,  the  flow  should  involve  an  all-liquid  core  which  can 
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extend  an  appreciable  distance  from  the  injector,  surrounded  by  a  growing  two- 
phase  shear  layer  which  originates  very  close  to  the  exit  of  the  injector  (Chehroudi 
et  aL  1985;  Hiroyasu  et  aL,  1982;  Ruff  et  aL,  1988;  Wu  et  aL,  1983.  1984).  The 
shear  layer  contains  irregularly-shaped  liquid  elements,  ligaments,  and  drops. 
Upon  merging  the  LHF  and  thin-flame  approximations,  a  thin  tlame  is  assumed  to 
cover  all  these  surfaces;  the  all-liquid  core,  the  irregularly-shaped  elements,  the 
ligaments,  and  the  drops.  Except  for  very  near  the  liquid  surface,  the  liquid  is  at 
the  same  state  as  in  the  injector:  while  beyond  the  outer  edge  of  the  thin  tlame,  the 
gas  has  uniform  properties  equivalent  to  adiabatic  flame  conditions.  Locally,  the 
relative  velocities  of  the  phases  (slip)  is  small. 

This  picture  is  analogous  in  many  ways  to  a  gaseous  premixed  flame  in  the 
multiple  reaction  sheet  regime  defined  by  Williams  (1985).  except  that  breakup  of 
the  liquid  yields  well-defined  islands  of  the  reactant,  rather  than  high  localized 
values  of  flame  stretch  in  a  single-phase  How.  Effects  of  flame  curvature  and  high 
flame  stretch,  analogous  to  gaseous  flames,  are  still  present  in  the  spray,  e.g, 
reaction  in  the  immediate  vicinity  of  small  drops  and  drops  having  very  high 
relative  velocities  may  be  extinguished.  However,  due  to  the  small  flame  thickness 
for  present  conditions,  drops  having  no  slip  will  be  very  small  ( ca .  1  ,nm  in 
diameter)  when  they  extinguish  and  will  not  have  much  effect  on  the  mixing 
properties  of  the  flow,  while  large  slip  velocities  are  precluded  under  the  LHF 
approximation. 

Other  major  assumptions  of  the  analysis  are  as  follows;  (1)  steady  (in  the  mean) 
axisymmetric  flow  with  no  swirl;  (2)  low  Mach  number  flow  with  negligible 
potential  and  kinetic  energy  changes,  and  negligible  viscous  dissipation,  in  the 
mean;  (3)  boundary-layer  approximations  apply;  (4)  negligible  effects  of  radiant 
energy  exchange;  (5)  equal  exchange  coefficients  of  all  species  and  heat;  and  (6) 
high  Reynolds  numbers,  so  that  laminar  transport  is  negligible  in  comparison  to 
turbulent  transport.  Most  of  these  assumptions  are  satisfied  by  the  conditions  of 
present  experiments.  In  particular,  the  boundary  layer  approximations  are  justified 
in  view  of  the  large  aspect  ratio  of  the  liquid-containing  region,  see  Figure  2.  The 
small  dimensions  of  the  liquid-containing  region  and  high  near-injector  velocities 
imply  low  radiation  numbers  and  thus  negligible  effects  of  radiant  energy 
exchange.  The  assumption  of  equal  exchange  coefficients  of  all  species  and  heat  is 
always  suspect  under  the  LHF  approximation,  since  rates  of  diffusion  of  even 
small  drops  are  small  (Faeth,  1983,  1987).  However,  laminar  exchange  is  not  likely 
to  be  very  important  at  Reynolds  numbers  on  the  order  of  104  and  this 
approximation  seems  no  worse  than  other  aspects  of  the  LHF  and  thin-flame 
limits. 

Under  these  assumptions,  Bray  (1978,  1980)  shows  that  flow  properties  can  be 
found  by  solving  governing  equations  for  conservation  of  mass,  momentum  and 
the  reaction  progress  variable,  in  conjunction  with  second  order  turbulence  model 
equations  for  turbulence  kinetic  energy  and  its  rate  of  dissipation.  The  governing 
equations  then  take  the  following  general  form  (Bray,  1980;  Bilger,  1976;  Jeng  and 
Faeth,  1984): 


1  d 


1  d  /  n,  d<f>\ 


—  (pdp)  +  -  -~[rpv<t>)  =  -  —  r—  —  +  S,, 


dx 


dr 


r  dr  \  o,  dr} 


U) 


where  \  (for  conservation  of  mass),  it,  c,  k  and  f.  The  source  terms,  5„ 
appearing  in  Eq.  (1)  are  summarized  in  Table  II,  along  with  the  empirical  constants 
used  during  present  computations.  The  turbulent  visosity  was  calculated  as  usual: 
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(2) 

The  boundary  conditions  for  Eq.  (1)  involve  symmetry  at  the  axis,  and  a  constant- 
property  ambient  environment,  e.g.: 

rm0,  ~  =  0\  r~°o.  ^=-0(for<*/c);  c-  1.  (3) 

dr 

Except  for  C^the  empirical  constants  appearing  in  Table  II  are  the  same  as  past 
analysis  of  single-  and  multi-phase  round  jets  in  this  laboratory  (Faeth.  1983. 
1987;  Jeng  and  Faeth,  1984).  These  values  were  selected  to  match  measurements 
for  a  variety  of  constant  and  variable  density  noncombusting  single-phase  round 
jets  (Faeth,  1983;  Jeng  and  Faeth,  1984).  Values  used,  however,  are  not  very 
different  from  those  used  during  early  work  with  k-f  turbulence  models 
(Lockwood  and  Naguib,  1975). 


TABLE  II 

Source  terms  in  the  governing  equations 


fi 

s. 

1 

0 

/i 

0 

c 

C,pc{  1  -  c)e/k 

k 

H,\dwdr\'-  -  pc 

e 

[C„pAdwdrf  -  C,.pt)c/k 

C. 

C„ 

Q 

C,  <7,  a, 

a. 

0.09 

t.-u 

1.87 

1.87  1.0  1.3 

0.7 

Values  of  C,  depend  on  the  shape  assumed  for  the  reacting-mode  probability 
density  function  and  several  turbulence  modeling  constants,  as  discussed  by  Bray 
(1980).  Considering  various  probability  density  functions  and  estimates  of 
empirical  constants  from  Bray  (1980),  Spalding  (1971),  Lockwood  and  Naguib 
(1975)  and  Jeng  and  Faeth  (1984),  yielded  values  of  C,  of  roughly  two.  Thus, 
C,=  Cf:=1.87  was  adopted  during  present  calculations,  as  a  baseline.  The 
influence  of  the  value  of  C,  on  predictions  was  evaluated  by  parameter  sensitivity 
calculations,  to  be  discussed  later. 

It  was  pointed  out  earlier  that  injector  exit  conditions  were  not  measured  and 
that  the  flow  involved  complicated  internal  injector  passages  which  precluded 
accurate  estimates  of  exit  conditions.  Furthermore,  computations  indicated  that 
spray  properties  were  strongly  inflenced  by  injector  exit  conditions.  In  order  to 
highlight  this  effect,  a  range  of  initial  conditions  was  studied.  This  included  limiting 
conditions  of  slug  flow,  with  a  uniform  exit  velocity  and  low  turbulence  intensities, 
and  fully-developed  pipe  flow.  In  addition,  calculations  were  carried  out  for  a 
range  of  length-to-diameter  ratios  of  the  injector  passage,  assuming  a  clean  entry 
with  no  vena  contracta.  Except  for  slug  flow,  profiles  of  mean  velocities  and 
turbulence  quantities  for  these  conditions  were  obtained  from  Hinze  (1959)  and 
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Schlichtlng  (1979),  given  the  flow  rate  and  the  diameter  of  the  passage.  The  initial 
shear  layer  for  the  slug  flow  calculations  was  assumed  to  have  a  thickness  of  0.1 
percent  of  the  passage  diameter,  while  properties  in  this  layer  were  specified 
similar  to  past  work  (Faeth.  1983). 

Under  the  thin-flame  approximation,  the  process  only  has  two  scalar  states, 
namely,  unbumed  reactant  liquid  and  completely  reacted  gaseous  products  of 
adiabatic  combustion.  The  properties  of  these  states  are  summarized  in  Table  t  for 
a  range  of  pressures.  Density  is  the  only  liquid  (reactant)  property  needed  for  the 
calculations:  the  values  given  in  Table  I  allow  for  liquid  compressibility  and  were 
provided  by  Friedman  (1986).  Combustion- product  properties  were  found  by 
assuming  thermodynamic  equilibrium  for  an  adiabatic  constant-pressure 
combustion  process  at  the  stated  pressure,  using  the  Gordon  and  McBride  (1971) 
computer  code. 

The  thin-flame  approximation  only  admits  a  double-delta  function  probability 
density  function  for  the  reaction  progress  variable.  Then,  mean  and  fluctuating 
scalar  properties  are  functions  of  c.  which  is  known  from  the  solution  of  Eq.  (1). 
The  functions  for  various  scalar  properties  are  as  follows  (Bray,  1978,  1980): 


0“  A»(l  -c)+  (4) 

f  =  ( $,p„(  1  -  c)  +  0„ac'Wp*(  1  -  c)  +  p{)c),  ( 5 ) 

f2~(k-tJc(  1-c),  (6) 

p\ fa-  )’c(  l  -  cMpyjoJ,  '  (7) 


where  p,  needed  to  solve  Eq.  (1),  can  be  found  by  setting  <j>  =  p  in  Eq.  (5).  Time- 
averaged  mean  and  fluctuating  liquid  volume  fractions  are  of  particular  interest, 
since  these  quantities  were  provided  by  the  experiments.  The  equations  for  these 
properties  are  as  follows: 

d;-  p[  1  -  c)/pn,  (8) 

^7-p;c(l-c)W.).  (9) 

The  calculations  were  performed  using  a  GENMIX  algorithm  due  to  Spalding 
(1978).  The  large  density  variations  in  the  flow  created  problems  of  stability  and 
numerical  accuracy,  requiring  a  much  finer  grid  than  is  usually  needed  for  single¬ 
phase  flows.  Results  reported  here  involved  360  crosstream  grid  nodes;  with 
streamwise  step  sizes  limited  to  0.3  percent  of  the  current  flow  width.  Doubling 
the  number  of  grid  nodes,  in  both  the  crosstream  and  streamwise  directions, 
changed  the  predictions  by  less  than  0.7  percent;  therefore,  the  numerical  accuracy 
far  exceeds  expected  uncertainties  of  other  aspects  of  the  analysis. 


RESULTS  AND  DISCUSSION 
Evaluation  of  Predictions 

Predictions  and  measurements  of  time-averaged  liquid  volume  fractions  along  the 
spray  axis  are  plotted  as  a  function  of  normalized  distance  from  the  injector  exit  in 


48 


COMBUSTING  MONOPROPELLANT  SPRAYS 


103 


Figure  3.  The  measurements  for  both  injector  diameters  are  very  similar  when 
plotted  in  this  manner.  By  itself,  this  finding  suggests  a  mixing-controlled  process 
for  the  conditions  of  the  experiments.  Similarly,  predictions  are  also  relatively 
insensitive  to  changes  of  injector  diameters  and  Reynolds  numbers  over  the  range 
of  the  tests:  therefore,  only  a  single  line  is  shown  for  each  injector  exit  condition.  As 
noted  earlier,  predictions  were  carried  out  for  a  pressure  of  L0  MPa,  in  order  to 
match  the  ambient  density  of  the  experiments  at  6-8  MPa.  However,  later  results 
will  show  that  the  effect  of  pressure  on  the  predicitions  for  this  range  of  conditions 
is  small;  therefore,  using  pressures  of  6-8  MPa  for  the  calculations  would  not 
materially  influence  the  results. 


2  4  6  8  10  20  4  0  6  0  80  100  200 

x/d 

FIGURE  3  Predicted  and  measured  variation  of  a,  along  axis. 

In  contrast  to  other  variables,  the  degree  of  flow  development  at  the  injector  exit 
has  a  strong  influence  on  the  predictions  illustrated  in  Figure  3.  For  example,  if  the 
characteristic  length  of  the  combustion  process  is  represented  by  the  condition 
where  a/r  =  0.l.  this  length  is  roughly  Five  times  longer  for  slug  flow  (Ud^O)  than 
for  fully-developed  flow  {L/d=  «>).  Due  to  uncertainties  concerning  experimental 
injector  exit  conditions,  this  large  effect  precludes  definitive  assessment  of 
predictions,  or  any  attempt  to  optimize  Cr  Nevertheless,  estimates  of  the  properties 
of  the  test  injectors  suggest  an  effective  L/d  in  the  range  2-10;  therefore,  it  is 
encouraging  that  predictions  at  these  limits  generally  bound  the  measurements.  The 
slopes  of  predictions  and  measurements  are  somewhat  different,  however,  in  the 
region  where  d*.  decreases  most  prominently.  This  may  be  due  to  line-of-sight 
biasing  and  difficulties  in  distinguishing  low  concentrations  of  gas  or  liquid  from 
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the  other  phase  on  the  photographs.  Effects  of  initial  flow  properties  may  also  be  a 
•  factor.  Further  study,  and  tests  with  better-defined  injector  flow  conditions, 
required  to  quantify  these  effects. 

The  strong  effect  of  initial  conditions  on  the  predicted  flow  properties  illustrated 
in  Figure  3  is  perhaps  not  surprising.  It  is  widely  recognized  that  changes  in  initial 
conditions  for  single-phase  turbulent  jets  can  significantly  influence  the 
development  of  the  flow  near  the  injector  (typically,  x/d<  20-30).  The  large 
density  of  the  liquid,  in  comparison  to  the  gas,  in  a  pressure-atomized  spray,  clearly 
provides  the  potential  for  carrying  effects  of  conditions  at  the  injector  exit  farther 
into  the  flow  field. 

Present  findings  concerning  the  strong  effect  of  initial  flow  development  on  the 
properties  of  combusting  monopropellant  sprays  are  also  supported  by  recent 
measurements  and  analysis  of  noncombusting  pressure-atomized  sprays  due  to 
Ruff  et  al.  (1987).  Their  measurements  show  that  fully-developed  flow  at  the 
injector  exit  causes  much  faster  mixing  than  slug  flow,  similar  to  the  predicted 
trends  illustrated  in  Figure  3.  Furthermore.  LHF  analysis,  similar  to  the  present 
approach  but  naturally  not  considering  premixed  combustion,  provided 
reasonably-good  estimates  of  the  effect  of  flow  development  in  the  injector  exit  for 
operation  in  the  atomization  breakup  regime— the  regime  of  interest  here. 

Direct  experimental  proof  of  effects  of  initial  flow  development  on  the  mixing 
properties  of  combusting  monopropellant  sprays  is  clearly  needed,  but  the 
preceding  discussion  indicates  that  there  is  some  evidence  that  the  predicted 
behavior  is  illustrated  in  Figure  3  is  real.  This  has  important  practical  implications, 
since  the  findings  suggest  that  fully-developed  injector  flows  can  substantially 
reduce  combustion  volumes.  This  can  be  used  to  provide  either  a  more-compact 
combustion  chamber,  or  to  reduce  tendencies  for  combustion  instabilities  by 
minimizing  the  amount  of  unreacted  propellant  within  the  combustion  chamber  at 
any  instant. 

Predicted  and  measured  time-averaged  fluctuating  liquid  volume  fractions  along 
the  axis  of  the  sprays  are  illustrated  in  Figure  4.  Both  predictions  and 
measurements  reach  a  maximum  value  of  o’,,.  =  0.5  at  afr  =  0.5,  which  is  a 
fundamental  requirement  of  an  intermittent  (on/off)  property  at  the  thin-flame 
limit.  Predictions  and  measurements  were  not  influenced  significantly  by  different 
injector  sizes  and  Reynolds  numbers,  similar  to  findings  discussed  in  connection 
with  Figure  3.  As  before,  predictions  for  L/d=  2  and  10  tend  to  bound  the 
measurements.  However,  the  measured  lengths  of  the  <?,<.  profiles  in  the  streamwise 
direction  are  narrower  than  predicted,  possibly  due  to  line-of-sight  and 
photographic-contrast  biasing,  as  discussed  earlier. 

Radial  profiles  of  predicted  and  measured  time-averaged  mean  and  fluctuating 
liquid  volume  fractions  are  illustrated  in  Figures  5  and  6.  The  results  are  plotted  as 
a  function  of  r/x,  the  similarity  variable  for  turbulent  jets,  at  various  distance  ( xld] 
from  the  injector.  The  experimental  scatter  of  the  tests  is  greater  for  the  radial 
profiles  than  the  results  plotted  in  Figures  3  and  4  for  properties  along  the  axis, 
particularly  near  the  end  of  the  liquid-containing  .region  [xJd-  20)  where  the 
streamwise  variation  of  flow  properties  is  quite  large,  see  Figure  3.  Nevertheless, 
effects  of  injector  diameter  and  Reynolds  number  are  seen  to  crudely  scale  as  a 
mixing-controlled  process, 

Predictions  are  illustrated  in  Figures  5  and  6  for  the  fully-developed  and  slug 
flow  limits.  The  predictions  scale  in  the  same  way  as  the  measurements;  however, 
the  measured  profiles  are  broader  than  either  of  the  limiting  predictions,  even 
though  LHF  predictions  generally  overestimate  the  width  of  the  liquid-containing 
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FIGURE  4  Predicted  and  measured  variation  of  a,  along  axis. 

region  (Faeth.  1983,  1987;  Mao  et  al,  1980,  1981;  Ruff  et  at..  1988).  There  are 
several  possible  explanations  for  this  behavior.  First  of  all,  the  predictions  shown 
are  not  expected  to  properly  bound  measurements  in  the  radial  direction,  based  on 
the  results  illustrated  in  Figure  3.  Slug  flow  mixes  much  more  slowly  than 
measured,  accounting  for  its  relatively  small  spread  in  the  radial  direction  in 
comparison  to  measurements.  Fully-developed  flow  mixes  much  more  rapidly  than 
measured,  and  volume  fractions  predicted  for  fully-developed  flow  are  much  lower 
than  measured  values  everywhere,  particularly  for  the  results  illustrated  in  Figure  5. 
Thus  use  of  intermediate  flow  development  conditions,  e.g.,  Ud=  2  or  10,  yields 
better  quantitative  agreement  with  measurements,  although  such  matching  was  not 
pursued  since  injector  exit  conditions  were  uncertain  in  any  event.  Next,  there  was 
evidence  of  unstable  flapping  of  the  jet  as  a  whole.  Such  flapping  does  not  have  a 
significant  effect  on  properties  in  the  streamwise  direction,  but  does  increase  the 
apparent  radial  spread  of  the  flow  through  a  mechanism  that  is  not  a  direct  effect  of 
turbulence.  Finally,  somewhat  broadened  measured  profiles  are  expected,  due  to 
line-of-sight  and  photographic-contrast  biasing,  as  noted  earlier.  In  view  of  these 
considerations,  the  order-of-magnitude  and  the  trends  of  the  predictions  are 
reasonable.  However,  a  wider  range  of  experimental  conditions  with  improved 
control  of  initial  conditions  are  needed  for  a  more  definitive  evaluation  of  the 
present  analysts 
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FIGURE  5  Predicted  and  measured  radial  variation  of  a,. 


Additional  Predictions 

Since  predictions  were  at  least  in  qualitative  agreement  with  the  measurements,  the 
analysis  was  exploited  to  examine  some  of  the  general  properties  of 
monopropellant  sprays  using  the  HAN-based  propellant.  The  influence  of  ambient 
pressure  and  the  state  of  flow  development  at  the  injector  exit  can  be  seen  from  the 
results  illustrated  in  Figures  7  and  8.  Predictions  of  o?fc  and  cc  for  combusting 
sprays  are  plotted  at  ambient  pressures  ranging  from  10-400  MPa.  which  covers 
the  range  of  pressures  encountered  for  most  applications  of  monopropellant  spray 
flames. 
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FIGURE  6  Predicted  and  measured  radial  variation  of  a',. 


The  findings  for  fully-developed  injector  exit  conditions  are  illustrated  in  Figure 
7.  The  Favre-averaged  reaction  progress  variable,  cp  is  a  measure  of  the  mass  of 
propellant  reacted  in  the  flow.  Pressure  clearly  has  a  strong  influence  on  cn  with  the 
reaction  nearing  completion  much  closer  to  the  injector  at  high  pressures.  This 
behaviour  is  caused  bv  higher  entrainment  rates  of  the  jet  at  higher  pressures, 
which  increases  the  rate  of  reaction  at  the  turbulent  mixing-controlled  limit 
considered  during  the  analysis.  Existing  measurements  suggest  that  the  relative 
entrainment  rates  of  turbulent  jets  are  proportional  to  (a/A=)i/2  >n  the  similarity 
region  of  the  flow  see  Ricou  and  Spalding  (1961)  for  single-phase  jets,  and  Faeth 
( 1983,  1987)  for  multiphase  jets  under  the  LHF  approximation.  Jet  development  in 
the  near-injector  region  tends  to  modify  this  behavior  for  the  results  illustrated  in 
Figure  7;  however,  the  trend  is  crudely  followed,  e.g.,  cc  tends  to  approach  unity  at 
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FIGURE  7  Mean  scalar  properties  along  axis  at  various  pressures  for  fully-developed  flow  at  injector 


FIGURE  8  Mean  scalar  properties  along  axis  at  various  pressures  for  slug  flow  at  injector  exit. 
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x/d*  15  and  60  for  ambient  pressures  of  400  and  40  MPa.  In  contrast,  the  effect  of 
pressure  on  dfr  distributions  is  relatively  small  for  pressures  in  the  range 
10-40  MPa.  This  is  due  to  compensation  of  increased  entrainment  rates  by 
reduced  phase-density  ratios  as  the  pressure  is  increased. 

The  most  obvious  difference  between  the  results  for  slug  flow  in  Figure  8  and  for 
fully-developed  flow  in  Figure  7  is  the  dominating  effect  of  an  extended  all-liquid 
potential  core  for  slug  flow.  Such  all-liquid  cores  have  been  observed  for  pressure- 
atomized  sprays  by  Hiroyasu  et  at.  ( 1982),  Chehroudi  et  at  (1985)  and  Ruff  et  al. 
(1988);  therefore,  this  behavior  is  not  unexpected.  The  presence  of  the  core  tends 
to  delay  the  development  of  the  flow,  particularly  at  low  pressures.  However,  in 
spite  of  the  core,  the  penetration  of  the  sprays,  represented  by  cc  approaching  unity, 
still  tends  to  scale  according  to  (pt/px)1  -,  similar  to  the  entrainment  properties  of 
fully-developed  jets.  In  contrast,  the  effect  of  pressure  on  d/r  differs  from  results  for 
fully  developed  flow,  since  the  presence  of  the  all-liquid  core  causes  afl.  to  vary 
significantly  with  pressure  over  the  whole  range  considered  in  Figure  8. 

Predicted  isochors  (for  df=0.01.  which  is  representative  of  minimum  liquid 
volume  fractions  generally  resolvable  from  photographs)  are  illustrated  in  Figure  9. 
Results  are  shown  for  fully-developed  and  slug  flow  initial  conditions  at  pressures 
in  the  range  10-400  MPa.  The  plots  are  distorted  to  improve  their  readability,  e.g., 
the  radial  scale  is  expanded  by  a  factor  of  ten  in  comparison  to  the  streamwise 
scale.  The  liquid-containing  region  is  relatively  narrow  in  comparison  to  its  length 
at  all  pressures,  suggesting  that  use  of  the  boundary-layer  approximations  is 
adequate  for  these  flows.  This  is  very  helpful,  since  the  turbulent  premixed  flame 
analysis  has  been  most  successful  for  boundary-layer  flows  (Bray,  1978,  1980). 
Increasing  pressure  acts  to  reduce  the  size  of  the  region  bounded  by  the  d|*0.01 
isochor,  for  both  boundary  conditions  and  all  pressures.  However,  the  effect  of 
pressure  on  the  variation  of  afn  at  higher  values  of  this  parameter,  is  relatively  small 
for  fully-developed  flow  at  pressures  in  the  range  of  10-40  MPa,  as  noted  earlier, 
see  Figure  7. 


FIGURE  9  Predicted  isochors  id,  =  0.01 1  at  various  pressures  (or  fully-developed  and  slug  flows  at  the 
injector  exit. 
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Parameter  Sensivitity 

In  addition  to  the  extent  of  flow  development  at  the  injector  exit,  uncertainties 
concerning  turbulence  properties  at  the  injector  exit,  and  for  other  parameters  of 
the  analysis,  contribute  to  uncertainties  of  predictions.  This  was  evaluated  by 
parametrically  varying  these  quantities  in  order  to  observe  their  effect. 

Table  III  is  a  summary  of  a  portion  of  the  parameter  sensitivity  results.  The 
percent  change  of  a  given  output  parameter  is  tabulated  for  either  a  50  percent 
reduction  or  a  100  percent  increase  of  various  input  parameters,  as  indicated  in  the 
table.  These  results  are  for  fully-developed  flow  at  10  MPa  and  x/d=2§.  This 
position  corresponds  to  conditions  near  the  end  of  the  liquid-containing  region,  see 
Figure  3.  The  most  sensitive  output  parameter  is  d,;,  with  mean  concentrations  of 
unreacted  liquid,  ( I  -  c)„  and  the  mean  velocity  ratio,  u/udh  being  relatively 
insensitive  to  the  parameter  variations.  The  effect  of  pressure  is  not  very  significant 
for  any  of  the  output  variables  at  this  pressure;  however,  this  would  not  be  the  case 
at  higher  pressure  levels,  see  Figures  7  and  8.  The  initial  turbulence  kinetic  energy 
has  the  strongest  influence  on  a,,,  with  a  50  percent  reduction  of  k,t  causing  almost 
a  300  percent  increase  of  d,r  This  highlights  the  importance  of  turbulent  mixing, 
which  increases  with  increasing  in  the  near-injector  region,  for  the  mixing- 
controlled  turbulent  reaction  limit  considered  in  the  calculations.  Predictably,  a 
factor-of-two  reduction  of  C,  causes  roughly  a  100  percent  increase  in  afn  implying 
that  this  constant  can  be  calibrated  using  measurements,  if  the  fundamental 
approximations  of  the  analysis,  like  the  effect  of  slip,  are  appropriate.  The  other 
parameters  considered,  q,  and  Re,  exhibit  relatively  small  effects  on  output 
properties  for  these  conditions. 


TABLE  HI 

Effect  of  parameter  variations-1 


Input 

parameter 

Outpui  parameter!0/,, increase) 

ii  -a, 

li,  ti,„ 

c 

110 

2 

_  1 

k, 

290 

4 

5 

e.i 

-37 

-3 

-4 

P~ 

_  1 

-5 

-0 

Re 

17 

1 

1 

For  fully-developed  injector  exit  conditions  al  10  MPa.  Output  parameter  evaluated  at  x/d  “  20 
"Parameter  variations  as  follows:  C,»  k „  -  fc,K  2.  f„  =  tv 2,  p„  =*  and  /?*,,"  2 Re„, 


CONCLUSIONS 

The  major  conclusions  of  the  investigation  are  as  follows: 

1)  Use  of  the  locally-homogenous  flow  and  thin  laminar  flamelet 
approximations  yielded  encouraging  agreement  with  measurements;  however, 
uncertainties  concerning  injector  exit  conditions  for  the  experiments  precluded 
definitive  assessment  of  the  analysis. 
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2)  Pressure-atomized  monopropellant  spray  flames  appear  to  be  unusually 
sensitive  to  the  degree  of  flow  development  and  the  turbulence  levels  at  the  injector 
exit;  fully-developed  flows  with  enhanced  turbulence  intensities  require 
significantly  smaller  combustion  volumes  than  slug  flows  with  low  initial  turbulence 
intensities. 

3)  Predictions  indicate  that  pressure-atomized  monopropellant  spray  flames 
generally  satisfy  the  boundary-layer  approximations,  except  at  very  high  pressures 
where  phase  densities  become  comparable.  This  is  helpful  since  premixed  turbulent 
flame  analysis  has  been  most  successful  for  flows  of  this  type,  and  the  degree  of 
empiricism  needed  for  predictions  is  minimized  (Bray,  1978, 1980). 

4)  Present  predictions  at  the  mixing-controlled  limit  suggest  that  combustion 
volumes  decrease  with  increasing  pressure  due  to  increased  entrainment  rates; 
however,  this  effect  is  modified  appreciably  by  the  degree  of  flow  development  at 
the  injector  exit. 
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NOMENCLATURE 

c  reaction  progress  variable 

C,  constants  in  turbulence  model.  Table  II 

d  injector  exit  diameter 

k  turbulence  kinetic  energy 

L  length  of  injector  passage 

p  pressure 

r  radial  distance 

Re  injector  exit  Reynolds  number 

St  source  source  term.  Table  II 

u  streamwise  velocity 

v  radial  velocity 

x  streamwise  distance 

Greek  Letters 
a 

E 

p 

o, 

<t> 

Subscripts 
c 

f 

3 
0 

Superscripts 

(-),  (-’)  time-averaged  mean  and  root-mean-squared  fluctuating  properties 
(”).("*)  Favre-averaged  mean  and  root-mean-squared  fluctuating  properties. 


centerline  value 
liquid-phase  property 
gas-phase  property 
injector-exit  condition 


phase  volume  fraction 

rate  of  dissipation  of  turbulence  kinetic  energy 

turbulent  viscosity 

density 

turbulent  Prandtl/Schmidt  number 
generic  property 
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Monopropellant  Sprays,"  by  T.-W.  Lee,  L.-K.  Tseng  and  G.M.  Faeth;  AIAA 
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Abstract 

The  drop  and  spray  combustion  properties  of  the  HAN- 
based  monopropeilant  LOP  1845  were  studied.  Drop  burning 
rates  were  measured  with  drops  supported  in  a  combustion  gas 
environment  at  pressures  of  0.2-7 D  MPa.  Some  internal  gasifi¬ 
cation  of  drops  —  causing  swelling,  partial  bursting,  and  mi¬ 
croexplosions  —  was  observed  throughout  this  region  but  these 
disturbances  decreased  with  increasing  pressure.  Effective  drop 
burning  rates  (including  effects  of  both  surface  gasification  and 
bursting)  were  relatively  constant,  ca.  10  mm/s,  and  were  con¬ 
sistent  with  earlier  strand  burning  rate  measurements  of  gelled 
propellant.  Pressure-atomized  combusting  sprays  were  studied 
in  combustion  gas  environments  at  pressures  of  3-9  MPa.  The 
liquid-containing  region  was  significantly  larger  than  earlier 
measurements  of  Birk  and  Reeves,  as  well  as  predictions  based 
on  the  locally-homogeneous-flow  approximation  of  multiphase 
flow  theory,  In  conjunction  with  drop  trajectory  calculations, 
based  on  present  measurements  of  drop  burning  rates,  these 
findings  suggest  significant  effects  of  separated  flow  in  com¬ 
busting  HAN -based  monopropeilant  sprays. 

Nomenclature 

Cq  =  drop  drag  coefficient 

d  =  injector  diameter 

dp  =  drop  diameter 

k  =  turbulence  kinetic  energy 

K  p  =  drop  bumi  ng  rate 

L  =  length  of  injector  passage 

Oh  =  Ohnesorge  number 

p  =  pressure 

r  =  radial  distance 

rv,  =  drop  radius 

Re  =  Reynolds  number 

t  =  time 

u  =  stream  wise  velocity 

We  =  Weber  number 

x  =  stream  wise  distance 

otf  =  liquid  volume  fraction 

e  =  rate  of  dissipation  of  turbulence  kinetic  energy 

p  =  viscosity 

p  =  density 

<J  =  surface  tension 

Subscripts 

c  =  centerline  value 

f  =  liquid- phase  property 

g  =  gas-phase  property 

p  =  drop  property 

o  =  injector  exit  condition 

«•  =  ambient  condition 

Superscripts 

(')  =  tune-averaged  property 

(")  =  Favre-averaged  property 
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Introduction 

Combusting  monopropeilant  sprays  have  applications  to 
regenerative  liquid-propellant  guns,  throttable  thrustors,  and  un¬ 
derwater  propulsion  systems.  The  objective  of  the  present  study 
was  to  experimentally  investigate  aspects  of  monopropeilant 
spray  combustion,  seeking  to  extend  earlier  theoretical  results 
obtained  in  this  laboratory-1-1  Two  spray  processes  were  con¬ 
sidered,  as  follows:  (1)  the  combustion  properties  of  individual 
drops  supported  in  combustion  gas  environments  at  pressures  of 
0.2-7.0  MPa;  and  (2)  the  structure  of  pressure-atomized  com¬ 
busting  sprays  in  combustion  gas  environments  at  pressures  of 
3-9  MPa.  The  new  measurements  were  used  to  assess  the  im¬ 
portance  of  separated-flow  phenomena  within  pressure-atomized 
combusting  monopropeilant  sprays,  i.e.,  effects  of  finite  relative 
velocities  and  transport  rates  between  the  phases.  Similar  to  our 
earlier  work,1'1  the  investigation  was  limited  to  a  hydroxyl-am¬ 
monium  nitrate  (HAN)-based  monopropeilant  (LGP  1845) 
which  is  of  interest  for  several  high-pressure  monopropeilant 
combustion  systems. 

Individual  drop  burning  rates  arc  needed  for  fundamental 
consideration  of  the  properties  of  combusting  monopropeilant 
sprays.  Earlier  studies  relevant  to  drop  burning  rates  of  HAN- 
based  monopropellants  have  included  measurements  of  strand 
burning  rates4-6  and  the  burning  rates  of  individual  drops  in 
heated  environments.7'11  McBratney4*5  measured  strand  burn¬ 
ing  rates  of  HAN-based  monopropellants  at  pressures  of  7-100 
MPa.  The  propellant  liquid  was  gelled  with  2  weight  percent 
Kelzan  in  order  to  stabilize  turbulent-like  disturbances  of  the  liq¬ 
uid  surface  that  are  normally  encountered  during  strand 
combustion  tests  at  high  pressures.  The  strand  burning  rates  of 
gelled  LGP  1845  were  high  (ca.  20  mm/s)  and  the  pressure  de¬ 
pendence  was  relatively  weak  (ca.  p0-1).  A  frothy  region  was 
observed  at  low  pressures,  where  the  thermal  disturbance  of  the 
combustion  wave  extends  an  appreciable  distance  into  the  un- 
bumed  propellant,  suggesting  significant  reaction  in  the  con¬ 
densed  phase  for  these  conditions.  While  these  results  are  valu¬ 
able,  however,  the  use  of  a  gelling  agent  raises  questions  con¬ 
cerning  its  influence  on  the  process.  Vosen6  measured  strand 
burning  rates  of  two  ungelled  HAN-based  monopropellants, 
LGP  1846  and  a  9.1  molar  solution  of  HAN  and  water,  at  pres¬ 
sures  of  7-30  MPa.  The  burning  rates  of  both  propellants  were 
very  high,  100-250  mm/s,  and  liquid  surfaces  were  clearly  dis¬ 
turbed,  indicative  of  turbulent- like  instability  of  burning  liquid 
strands  normally  seen  at  high  pressures;  therefore,  these  results 
are  difficult  to  interpret  to  find  the  fundamental  combustion 
properties  of  the  propellants. 

Zhu  and  Law7  studied  the  drop  combustion  properties  of 
LGP  1845  and  other  HAN-based  propellants,  in  combustion 
gases  at  1 170  K  and  1  atm.  The  drops  were  observed  to  heat  up 
with  no  radius  change  at  first,  then  gasify  from  the  surface  for  a 
time  (with  surface  regression  rates  of  ca.  0.2  mm/i),  and  finally 
burst  when  the  drop  diameter  had  decreased  by  roughly  15  per¬ 
cent  Beyer8'9  and  Beyer  and  Teague10  studied  the  combustion 
of  LGP  1 846  drops  supported  in  nitrogen  at  temperatures  of 
570-920  K  and  pressures  of  0. 1-8.2  MPa.  These  observations 
yielded  results  similar  to  Zhu  and  Law:7  after  a  heat-up  time  and 
a  period  of  relatively  slow  surface  gasification  (0.2  mm/s  at  730 
K  and  1  MPa)  the  drops  often  burst  --  particularly  the  larger 
drops.  Both  sets  of  drop  experiments  suggest  that  bulk  liquid 
reaction  and  microexplosioni  may  be  important  for  combustion 
of  HAN-based  monopropellants  but  drop  environment  tempera  - 
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tur's  were  low  in  comparison  to  the  adiabatic  combustion  'em- 
d  raturl  of  the  monopropall  am,  (ca.  2150  K);  therefore,  the 
drops  may  not  have  ignited  in  a  manner  representative  of  spray 
combustion. 

Earlier  theoretical  work  in  this  laboratory  addressed  liq¬ 
uid  surface  and  spray  properties  of  combusting  HAN-based 
monopropellants.  >'3  Analysis  of  liquid  surface  properties, !-2 
indicated  relatively  high  liquid  surface  temperatures  fin  the  range 
800  -  1050  K  for  pressures  greater  than  10  MPa)  and  unusually 
high  pressures  for  the  liquid  surface  to  reach  its  thermodynamic 
critical  point  (250  MPa  with  an  estimated  uncertainty  of  50  per¬ 
cent).  The  high  surface  temperatures  of  the  liquid  surface  pro¬ 
vides  greater  potential  for  significant  effects  of  chemical  reaction 
in  the  bulk  liquid  than  most  monopropdlants,  helping  to  explain 
observations  of  microcxplosions  reported  in  Refs.  7-10.  Fur¬ 
thermore,  the  high  critical  combustion  pressure  suggests  that 
spray  combustion  of  HAN-based  monopropellants  involves 
subcritical  combustion  with  a  drop-containing  combusting  spray 
for  most  applications. 

The  earlier  analysis  of  combusting  HAN-based  mono¬ 
propellant  sprays,1'3  was  based  on  the  locally-homogeneous- 
flow  (LHF)  approximation  of  multiphase  flow  theory,  i.e.,  the 
assumption  that  velocity  differences  between  the  phases  are 
negligible  at  each  point  in  the  flow;11-13  and  the  thin  laminar 
flamelet  approximation  of  turbulent  premixed  flame  theory,  pro¬ 
posed  by  Bray.14'15  Turbulent  mixing  was  estimated  using  a 
Favre-averaged  k-e  turbulence  model,  with  empirical  constants 
established  from  measurements  in  noncombusting  variable-den¬ 
sity  round  jets,16'17  however,  the  constants  used  are  very  simi¬ 
lar  to  early  proposals  based  on  constant-density  turbulent 
flows.18  The  performance  of  the  analysis  was  evaluated  using 
the  measurements  of  Birk  and  Reeves19  for  pressure  atomized 
combusting  LOP  1846  sprays  at  pressures  of  6-8  MPa.  There 
was  encouraging  agreement  between  predictions  and  measure¬ 
ments,  however,  predictions  were  very  sensitive  to  the  degree  of 
flow  development  at  the  injector  exit  which  was  not  known  very 
well;  therefore,  this  assessment  was  not  definitive.  Later  mea¬ 
surements  of  noncombusting  pressure-atomized  sprays  by  Ruff 
et  al.13  established  the  strong  sensitivity  of  spray  properties  to 
the  degree  of  flow  development  at  the  jet  exit  and  observed  rea¬ 
sonably  good  performance  of  LHF  analysis  in  the  dense-spray 
region  (liquid  volume  fractions  greater  than  0.2)  near  the  injector 
exit  for  atomization  breakup.  However,  these  measurements 
also  disclosed  significant  deficiencies  of  LHF  analysis  for  other 
breakup  regimes  and  in  the  dilute  portion  of  the  spray  -  the  last 
being  in  general  agreement  with  other  recent  evaluations  of  the 
LHF  approach  for  dilute  sprays.11’12 


Tilt  present  investigation  sought  to  extend  past  work 
concerning  both  drop  and  spray  combustion  of  HAN-b=«scd 
monopropellants.  Drop  combustion  was  observed  using  an  ap¬ 
proach  similar  to  Beyer10  for  pressures  of  0.2-7  MPa,  however, 
the  drop  environment  more  closely  matched  the  gas  temperature 
of  a  combusting  monopropellant  spray.  Measurements  of  spray 
properties  were  undertaken  seeking  to  confirm  the  measurements 
of  Birk  and  Reeves,19  while  considering  a  broader  range  of  ex¬ 
perimental  conditions  at  pressures  of  3-9  MPa.  The  new  spray 
measurements,  in  conjunction  with  both  LHF  analysis  and  drop 
trajectory  calculations  based  on  the  present  drop  burning  rate 
measurements,  were  used  to  assess  the  importance  of  separated- 
flow  phenomena  for  these  flows. 

Drop  Combustion 

Experimental  Methods 

Apparatus.  Figure  1  is  a  sketch  of  the  drop  combustion 
test  apparatus.  The  supported-drop  technique  was  used  with  the 
drops  exposed  to  gases  in  the  post-flame  region  of  a  premixed 
burner  which  was  operated  within  a  pressure  vessel.  The  pres¬ 
sure  vessel  had  an  inside  diameter  and  length  of  130  and  430 
mm  and  was  fitted  with  two  25  mm  diameter  quartz  windows  so 
that  the  drops  could  be  observed. 

Tlie  premixed  burner  had  a  diameter  of  10  mm  with  a 
stainless  steel  screen  (0.17  mm  diameter  wires,  2000  wires/m, 
square  pattern)  to  help  stabilize  the  flame.  The  gas  flow  rates  of 
the  premixed  burner  were  metered  and  controlled  with  critical 
flow  orifices  and  pressure  regulators.  Burner  operating  times 
were  short,  just  sufficient  to  stabilize  the  premixed  flame  and 
complete  the  drop  combustion  test.  Burner  gas  flows  were 
initiated  and  terminated  with  solenoid  valves  while  the  burner 
was  ignited  with  an  exploding  wire.  The  pressure  rise  of  the 
chamber  (measured  with  a  pressure  transducer)  was  small  in  the 
period  when  the  burner  was  operating,  ca.  5  percent;  therefore, 
the  chamber  pressure  was  set  by  backfilling  it  with  air.  The 
properties  of  the  post-flame  region  of  the  premixed  burner 
roughly  approximated  the  temperatures  of  adiabatic  combustion 
of  the  monopropellam,  but  contained  significantly  lower  con¬ 
centrations  of  water  vapor,  see  Table  1  for  the  combustion 
product  properties  of  LGP  1845  and  the  burner  gases  (denoted 
burner  1  and  2). 

The  drop  support  assembly  is  illustrated  in  Fig.  2.  The 
drops  were  mounted  on  quartz  fibers,  50-150  pm  in  diameter, 
with  the  bottom  end  of  the  fiber  flame  polished  to  a  bead  of 
somewhat  larger  diameter  to  help  support  the  drop.  The  drop 
was  surrounded  with  a  retractable  shield  to  protect  it  from  tran¬ 
sients  when  the  premixed  flame  was  ignited.  Once  the  premixed 


Table  1.  Combustion  Product  Properties* 


Mixture 

LGP  1 845b 

Burner  ld 

Simulent  Gases 

Burner  2d 

Spray 

Temperature  (K) 

2150 

2295 

2230 

2790 

Composition  {%  by  volume)0 
H20 

69.2 

18.8 

18.0 

19.8 

CO2 

12.9 

9.1 

8.9 

-- 

n2 

17.4 

71,2 

72.6 

38.3 

Ar 

"" 

40.9 

*  Computed  for  10  MPa  using  the  Gordon  and  McBride20  algorithm,  but  effects  of  dissociation  are  small. 
b  Reactant  composition  (%  by  mass);  HAN,  63.2, 1FAN,  20;  and  IHO,  1 6.8. 
c  Major  species  only.  Minor  species  include  CO,  Hi,  NO,  OH  and  Or 
d  Volume  flow  rate  of  burner  gases  (cold)  of  6.28  x  10" 5  m3/s. 
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n  stabilized  the  shield  was  rapidly  retracted  by  fusing  its 

'  .  c  ■  a  jo  that  the  unbalanced  pressure  force  on  the  shield 
j^fed  it  ,0  onc  side  of  the  pressure  vessel  where  it  was  stopped 

by^  rubber  cushion. 


[ttslfnmentation,  Drop  diameter  was  measured  as  a 
function  of  time  using  backlighted  high-speed  motion  picture 
photographs.  The  arrangement  of  the  illuminating  and  camera 
system  is  illustrated  in  Fig,  1.  The  drops  were  backlighted  by  a 
continuous  arc  source,  using  a  condensing  lens  to  direct  the  light 
to  a  diffusion  screen  located  at  one  of  the  windows,  The  pho¬ 
tographs  were  obtained  using  a  high-speed  motion  picture  cam¬ 
era  operating  ai  roughly  1000  pictures  per  second  which  incor¬ 
porated  an  internal  timing  marker. 


Results  and  Discussion 


Drop  combustion  at  low  pressures  yielded  very  irregular 
variations  of  drop  diameter  as  a  function  of  time  due  to  bubble 
formation  and  bursting  within  the  drops.  Some  typical  results  at 
low  pressures  are  illustrated  in  Fig,  3.  Drop  diameters  are  plot¬ 
ted  as  a  function  of  time  for  five  tests  at  0.51  MPa  with  initial 

drop  diameters  in  the  range  580-770  )im  and  a  300  pm  diameter 
bead  on  the  quartz  fiber  to  help  support  the  drop,  The  origin  of 
these  plots  is  somewhat  arbitrary  since  the  motion  of  the  re¬ 
tractable  shield  disturbed  the  premixed  flame  causing  it  to  flap 
for  a  time;  therefore,  the  time  when  the  drop  was  finally  sub¬ 
merged  in  the  post  flame  gases  was  uncontrolled  and  variable. 
The  results  in  Fig.  3  show  swelling  of  the  drop  due  to  bubbles 
in  the  liquid  in  every  case.  The  bubbles  would  periodically 
burst,  canying  off  some  of  the  liquid,  and  occasionally  the 
bursting  of  a  bubble  was  sufficiently  severe  to  carry  off  all  of  the 
liquid.  At  these  low  pressure  conditions,  internal  reaction  and 
bursting,  with  some  mechanical  removal  of  liquid  caused  by  the 
bursts,  appears  to  be  the  main  mechanism  for  the  reduction  of 
the  drop  diameter. 


The  degree  of  drop  swelling  due  to  the  presence  of  bub¬ 
bles  in  the  bulk  liquid,  and  the  severity  of  drop  bursting,  de¬ 
creased  as  the  pressure  was  increased.  Some  typical  results  at 
higher  pressures  are  illustrated  in  Fig.  4.  Drop  diameter  is  plot¬ 
ted  as  a  function  of  time  for  five  representative  tests  at  2.1  MPa 

with  initial  drop  diameters  in  the  range  520-680  nm  and  a  200 

(am  diameter  bead  on  the  quartz  fiber.  As  before,  the  time  origin 
is  arbitrary  due  to  effects  of  initial  disturbances  on  the  premixed 
flame.  All  these  conditions  exhibited  some  degree  of  internal 
bubble  formation,  however,  effects  of  bubbles  bursting  were 
relatively  mild  and  complete  bunting  of  drops  was  not  observed 
at  pressures  of  2. 1  MPa  and  higher. 

Reduced  effects  of  internal  bubbles  at  high  pressures  ap¬ 
pears  to  be  largely  caused  by  increased  gas  density  so  that  a 
given  degree  of  bulk  liquid  reaction  yields  a  lower  volume  of 
gas;  this  reduces  bubble  sizes  and  growth  rates  which  tends  to 
reduce  the  severity  of  bunting  phenomena,  Counter  to  this  is 
the  fact  that  liquid  surface  temperatures  tend  to  increase  with  in¬ 
creasing  pressure  for  the  present  range  of  conditions  (reaching  a 
maximum  at  roughly  25  MPa):*>2  this  is  expected  to  increase 
rates  of  bulk  liquid  reaction. 

The  time  period  of  drop  swelling,  or  relatively  constant 
drop  diameten,  was  irregular  due  to  uncertainties  concerning  the 
time  when  the  drop  was  submerged  in  the  combustion  gas  envi¬ 
ronment.  However,  the  period  when  the  drop  diameter  de¬ 
creased  was  analyzed  to  obtain  effective  drop  burning  rates. 
Plots  of  drop  diameter  as  a  function  of  time  in  the  period  where 
the  drop  diameter  is  decreasing  are  illustrated  for  the  present  test 
conditions  in  Fig.  5.  The  origins  of  these  plots  are  arbitrary 
since  the  data  has  been  plotted  to  overlap  in  the  region  where  the 
drop  diameter  is  decreasing.  In  addition,  conditions  where  the 
drops  burst  completely  at  low  pressures  have  been  excluded, 
Results  at  low  pressures  show  wide  variations  due  to  significant 
effects  of  bubble  swelling  and  bursting  but  the  diameter  traces 
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become  more  regular  and  repeatable  at  high  pressures.  These 
data  were  fitted  to  determine  effective  burning  rates  for  the 
drops,  Kp  =  -drp/dt;  the  fits  arc  also  illustrated  in  Fig.  5. 

Present  effective  burning  rates  are  plotted  as  a  function 
of  pressure  in  Fig.  6.  These  results  are  for  drop  diameters  in  the 
range  300-1200  and  include  effects  of  both  internal  reaction 
forming  bubbles  which  burst,  mechanically  removing  some  liq¬ 
uid,  as  well  as  conventional  gasification  at  the  surface  of  the 
drop,  This  combination  of  effects  causes  ihe  effective  burning 
rate  to  be  highest  at  the  lowest  pressure,  where  bursting  domi 
nates  the  process,  and  then  to  show  relatively  little  change  with 
pressure  over  most  of  the  region  considered  during  present  tests. 
The  strand  burning  results  of  McBratney4'5  and  Vosen6  are  also 
illustrated  in  Fig.  6.  The  present  results  are  a  crude  extension  of 
McBratney's4'5  measurements  of  gelled  propellants  at  higher 
pressures.  The  results  of  Vosen6  are  much  higher  than  the  rest 
of  the  measurements  due  to  effects  of  liquid  surface  disturbances 
of  burning  liquid  strands  at  high  pressures,  noted  earlier. 

Snrav  Combustion 
Experimental  Methods 

Apparatus.  The  present  spray  combustion  test  apparatus 
was  similar  to  the  arrangement  used  by  Birk  and  Reeves.19  A 
sketch  of  the  apparatus  appears  in  Fig.  7.  The  experiments  were 
conducted  in  the  same  chamber  as  the  drop  combustion  tests. 
The  combustion  environment  was  produced  by  filling  the  cham¬ 
ber  with  a  combustible  mixture  and  then  igniting  it  with  two 
sparks  to  achieve  the  combustion  gas  properties  summarized  in 
Table  1  (denoted  spray).  The  pressure  of  the  spray  tests  was 
adjusted  by  varying  the  initial  pressure  of  the  combustible  gas 
mixture  since  combustion  of  this  gas  approximated  a  constant 
volume  process.  The  combustible  gas  mixture  had  temperatures 
that  were  somewhat  greater  than  adiabatic  combustion  tempera¬ 
ture  of  the  monopropellant. 

The  spray  was  pressure  atomized  using  injectors  having 
exit  diameters  of  0.31,  0,58,  1.08  and  1,17  mm.  The  inlet  of 
the  injectors  had  baffles,  to  control  any  swirl  in  the  liquid,  and 
smooth  entries,  to  reduce  effects  of  cavitation.  Injectors  having 
length-to-diameter  ratios  of  2,  17  and  42  were  considered  since 
earlier  work  indicated  that  the  degree  of  flow  development  at  the 
injector  exit  influenced  spray  mixing  properties. 2.3.13  The  in¬ 
jectors  were  directed  vertically  upward. 

A  test  was  run  by  placing  a  propellant  sample  (3-4  ml)  in 
the  fuel  delivery  tube  and  filling  the  injector  passage  up  to  its 
exit.  A  cap  was  then  placed  over  the  exit  to  prevent  gas  inflow 
when  the  chamber  was  filled  with  the  combustible  gas  mixture 
and  further  pressurized  as  this  gas  burned.  The  propellant  flow 
was  initiated  by  venting  nitrogen  from  an  accumulator  into  the 
fuel  delivery  tube  by  opening  a  solenoid  valve.  Once  the  pres¬ 
sure  of  the  propellant  was  greater  than  the  chamber  pressure,  the 
cap  popped  off  and  the  resulting  propellant  flow  generated  a 
spray  in  the  hot  gas  mixture.  The  process  ended  when  all  the 
propellant  was  consumed,  The  injector  passage  continued  to  be 
purged  by  the  nitrogen  flow  from  the  accumulator  for  a  time  be¬ 
fore  the  accumulator  flow  was  ended. 

Instrumentation.  The  combusting  sprays  were  observed 
using  motion  picture  shadowgraphs  as  illustrated  in  Fig.  7. 
Backlighting  was  provided  by  a  flash  lamp  source,  ca.  1  ps 
flash  duration,  which  was  synchronized  with  the  camera; 
therefore,  the  image  of  the  spray  was  effectively  stooped  on  the 
film.  The  shadowgraphs  were  recorded  with  a  lo  mm  high¬ 
speed  camera  operating  at  roughly  1000  pictures  per  second, 
using  Tri-X  negative  film.  The  camera  optics  yielded  a  25  mm 
diameter  field  of  view;  therefore,  it  was  necessary  to  adjust  the 
position  of  the  injector  to  observe  the  full  length  of  the  liquid 
containing  region. 
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Since  propellant  combustion  does  not  produce  particu¬ 
lars  and  gas  temperatures  are  relatively  uniform  in  monopro- 
Ivliant  soray  flames,  the  boundaries  of  the  spray  were  reason- 
b lv  well  defined  --  similar  to  past  measurements  of  Birk  and 
Reeves  19  The  fi'ms  were  to  yield  mean  and  fluctuat¬ 

es  time  averaged  liquid  volume  fractions,  assigning  dark  zones 
to  unhurried  liquid  reactant  and  light  zones  to  gaseous  combus¬ 
tion  products.  For  each  test,  10-30  frames  were  available  for 
analysis  during  the  steady  flow  portion  of  the  spray  combustion 
process.  Separating  dark  and  light  zones  was  somewhat 
subjective;  and  since  the  measurements  correspond  to  line-of- 
sight  projections,  they  are  biased  downstream  and  radially  out¬ 
ward  from  correct  point  measurements  of  mean  and  fluctuating 
liquid  volume  fractions.  Predictions  were  analyzed  to  estimate 
the  line-of-sight  biases,  as  discussed  in  the  next  section, 

Test  Conditions.  Test  conditions  for  the  spray  combus¬ 
tion  tests  are  summarized  in  Table  2.  Most  of  the  injector  flows 
correspond  to  fully  developed  flow  at  the  injector  exit,  which 
due  to  the  Reynolds  number  of  the  passage  flow  corresponds  to 
turbulent  pipe  flow.  Injection  velocities  were  in  the  range  49-65 
m/s:  these  conditions  correspond  to  the  atomization  break-up 
regime,  i.e.,  a  drop-containing  shear  layer  begins  to  develop  at 
the  liquid  surface  immediately  at  the  injector  exit.12-11 

Theoretical  Methods 

Present  measurements  were  compared  with  the  mono¬ 
propellant  spray  combustion  analysis  developed  earlier  in  this 
laboratory.1'3  Drop  trajectory  calculations  were  also  carried  out 
in  order  to  help  assess  effects  of  separated-flow  phenomena. 
Both  analyses  are  described  in  the  following. 

Spray  Analysis.  The  main  features  of  the  spray  analysis 
will  be  only  briefly  described  in  the  following,  original  sources 
should  be  consulted  for  details. 

The  analysis  involves  use  of  the  LHF  approximation  of 
multiphase  flow  theory1 1-13  and  the  thin  laminar  flamelet  ap¬ 
proximation  of  premixed  turbulent  flame  theory. t4-15  Turbulent 
mixing  was  treated  using  a  Favre-averaged  turbulence 
model. 14-12  This  approach  provides  a  useful  limit  since  both 
multiphase  and  chemical  reaction  phenomena  are  controlled  by 
turbulent  mixing  which  minimized  the  empiricism  needed  for 
predictions,  e.g.,  initial  drop  size  and  velocity  distributions, 
chemical  kinetic  properties,  etc.,  are  not  needed  to  define  the 
problem.  The  main  limitation  of  the  LHF  approximation  is  that 
its  use  generally  tends  to  overestimate  the  rate  of  development  of 
sprays,  particularly  in  dilute-spray  regions  far  from  the  injector 


exit. 11-13  However,  Ruff  et  al.13  find  lhai  the  LHF  approach, 
using  the  present  turbulent-mixing  model,  provided  reasonably 
good  estimates  of  mixing  properties  in  the  near-injector  dense- 
spray  region  of  nonevaporating  pressure-atomized  sprays  in  the 
atomization  breakup  regime  -•  conditions  that  are  representative 
of  present  tests. 

The  thin  laminar  flamelet  approximation  implies  that  het¬ 
erogeneous  monopropellant  flames  cover  all  liquid  surfaces. 
Except  for  very  near  the  liquid  surface,  the  liquid  is  at  the  same 
state  as  in  the  injector  while  beyond  the  outer  edge  of  the  thin 
flame  the  gas  has  uniform  properties  equivalent  to  adiabatic 
flame  conditions  noted  in  Table  1  Under  the  LHF  approxima¬ 
tion,  relative  velocities  between  the  phases  (slip)  are  neglected. 

Other  major  assumptions  of  the  analysis  are  as  follows; 
(1)  steady  (in  the  mean)  axisymmetnc  flow  with  no  swirl,  t2) 
low  Mach  numbers  with  negligible  potential  and  kinetic  energy 
changes,  and  negligible  viscous  dissipation;  (3)  boundary- layer 
approximations  apply;  (4)  negligible  effects  of  radiant  energy 
exchange;  (5)  equal  exchange  coefficients  of  all  species  and  heat; 
and  (6)  high  Reynolds  numbers,  so  that  laminar  transport  is 
negligible  in  comparison  to  turbulent  transport.  Justification  of 
these  assumptions  is  presented  in  Refs.  1  and  3. 

Under  these  assumptions,  flow  properties  can  be  found 
by  solving  governing  equations  for  conservation  of  mass,  mo¬ 
mentum  and  reaction  progress  variable,  in  conjunction  with  sec¬ 
ond-order  turbulence  model  equations  for  turbulence  kinetic  en¬ 
ergy  and  its  rate  of  dissipation.14’1 3  The  formulation,  all  em¬ 
pirical  constants  used  in  the  turbulence  model,  and  the  method  of 
solution,  can  be  found  in  Ref.  3. 

The  predictions  were  also  used  to  estimate  potential  ef¬ 
fects  of  line-of-sight  biasing  on  the  measured  distributions  of 
liquid  volume  fractions  using  a  stochastic  approach  developed 
for  radiation  calculations  in  this  laboratory.21  Knowing  the 
time-averaged  probability  density  function  of  the  reaction 
progress  variable  along  paths  through  the  flow,  the  reaction 
progress  variable  was  simulated  for  a  series  of  statistic  ally-inde¬ 
pendent  eddies  along  the  path.  Counting  the  presence  of  any 
liquid  in  the  path  as  a  condition  which  would  block  the  light, 
giving  a  daifc  image  on  the  film,  yielded  estimates  of  time  aver¬ 
aged  mean  and  fluctuating  liquid  volume  fractions  for  the  path. 
"Hiis  procedure  has  not  been  calibrated  using  known  flows, 
however,  it  does  provide  at  least  a  qualitative  indication  of  po¬ 
tential  effects  of  line-of-sight  bias. 


Table  2.  Summary  of  Combusting  Spray  Test  Conditions 


Diameter 

(mm) 

L/d 

Flow 

Type4 

Amb.  Pres. 
(MPa) 

Pressure 

Drop  (MPa) 

Inj.  Velocity® 
(m/s) 

Re® 

Ohd 

Wet® 

Radial  Measurements; 
1.17  17 

FDF 

3.11 

2,02 

52.7 

12600 

0.021 

70700 

1.17 

17 

FDF 

7.07 

1.79 

49.3 

11800 

0.021 

62400 

1.17 

17 

FDF 

9.10 

1.81 

49.7 

11900 

0.021 

63300 

1.08 

2 

SF 

6.83 

1.93 

51.4 

11400 

0022 

62200 

0.58 

42 

FDF 

3.19 

2.07 

53.2 

6300 

0.030 

35800 

0.58 

42 

FDF 

6.15 

2.76 

61.6 

7300 

0.030 

47900 

Axial  Measurements; 
0.58  42 

FDF 

3.22 

2.33 

56.6 

6700 

0.030 

40400 

0.58 

42 

FDF 

6.16 

1.81 

49,5 

5900 

0.030 

30900 

0.58 

42 

FDF 

8.99 

2.14 

54.3 

6400 

0.030 

37100 

0.31 

42 

FDF 

6.51 

3.05 

64.8 

4000 

0.041 

27800 

bUnity  flow  coefficient. 

<JOh  =*  jif /  ■'/pfdo ,  o  *  0.0669  kg/s2. 


•FDF  =  fully  developed  flow;  SF  =  slug  flow. 
cRe  ■  pfU<3d/Uf,  Hf  =  0.007 1  kg/ms. 

«Wef  =  ppio^o. 
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ppp  Trajectory  Analysis.  Direct  assessment  of  the  ap¬ 
proximations  of  the  LHF  approach  for  the  monopropellant 
sprays  was  undertaken  using  drop  trajectory  calculations,  similar 
to  the  approach  used  by  Shearer  et  al.22  and  Mao  et  al.23  for 
nonpremixed  spray  flames.  These  calculations  were  limited  to 
drops  moving  along  the  axis  of  the  spray.  The  drops  were  as¬ 
sumed  to  be  always  in  contact  with  the  gas  phase  which  was 
taken  to  have  the  properties  summarized  in  Table  1.  Estimates 
of  the  gas  velocity  variation  along  the  axis  were  obtained  from 
the  LHF  analysis. 

Drop  trajectory  calculations  were  limited  to  deterministic 
calculations,  ignoring  effects  of  turbulence/drop  interactions; 
therefore,  mean  gas  velocities  from  the  LHF  analysis  were  used 
in  the  governing  equations  for  drop  motion.  Drops  were  as¬ 
sumed  to  be  surrounded  by  gas  immediately  at  the  injector  exit, 
ignoring  the  all-liquid  core  present  in  these  sprays.13  Effecis  of 
drop  heat-up  were  also  ignored:  the  drop  radius  was  assumed 
to  decrease  throughout  the  entire  trajectory  at  10  mm/s  -  based 
on  the  results  of  Fig.  6  for  the  present  test  range.  This  high 
burning  rate  implies  that  the  decomposition  flame  is  located  near 
the  drop  surface,  well  within  the  boundaries  of  the  flow  field 
around  the  drop;  therefore,  gas-phase  properties  used  to  estimate 
drop  drag  were  taken  to  be  ambient  gas  properties  and  effects  of 
forced  convection  or  drop  burning  rates  were  ignored.  Other 
aspects  of  the  analysis  were  similar  to  Refs.  22  and  23:  the  flow 
field  around  the  drop  was  assumed  to  be  quasi-steady;  virtual 
mass,  pressure-gradient,  Basset  history  and  gravitational  forces 
were  ignored;  swelling  of  the  drops  was  ignored;  and  drop  drag 
was  estimated  using  the  standard  drag  correlation  for  solid 
spheres. 

Under  these  assumptions,  the  governing  equations  of 
drop  motion  along  the  axis  are  as  follows: 


dxp/dt  =  Up 

(1) 

ddp/dt  =  -2  Kp 

(2) 

dup/dt  =  -3pg  Cp  1  Up  -  u  1  (up  -  u)/4 

(3) 

where 

CD=24  (l+Rcp^yRcp,  Rep  51000;  CD=0.44,  Rep>1000  (4) 

The  initial  condition  is  up  =  u0,  d„  =  dpo  and  xp  =  0  at  t  =  0. 
Equations  (1)  -  (3)  were  integrated  using  a  Runge-Kutta  algo¬ 
rithm. 

Results  ani  Discussion 

Due  to  strong  background  lighting,  it  was  not  possible  to 
determine  that  ignition  had  taken  place  from  flame  luminosity. 
Nevertheless,  ignition  was  readily  identified  from  the  pressure 
trace:  inert  liquid  (like  water  or  unignited  monopropellant  spray 
liquid)  caused  the  hot  combustion  products  of  the  premixed  gas 
flame  to  be  quenched  which  resulted  in  a  rapid  reduction  of  the 
chamber  pressure;  in  contrast,  energy  release  from  the  com¬ 
busting  monopropellant  spray  caused  an  increase  of  chamber 
pressure  in  the  period  when  the  propellant  was  flowing.  It  was 
possible  to  consistently  ignite  the  spray  at  pressures  as  low  as 
2.7  MPa,  however,  the  bulk  of  spray  measurements  were  ob¬ 
tained  at  pressures  of  3-9  MPa. 

Measured  and  predicted  time- averaged  liquid  volume 
fractions  along  the  axis,  5fc,  are  plotted  as  a  function  of  nor¬ 
malized  distance  from  the  injector  exit,  x/d,  in  Fig,  8.  Both  pre¬ 
sent  measurements  and  those  of  Birk  and  Reeves19  are  shown 
on  the  plot  Predictions  include  direct  values  of  5fc  u  well  as 
results  allowing  for  line-of-sight  biu,  u  noted  earlier. 


Test  conditions  used  by  Birk  and  Reeves19  were  similar 
to  present  test  conditions,  except  that  injector  L/d  were  in  the 
range  1. 2-2.4  and  the  injector  inlet  was  not  rounded  (see  Lee  et 
al.3  for  a  sketch  of  the  injectors),  and  injector  pressure  drops 
were  1. 5-2.0  times  higher  than  the  present  study.  The  motion 
picture  shadowgraphs  of  both  investigations  were  obtained  in  a 
similar  manner  and  were  analyzed  in  this  laboratory.  Each  set  of 
experimental  results  also  exhibits  a  significant  degree  of  internal 
consistency  and  repeatability  when  piolted  in  the  manner  of  Fig. 
8.  Finally,  pressure  traces  indicated  that  measurements  were 
obtained  for  combusting  sprays  for  both  studies.  Nevertheless, 
present  measurements  exhibit  a  much  longer  liquid-containing 
region  than  those  of  Birk  and  Reeves,19  e.g.,  5fc  =  0.5  at  x/d 
roughly  150  and  25  for  the  two  sets  of  measurements.  Specific 
reasons  for  these  differences  are  not  obvious  since  so  many 
features  of  the  two  studies  were  the  same,  however,  changes  in 
injection  properties  offer  the  most  plausible  explanation.  In  par¬ 
ticular,  the  sharper  injector  inlet  used  by  Birk  and  Reeves19 
could  have  caused  cavitation  in  the  injector  passage  resulting  in  a 
more  finely  atomized  spray  with  a  rapid  rate  of  radial  spread. 
Similarly  ,  the  injector  used  by  Birk  and  Reeves19  did  not  have  a 
flow  straightener  and  swirl  induced  in  the  injector  (low  passage 
could  have  resulted  in  unusually  high  radial  spread  rates;  al¬ 
though  the  fuel-injection  system  only  involved  rectilinear  motion 
and  doesn't  appear  to  be  fundamentally  prone  to  induce  swirl. 
Finally,  Birk  and  Reeves19  employed  somewhat  higher  injector 
pressure  drops  which  would  be  expected  to  yield  smaller  drop 
sizes  in  the  spray;  nevertheless,  spray  conditions  for  both  inves¬ 
tigations  were  in  the  atomization  breakup  regime  and  the  pres¬ 
sure  drop  increase  doesn't  appear  to  be  sufficient  to  explain  the 
differences  seen  in  Fig.  8  based  on  the  relatively  small  effect  of 
pressure  drop  variations  observed  during  this  investigation.  In 
any  event,  extensive  rechecking  of  measurements  using  the  pre¬ 
sent  injectors  could  not  reproduce  the  results  of  Birk  and 
Reeves.19 

Present  measurements  in  Fig.  8  are  roughly  similar  (in 
terms  of  x/d)  at  all  test  conditions,  with  the  downstream  limit  of 
the  liquid-containing  region  at  x/d  ca.  350.  Since  these  results 
involve  a  range  of  pressures  and  injector  diameters,  this  behav¬ 
ior  suggests  a  mixing-controlled  process  supporting  the  use  of 
LHF  analysis  -  a  conclusion  reached  in  Ref.  3,  based  on  the 
measurements  of  Birk  and  Reeves,19  Closer  examination  of  the 
data,  however,  reveals  trends  that  suggest  significant  separated- 
flow  effects.  First  of  all,  results  for  the  0.31  mm  diameter  in¬ 
jector  consistently  exhibit  higher  values  of  ctfe  at  a  particular  x/d 
than  the  0.58  mm  diameter  injector.  This  is  a  separated- flow 
property  since  drop  diameters  are  not  strongly  affected  by  injec¬ 
tor  diameters  while  drops  of  a  particular  size  must  penetrate  a 
certain  distance  in  order  to  disappear:  this  causes  in  a  tendency 
for  penetration  distances,  x,  to  be  constant  for  separated  flows 
rather  than  x/d,11  Another  effect  is  that  cifc  at  a  particular  x/d  is 
lower  for  a  chamber  pressure  of  8.99  MPa  than  the  other  pres¬ 
sures  considered  for  the  0.58  mm  diameter  injector  this  behav¬ 
ior  parallels  the  effective  burning  rate  results  of  Fig.  6  where 
drop  burning  rates  at  9  MPa  are  higher  than  for  pressures  in  the 
range  3-6  MPa,  which  are  roughly  the  same.  A  final  effect  is 
that  use  of  long  and  short  L/d  injectors  yielded  roughly  the  same 
results  while  mixing-controlled  flows  would  result  in  much 
faster  mixing  rates  for  the  long  L/d  injector,13 

Predictions  illustrated  in  Fig.  8  are  for  fully -developed 
flow  at  the  injector  exit,  which  corresponds  to  the  bulk  of  pre¬ 
sent  test  conditions.  Effects  of  ambient  pressure,  injector  diam¬ 
eter  and  injector  Reynolds  number  had  little  effect  on  the  predic¬ 
tions;  therefore,  only  single  lines  are  shown  for  results  with  and 
without  the  line-of-sight  bias  correction.  Comparing  predictions 
with  and  without  the  line-of-sight  bias  correction  indicates  sig¬ 
nificant  effects  of  bias  for  intermediate  values  of  5ft;  however, 
predictions  of  the  downstream  end  of  the  liquid-containing  re¬ 
gion  are  not  strongly  influenced  by  the  bias, 
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In  view  of  the  bias  uncertainties,  the  predictions  illus- 
oated  in  Fig.  8  are  in  fair  agreement  with  the  measurements  of 
Birk  and  Reeves.19  This  observation  prompted  earlier  encour¬ 
agement  concerning  the  value  of  the  LHF  and  thin  laminar 
flamelet  approximations  for  analyzing  flows  of  this  type,  How¬ 
ever,  comparison  of  predictions  with  present  measurements  im¬ 
plies  that  use  of  the  LHF  approximation  causes  the  rate  of 
development  of  the  spray  to  be  substantially  overestimated,  in 
agreement  with  most  other  evaluations  of  the  LHF  approxima¬ 
tion  for  sprays.11'13 

Radial  profiles  of  time-averaged  liquid  volume  fractions 
at  various  distances  from  the  injector  are  illustrated  in  Fig.  9.  All 
measurements  shown  in  the  figure  were  obtained  during  the 
present  investigation  Predictions  shown  on  the  figure  account 
for  line-of-sight  bias  and  are  for  fully-developed  flow  at  the  in¬ 
jector  exit.  Similar  to  results  along  the  axis,  predictions  were 
relatively  independent  of  test  conditions  and  only  a  single  line  is 
shown  for  each  streamwise  position.  Results  ignoring  line-of- 
sight  bias  are  narrower  than  the  plots  illustrated  in  Fig,  9,  how¬ 
ever,  the  outer  extent  of  the  liquid-containing  region  is  about  the 
same. 

Similar  to  results  along  the  axis,  the  measured  radial 
profiles  are  crudely  similar  for  all  the  test  conditions  when  plot¬ 
ted  in  the  manner  of  Fig.  9.  In  terms  of  r/x,  the  radial  similarity 
variable  of  turbulent  jets,  the  liquid-containing  region  extends  tc 
0.05-0.07,  rather  than  0.15  which  is  the  typical  width  based  on 
scalar  properties  in  turbulent  jets,  Predictions  provide  a  fair  es¬ 
timate  of  flow  widths  near  the  injector  exit  but  progressively  fail 
with  increasing  distance  from  the  injector  exit  —  tending  to  over¬ 
estimate  the  rate  of  development  of  the  flow,  This  behavior  is 
similar  to  other  evaluations  of  the  use  of  the  LHF  approximation 
for  both  nongasifying  and  gasifying  sprays. 

Potential  effects  of  separated  flow  are  examined  directly 
by  the  drop  trajectory  computations  illustrated  in  Fig.  10.  Drop 
velocities  and  diameters  along  the  axis  are  plotted  as  a  function 
of  distance  from  the  injector  for  an  ambient  pressure  of  10  MPa. 
Predictions  of  velocities  along  the  axis  from  the  LHF  analysis 
are  also  illustrated  on  the  plot,  as  a  reference.  Results  are  shown 
for  initial  drop  diameters  of  10,  20,  100  and  200  pm;  drops 
much  larger  than  200  pm  would  be  subject  to  secondary  breakup 
due  to  excessively-high  drop  Weber  numbers.12  Unlike  LHF 
predictions,  drop  trajectory  calculations  depend  on  the  initial  in¬ 
jector  diameter,  as  noted  earlier:  the  results  illustrated  in  Fig,  10 
are  for  an  injector  diameter  of  1.00  mm. 

The  results  illustrated  in  Fig.  10  clearly  show  significant 
effects  of  separated  flow.  The  LHF  predictions  exhibit  a  decay 
of  velocity  beyond  the  potential-core-like  region  which  is 
roughly  inversely  proportional  to  pressure  -  similar  to  single¬ 
phase  jets.  Due  to  the  small  diameter  of  the  injector,  this  results 
in  a  rapid  deceleration  rate.  Only  the  smallest  drops  (initial  di¬ 
ameters  of  10  pm  or  less)  have  sufficiently  fast  response  to  ap¬ 
proach  the  velocities  of  the  continuous  phase  throughout  most  of 
their  trajectory.  With  increasing  drop  size,  the  drops  progres¬ 
sively  overshoot  the  velocity  of  the  continuous  phase  and  only 
approach  it  again  toward  the  end  of  their  life,  when  they  become 
very  small.  Similarly,  the  drops  pass  beyond  the  end  of  the  liq¬ 
uid-containing  region  estimated  by  the  LHF  analysis  (taken  to  be 
Ofc  >  1CH  since  liquid  volume  fraction  never  formally  reaches 
zero  in  the  LHF  analysis  due  to  its  statistical  treatment),  Use  of 
the  drop  burning  rate  estimates  of  Fig.  6,  however,  yields  drop 
trajectories  extending  to  xAl  ca.  300  for  initial  drop  diameters  of 

200  pm.  This  is  comparable  to  present  measurements  of  the 
extent  of  the  liquid-containing  region  suggesting  that  the  results 
of  the  drop  trajectory  calculations  are  at  least  reasonable. 

Taken  together,  the  results  of  Figs.  8-10  suggest  signifi¬ 
cant  effects  of  separated  flow  for  combusting  HAN-based 
monopropellant  sprays  over  the  present  range  of  test  conditions. 


In  view  of  the  relatively  modest  variation  of  burning  rate  with 
pressure  seen  in  Fig,  6,  the  insensitivity  of  drop  drag  properties 
to  pressure,11-12  and  the  relatively  high  critical  combustion 
pressure  of  HAN-based  monopropellants  (ca.  250  MPa1'2)  it  is 
likely  that  separated-fiow  phenomena  are  important  for  com¬ 
busting  HAN-based  monopropellant  sprays  for  most  of  their 
range  of  application, 

Conclusions 

The  present  study  involved  measurements  of  the  com¬ 
bustion  properties  of  the  HAN-based  monopropellant  LGP 
1 845,  both  as  drops  and  sprays  in  combustion  gas  environments 
at  pressures  of  0.2-9  MPa.  The  spray  measurements,  and  drop 
trajectory  calculations  based  on  the  present  drop  burning  rate 
measurements,  were  used  to  evaluate  earlier  analysis  of  com¬ 
busting  monopropellant  sprays  based  on  the  locally  homoge¬ 
neous  flow  and  thin  laminar  flamelet  approximations,  due  to  Lee 
et  al.1*3  Major  conclusions  of  the  study  are  as  follows: 

1 .  Measurements  yielded  effective  drop  burning  rates  of  ca. 
10  mm/s  for  drop  diameters  of  300-1200  pm  and  pres¬ 
sures  of  0.2-7  MPa.  The  effective  drop  burning  rate  in¬ 
volved  both  reaction  within  the  bulk  liquid  causing  bub¬ 
ble  formation  and  bursting,  dominating  the  process  at 
low  pressures;  and  conventional  gasification  from  the 
drop  surface,  dominating  the  process  at  high  pressures; 
taken  together,  these  effects  cause  burning  rates  to  be 
relatively  independent  of  pressure  over  the  present  test 
range. 

2.  Present  measurements  of  drop  burning  rates  at  pressures 
of  0.7-7  MPa  are  generally  consistent  with  earlier  strand 
burning  rate  measurements  of  gelled  propellants  due  to 
McBratney4'5  at  pressures  greater  than  10  MPa. 

3.  Present  measurements  exhibited  a  much  larger  liquid 
containing  region  for  combusting  sprays  at  pressures  of 
3-9  MPa  than  the  earlier  measurements  of  Birk  and 
Reeves19  even  though  test  conditions  and  methods  of 

data  analysis  were  similar,  e.g.,  otfC  =  0.5  at  x/d  roughly 
150  and  25  for  the  two  sets  of  experiments.  Reasons  for 
these  differences  have  not  been  firmly  established  but 
different  injector  passage  conditions,  possibly  leading  to 
effects  of  cavitation, swirl  and  finer  atomization  for  the 
measurements  of  Ref.  19,  have  been  advanced  as  a  pos¬ 
sible  explanation. 

4.  While  earlier  evaluation  of  analysis  using  the  LHF  and 
thin  laminar  flamelet  approximations  appeared  promising 
based  on  the  measurements  of  Birk  and  Reeves;19  cur¬ 
rent  findings  suggest  that  this  approach  substantially 
overestimates  the  rate  of  development  of  the  flow  which 
is  consistent  with  recent  findings  for  other  pressure-at¬ 
omized  spray  processes.11'13  Separated  flow  phenom¬ 
ena  appear  to  be  important  combusting  for  HAN-based 
monopropellant  sprays  over  much  of  their  range  of  ap¬ 
plication. 
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Structure  and  Mixing  Properties  of  Combusting  Monopropellant  Sprays 


T.-W.  Lee*  and  G.  M.  Faeih* 
Depanmem  of  Aerospace  Engineenng 
The  University  of  Michigan,  Ann  Atbor,  Michigan 


Abstract 

The  spray  combustion  properties  of  the  HAN-based 
monopropeUant  LGP  1845  were  studied  both  theoretically  and 
experimentally.  Drop  size,  liquid  mass  flow  rate,  and  liquid 
mass  flux  distributions  were  measured  for  pressure-atomized 
sprays  in  the  atomization  breakup  regime,  burning  within  a 
combustion  gas  environment  at  pressures  of  4.5-5. 0  MPa. 
Two  separated  flow  models  were  evaluated  using  the  new 
measurements:  a  deterministic  separated  flow  model  where 
drop/turbulence  interactions  were  ignored,  and  a  stochastic 
separated  flow  model  where  drop/turbulence  interactions  were 
considered  using  random-walk  computations  for  drop  motion. 
When  based  on  burning  rates  found  from  earlier  single  drop 
experiments,  both  models  were  in  reasonably  good  agreement 
with  the  measurements.  Separated  flow  effects  are  quite 
important  for  these  sprays,  with  the  length  of  the  liquid- 
containing  region  being  relatively  independent  of  injector 
diameter  and  extending  roughly  300  mm  from  the  injector  exit 
for  injector  diameters  of  0.3  arid  0.6  mm. 

Nomenclature 

Cb  =  coefficient  in  drop  breakup  correlation 

Cd  =  drop  drag  coefficient 

Ci  =  constants  in  turbulence  model 

d  =  injector  diameter 

dp  =  drop  diameter 

dpavg  =  average  drop  diameter 

k  =  turbulence  kinetic  energy 

Kp  =  drop  burning  rate 

L«  =  characteristic  eddy  size 

tin  =  liquid  mass  flow  rate 

tit"  =  liquid  mass  flux 

mi  =  drop  mass  of  group  i 

rq  =  number  of  drops  per  unit  time  in  drop  group  i 

Oh  =  Ohnesorge  number,  PfUgci/jif 

r  =  radial  distance 

Re  =»  Reynolds  number,  pf  u^d/pr 

SMD  a  Sauter  mean  diameter 

S$  =  source  term 

t  =  time 

t«  =  characteristic  eddy  lifetime 

u  =  stream  wise  velocity 

ui  =  velocity  in  coordinate  direction  i 

v  =  radial  velocity 

vj  =>  volume  of  computational  cell  j 

2 

We^  =  liquid  Weber  number,  p  fu^d/o 

x  =  stream  wise  distance 

Xj  =  distance  in  coordinate  direction  i 

at  =  liquid  volume  fraction 
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£  =  rate  of  dissipation  of  turbulence  kinetic  energy 

p  =  absolute  viscosity 

Pt  =  turbulent  viscosity 

p  =  density 

0  =  surface  tension 

Oi  =  turbulent  Prandti/Schmidt  number 

$  =  generic  property 

Subscripts 

c  =  centerline  quantity 

f  =  liquid  property 

p  =  drop  property 

0  =  injector  exit  condition 

°°  =  ambient  condition 

Superscripts 

(_)  =  dme-averaged  quantity 

(  )'  =  time-averaged  root-mean-squared  fluctuating 

quantity 

Introduction 

Combusting  monopropellant  sprays  represent  an 
important  fundamental  multiphase  combusting  flow,  since 
they  are  the  premixed  counterpart  of  more  common  spray 
diffusion  flames.  Additionally,  combusting  monopropeUant 
sprays  have  applications  to  regenerative  liquid-propellant 
guns,  throttable  thrustors,  and  underwater  propulsion 
systems.  Motivated  by  these  considerations,  the  present 
investigation  considered  the  structure  of  combusting  pressure- 
atomized  monopropeUant  sprays  both  theoretically  and 
experimentally,  seeking  to  extend  earlier  work  in  this 
laboratory.  '-3  Similar  to  the  earlier  studies,1'3  the  inves¬ 
tigation  was  limited  to  a  hydroxyl-ammonium-nitrate  (HAN> 
based  monopropeUant  (LGP  1845)  which  is  of  interest  for 
several  high-pressure  monopropeUant  combustion  systems. 

The  importance  of  multiphase  flow  phenomena  is  a 
significant  issue  for  combusting  monopropeUant  sprays.  In 
particular,  most  applications  involve  high-pressure  combus¬ 
tion  (greater  than  10  MPa)  so  that  drops  can  approach  their 
thermodynamic  critical  point  where  multiphase  flow  charac¬ 
teristics  would  disappear.  Thus,  combustor  pressures 
required  to  reach  near-critical  conditions  have  been  theo¬ 
retically  studied  for  HAN-based  monopropellants. 1  It  was 
found  that  these  pressures  are  unusually  high,  ca.  250  MPa 
with  an  uncertainty  of  50  percent  due  to  uncertainties  of 
thermochemical  properties.  As  a  result,  multiphase  effects  are 
important  for  most  applications  involving  these  mono- 
propellants. 

The  importance  of  finite  interphase  transport  rates 
(separated-flow  effects)  is  also  a  significant  issue  for  high- 
pressure  combusting  sprays.  This  has  been  examined  during 
earlier  work  in  this  laboratory.2-3  Initial  work  involved 
development  of  a  locally-homogeneous  flow  (LHF)  model  for 
combusting  HAN-based  monopropeUams,  where  interphase 
transport  rates  are  assumed  to  be  infinitely  fast.2  This 
approach  yielded  reasonably  good  agreement  with  flow 
visualization  measurements  of  Birk  and  Reeves4  for  com¬ 
busting  LGP  1845. 2  However,  later  evaluation  in  this 
laboratory  suggested  significant  separated-flow  effects  using 
similar  measurements  as  well  as  drop  trajectory  calculations.3 
Thus,  the  need  for  a  separated-flow  treatment  of  these  sprays 
has  been  established  although  models  of  this  type  have  not 
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been  reported. 

Information  concerning  the  combustion  properties  of 
individual  drops  is  needed  to  treat  effects  of  separated  flow; 
therefore,  several  studies  relevant  to  drop  combustion  have 
been  undertaken  for  the  HAN-based  monopropellants.5-5-11 
This  includes  measurements  of  strand  burning  rates,5-7  as  well 
as  observations  of  individual  drop  combustion  in  heated  or 
combustion  gas  environments.5'8-11  Drop  combustion  at  low 
pressures  (less  than  2  MPa)  involves  considerable  subsurface 
liquid  reaction  with  bubble  formation  and  mechanical  removal 
of  liquid  by  microexplosions  being  the  main  mechanism  for 
the  reduction  of  the  drop  diameter.5-8-11  This  behavior 
follows  since  the  liquid  surface  temperatures  of  HAN-based 
monopropellants  are  unusually  high,  ca.  1000  K,  tending  to 
promote  liquid  phase  reactions.1  At  higher  pressures  (greater 
than  2  MPa)  within  combustion  gas  environments,  however, 
conventional  surface  gasification  as  a  heterogeneous  premixed 
flame  is  more  dominant.5-5,6  The  combined  effect  of  the  two 
mechanisms  yields  relatively  high  monopropellant  drop 
burning  rates,  10-40  mm/s,  with  a  weak  pressure  dependence, 
ca.  pl/5,  over  the  pressure  range  of  0.4-100  MPa.5 

Finally,  recent  studies  of  the  near-injector  (dense 
spray)  region  of  pressure  atomized  sprays  are  helpful  for 
identifying  initial  conditions  for  separated-flow  predictions  of 
spray  properties.12-16  It  is  generally  agreed  that  this  region 
involves  a  liquid  core,  somewhat  like  the  potential  core  of 
single-phase  jets,  surrounded  by  a  multiphase  mixing  layer. 
Flow  properties  approximate  estimates  based  on  the  LHF 
approximation  at  high  liquid  volume  fractions,  since  the  liquid 
dominates  mixing  properties  due  to  its  larger  density, 
however,  the  bulk  of  the  multiphase  mixing  layer  exhibits 
significant  effects  of  separated  flow.12-16  Finally,  when 
effects  of  liquid  turbulence  are  small,  average  drop  sizes  near 
the  liquid  surface  approximate  estimates  based  on  aerodynamic 
breakup  theory  due  to  Reitz  and  Bracco;12  while  larger  drops 
observed  when  liquid  turbulence  is  present  eventually  break 
up  to  yield  similar  drop  sizes.16 

The  present  investigation  sought  to  extend  past  work 
on  the  structure  of  combusting  HAN-based  monopropellant 
sprays.  Spray  combustion  was  observed  in  combustion  gas 
environments  near  5  MPa,  using  slide  impaction  to  measure 
drop  size,  liquid  mass  flow  rate  and  liquid  mass  flux 
distributions.  These  measurements  were  used  to  evaluate 
separated-flow  models  of  the  process  considering  two  limits, 
analogous  to  methods  used  for  evaporating  sprays:17 
deterministic  separated  flow  (DSF)  where  drop/turbulence 
interactions  arc  ignored,  and  stochastic  separated  flow  (SSF) 
where  drop/turbulence  interactions  are  considered  using 
random-walk  computations  for  drop  motion.  The  present 
description  of  the  study  is  brief,  more  details  and  tabulations 
of  data  can  be  found  in  Lee. 18 

Experimental  Methods 

Apparatus 

The  present  spray  combustion  apparatus  was  similar  to 
Lee  et  al 5  and  Birk  and  Reeves.4  A  sketch  of  the  apparatus 
appears  in  Fig.  1.  The  experiments  were  conducted  within  a 
windowed  pressure  vessel  having  an  inside  diameter  and 
length  of  130  and  430  mm.  The  sprays  were  injected 
vertically  upward  in  a  combustion  gas  environment.  The 
combustion  gas  environment  was  produced  by  filling  the 
chamber  with  a  combustible  gas  mixture  and  igniting  it  with  a 
spark  to  achieve  the  combustion  gas  properties  summarized  in 
Table  1  (these  properties  were  computed  assuming  thermo¬ 
dynamic  equilibrium).19  The  pressure  of  the  combusting 
spray  was  set  by  adjusting  the  initial  pressure  of  the 
combustible  gas  mixture,  since  combustion  of  this  gas 


approximated  a  constant  volume  process.5  The  adiabatic 
combustion  temperature  of  the  gas  mixture  was  somewhat 
greater  than  the  adiabatic  constant-pressure  combustion 
temperature  of  the  roonopropellant  (see  Table  l),  however, 
temperature  variations  in  this  range  do  not  have  a  large  effect 
on  drop  burning  rates.5 

The  spray  was  pressure-atomized  using  injectors 
having  the  diameters  listed  in  Table  1.  The  inlets  of  the 
injectors  were  baffled,  to  eliminate  swirl,  and  smooth,  to 
reduce  effects  of  cavitation.  Length-to-diameter  ratios  were 
42,  yielding  nearly  fully-developed  turbulent  ptp«  flow  at  the 
exit. 

The  apparatus  operated  by  placing  a  3-4  ml  sample  of 
the  propellant  in  the  injector  feed  bne  so  that  the  passage  was 
filled  right  up  to  the  exit.  A  cap  was  placed  over  the  exit  to 
prevent  gas  inflow  when  the  chamber  was  pressured  by  filling 
or  burning  the  combustion  gas  mixture,  Injection  was  initiated 
by  venting  nitrogen  from  an  accumulator  into  the  delivery 
tube.  Once  the  pressure  of  the  propellant  exceeded  the 
chamber  pressure,  the  cap  popped  off  and  the  spray  flowed 
into  the  combustion  gas  environment  The  process  ended 
when  all  the  propellant  was  injected,  after  which  the  injector 
passage  continual  to  be  purged  by  nitrogen. 

Instrumentation 

The  operation  of  the  apparatus  was  monitored  by 
measuring  injector  inlet  and  chamber  pressures  using  pressure 
transducers.18  Earlier  measurements  involved  flash  shadow¬ 
graph  motion  pictures  of  the  spray,  however,  present 
measurements  were  limited  to  slide  impaction. 

A  sketch  of  the  slide  impaction  system  appears  in  Fig. 
2.  The  arrangement  involves  a  translating  shield  with  a  shutter 
slot,  and  a  stationary  probe  on  which  a  MgO-layered  glass 
slide  (12  x  25  mm)  was  mounted  facing  downward  toward  the 
injector.  Initially,  the  translating  shutter  was  to  the  left  of  the 
slide  (when  viewed  in  Fig,  2)  so  that  the  slide  was  protected 
from  initial  disturbances  by  the  shield.  Once  the  spray  was 
burning  steadily,  the  wire  retainer  of  the  shield  was  fused  by 
an  electrical  current,  so  that  the  unbalanced  pressure  force 
accelerated  the  shutter  across  the  probe.  The  shutter  opening 
(3,2  mm  wide)  traversing  across  the  probe  allowed  drops  to 
impact  on  the  MgO  layer,  leaving  cavities  in  the  layer.  The 
motion  of  the  shutter  was  stopped  by  a  robber  cushion  after  it 
had  completely  traversed  the  slide,  leaving  the  slide  protected 
once  again  by  the  shield.  The  motion  of  the  shutter  was 
recorded  by  a  linear  displacement  transducer  (not  shown  in 
Fig.  2)  so  that  the  time  of  passage  across  various  portions  of 
the  slide  was  known, 

The  MgO  layer  was  produced  by  passing  the  slide 
back  and  forth  through  the  plume  of  an  Mg  flame  from  a  3  x 
25  mm  Mg  strip  burning  in  air.  Wetting  the  slide  with  water 
prior  to  coating  it  was  found  to  produce  a  layer  with 
reasonably  good  resistance  to  shock  disturbances  and  the 
high-temperature  environment.  The  relationship  between  drop 
and  impact  crater  size  was  found  by  calibrations  using  a 
Berglund-Liu  monodisperse  drop  generator.  The  calibrations 
yielded  a  ratio  of  roughly  2.0  between  the  impact  crater  and 
drop  diameters  for  drops  up  to  600  pm  in  diameter.  Crater 
diameters  were  measured  using  a  Unitron  inverted  metal¬ 
lurgical  microscope  at  X50  magnification,  with  a  linear 
micrometer  reticle  scale  on  the  eyepiece.  Only  drops  having 
diameters  larger  than  10  pm  could  be  collected  and  sized; 
therefore,  measured  distributions  are  biased  toward  larger 
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drop  diameters,  although  this  does  not  affect  estimates  of 
SMD  appreciably. 

The  slide  impactor  was  positioned  at  various  distances 
from  the  injector  exit,  by  moving  the  injector.  Due  to 
excessively-high  liquid  fluxes,  which  resulted  in  overlapping 
craters  and  washing  out  of  the  slide  at  the  shortest  feasible 
shutter-opening  times,  it  was  not  possible  to  make  measure¬ 
ments  closer  than  135  mm  from  the  injector.  By  accumulating 
data  for  several  slides  at  a  particular  condition,  50-200  drops 
were  sized  to  find  drop  size  distributions  and  the  SMD. 
Knowing  the  drop  sizes  and  slide  exposure  times,  these 
measurements  could  be  used  to  find  liquid  fluxes  and  the  total 
liquid  flow  rate  over  the  crossection  of  the  spray. 

Experimental  uncertainties  of  the  slide  impaction 
measurements  were  dominated  by  difficulties  of  overlapping 
impaction  craters  in  the  dense  portion  of  the  spray  near  the 
axis,  and  limitations  of  finite  sample  sizes.  In  particular, 
excessive  liquid  fluxes  often  required  extrapolation  of  findings 
away  from  the  axis  to  infer  spray  properties  over  the 
crossection.18  The  experimental  uncertainties  (95  percent 
confidence)  were  large  as  a  result  and  are  estimated  as  follows: 
SMD  less  than  25  percent,  liquid  fluxes  less  than  40  percent, 
and  liquid  flow  rates  only  within  a  factor  of  two  since  flow 
rates  near  the  axis  involved  extrapolation.  These  uncertainties 
are  substantially  greater  than  earlier  work  in  evaporating  air 
atomized  sprays  at  atmospheric  pressure,17  reflecting  the 
increased  problems  of  the  slide  impaction  technique  in  high- 
pressure  and  high- temperature  environments  within  dense 
pressure-atomized  sprays. 

Test  Conditions 

Test  conditions  are  summarized  in  Table  1.  Injector 
Reynolds  numbers  are  high  enough  to  yield  turbulent  flow  at 
the  injector  exit,  well  within  the  atomization  breakup  regime 
where  the  multiphase  mixing  layer  begins  to  form  right  at  the 
injector  exit  l2‘ 16  Combustor  pressures  were  in  the  range  4.6- 
4,9  MPa,  which  is  high  enough  for  drop  combustion  to  be 
dominated  by  surface  gasification  rather  than  microexplosions, 
based  on  single-drop  combustion  studies.3 

Theoretical  Mcihoda 
flas.-Ehasg  Formulation 

In  view  of  the  large  monopropellant  burning  rates  for 
present  est  conditions,  ca.  10  mm/s,  the  heterogeneous 
premixed  flame  that  covers  ail  liquid  surfaces  is  relatively  thin, 
ca.  1  pro.2  Thus  the  flow  largely  involves  two  states, 
unreacted  liquid  at  its  injection  temperature,  and  completely 
reacted  gaseous  combustion  products,  if  the  ambient  state  aiso 
consists  of  products  of  adiabatic  combustion.  Since  drop 
combustion  properties  were  not  strongly  influenced  by 
ambient  temperatures,3  this  approximation  is  made  in  the 
following.  Then,  both  phases  have  constant  densities  and 
there  is  no  need  to  treat  scalar  mixing,  which  substantially 
simplifies  the  formulation.  The  present  approach  involves 
Eulerian  calculations  for  the  properties  of  the  gas  phase,  and 
Lagrangian  calculations  of  drop  trajectories  within  the  gas 
phase.  Drop  source  times  in  the  governing  equations  for  the 
gas  phase  account  for  interphase  transport  of  mass  and 
momentum. 

Major  assumptions  of  the  gas  phase  formulation  are  as 
follows:  axisymmetric  and  steady  (in  the  mean)  flow  with  no 
swirl,  boundary-layer  approximations  apply,  constant  density 
gas,  and  negligible  effects  of  buoyancy.  The  first  of  these  is  a 
condition  of  the  experiments  while  pressure-atomized  sprays, 


like  single-phase  jets,  generally  satisfy  ihe  boundaiy  layer 
approximations.  The  constant  density  approximation  follows 
from  the  assumption  of  thin  heterogeneous  flames  as 
discussed  earlier.  Finally,  injector  velocities  are  ca.  70  m/s 
while  the  multiphase  flow  region  was  only  roughly  300  mra 
long;  therefore,  it  is  reasonable  to  neglect  effects  of  buoyancy. 
Similar  to  Solomon  et  al.,17  the  gas  phase  velocity  field  was 
found  using  a  k-e  turbulence  model  that  has  been  well 
calibrated  for  constant  density  turbulent  round  jets.20-22 

Boundary  conditions  along  the  liquid  core  are  specified 
as  described  later.  Within  the  multiphase  flow  region,  drop 
properties  were  found  by  solving  Lagrangian  equations  of 
motion  for  the  size  and  trajectories  of  a  statistically  significant 
sample  of  individual  drops  and  then  computing  source  terms 
for  mass  and  momentum  exchange  from  drops  in  the 
governing  equations  for  the  gas  phase.  This  involves  dividing 
the  drops  into  n  groups  (defined  by  initial  position,  size  and 
velocity)  and  then  computing  their  subsequent  motion  in  the 
flow.  Based  on  recent  findings  for  pressure-atomized  non¬ 
evaporating  sprays,4*16  effects  of  collisions  don’t  appear  to  be 
important  even  in  the  dense  portion  of  pressure -atomized 
sprays  and  were  ignored.  Similarly,  effects  of  adjacent  drops 
are  not  likely  to  influence  the  heterogeneous  flames,  since  they 
are  thin,  and  were  ignored  as  well.  Effects  of  secondary  drop 
breakup  may  be  important  for  present  sprays,16  however,  this 
still  was  neglected  due  to  the  lack  of  information  concerning 
this  phenomenon.  Since  initial  drop  size  distributions  were 
based  on  measurements  downstream  of  the  region  where 
secondary  breakup  would  be  important,  the  effect*  was 
accommodated  in  any  event.  Additionally,  direct  effects  of 
drops  on  turbulence  properties  (called  turbulence  modulation) 
were  ignored  since  the  bulk  of  the  flow  is  relatively  dilute  and 
methods  to  treat  these  phenomena  are  not  well  established.22 
Finally,  drop/turbulence  interactions  are  known  to  influence 
the  structure  of  noncombusling  sprays,1722  therefore,  this 
was  evaluated  for  combusting  monopropellant  sprays  by 
carrying  out  both  DSF  and  SSF  computations,  which  ignore 
and  consider  drop/turbulence  interactions  as  noted  earlier. 

FomiulatiQQ* 

Under  present  assumptions,  the  gas  phase  governing 
equations  are  the  same  for  the  DSF  and  SSF  formulations,  and 
can  be  written  as  follows:1722 

r  5/dx(p  u  0)  +  5/9r(rp  v  <j>)  =  d/9r((rpi/a$)ity/<fr’)  +  rS$  (4) 

where  there  is  no  need  to  consider  Favre-averaged  variables, 
since  the  density  of  the  gas  phase  is  constant.  The  parameters 
$  and  S$  appearing  in  Eq.  (1)  are  summarized  in  Table  2, 
along  with  all  empirical  parameters  used  in  the  calculations. 
Since  the  Reynolds  numbers  of  present  flows  are  relatively 
high,  laminar  transport  has  been  ignored.  The  source  terms 
for  =  1  (conservation  of  mass)  and  u  in  Table  2  are  the  drop 
source  terms  found  by  computing  the  net  change  of  mass  and 
momentum  of  each  drop  group  i  passing  through  compu¬ 
tational  cell  j,  having  a  volume  Vj.  The  boundary  conditions 
for  Eq,  (1)  involve  a  constant  property  ambient  environment 
and  symmetry  at  the  axis  (aside  from  the  liquid  core  region 
which  will  be  taken  up  later),  as  follows: 

r  — » »,  4>  =  0,  r  =  0,  Stydr  =  0  (2) 

Drop  Formulation 

Drop  trajectories  were  computed  similar  to  earlier 
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work.3  Drops  were  assumed  to  be  surrounded  by  gas  in  the 
multiphase  flow  region.  Effects  of  drop  heatup  were  ignored 
and  the  drop  radius  was  assumed  to  decrease  at  10  mm/s, 
based  on  drop  burning  rate  measurements  at  the  pressure  of 
present  tests.3  The  high  monopropellant  burning  rate  implies 
a  thin  decomposition  flame  near  the  surface,  well  within  the 
boundaries  of  the  drop  flow  field:  therefore,  ambient  gas 
properties  were  used  to  estimate  drop  drag  and  effects  of 
forced  convection  on  drop  burning  rates  were  ignored. 
Similarly,  the  eventual  loss  of  the  linear  burning  rate  law  for 
monopropellants  (Kp  =  -  drp/dt  =  constant)  for  small  drops 
was  ignored,  since  this  would  occur  for  drop  diameters  of  ca. 

I  jim.3  Other  assumptions  were  similar  to  earlier  separated 
flow  models:17'22  the  flow  field  around  the  drops  was 
assumed  to  be  quasisteady;  virtual  mass,  pressure-gradient. 
Basset  history  and  gravitational  forces  were  ignored;  drop 
swelling  was  ignored;  and  drop  drag  was  estimated  using  the 
standard  drag  correlation  for  solid  spheres. 

Under  these  assumptions,  the  governing  equations  for 
drop  motion  arc  as  follows:22 

d  Xpi/dt  =  upi,  i=l,2,3  (3) 

d  dp/dr  =  -  2  Kp  (4) 

dupi/dt  =  (3pCD/4  dpPp)(uj-uPi)  lu-Upl,  i=l,2,3  (5) 

where  a  Cartesian  coordinate  system  has  been  used  and  lu-upl 
denotes  the  magnitude  of  the  relative  velocity  between  the  drop 
and  the  gas.  The  drag  coefficient  expression  was  identical  to 
past  work.22 

Equations  (3)  and  (5)  were  solved  using  time-averaged 
gas  phase  velocities  for  the  DSF  model.  In  contrast,  the  SSF 
model  involves  an  approximate  approach  to  treat 
turbulence/drop  interactions  that  has  been  reasonably 
successful  during  past  evaluations.17-22  This  involves 
computing  drop  trajectories  as  they  encounter  a  succession  of 
turbulent  eddies.  The  properties  of  each  eddy  are  assumed  to 
be  uniform  and  to  change  randomly  from  one  eddy  to  the  next 
The  gas  velocities  within  the  eddies  were  found  by  making  a 
random  selection  from  the  probability  density  function  (PDF) 
for  velocity,  assuming  an  anisotropic  Gaussian  PDF  with 
stream  wise  and  crosstream  variances  of  k  and  k/2.  A  drop 
was  assumed  to  interact  with  an  eddy  as  long  as  the  time  of 
interaction  did  not  exceed  the  eddy  lifetime  i«,  or  the 
displacement  of  the  drop  did  not  exceed  the  characteristic  eddy 
size,  L«.  These  quantities  were  estimated,  as  follows:22 

U  =*  te  -  IV(2k/3)W  (6) 

All  quantities  needed  for  these  procedures  can  be  found  from 
the  gas  phase  turbulence  model 

Initial  Conditions 

Initial  conditions  were  specified  in  the  region  of  the 
liquid  core  as  illustrated  in  Fig.  3.  The  present  approach  is 
somewhat  similar  to  Chatwani  and  Bracco.13  Since  the  LHF 
approach  is  expected  to  give  reasonable  estimates  of  flow 
properties  at  large  liquid  volume  fractions, 14  the  earlier  LHF 
raonopropeliant  combustion  model  was  used  to  estimate  the 
flow  field  in  the  near-injector  region,3  The  boundary  of  the 
liquid  core  was  then  defined  by  locations  where  the  time- 
averaged  Uquid  volume  fraction  was  0.997,  with  mean  and 
fluctuating  streamwise  and  radial  velocities  identified  at  40 
locations  from  the  LHF  predictions  (see  Table  3  for 


representative  data  at  these  boundary  points).  The  velocities 
of  drop  groups  emanating  from  these  positions  were  then 
estimated  by  either  the  mean  velocities  (DSF  approach)  or  by  a 
random  selection  from  the  velocity  PDFs  (SSF  approach). 
The  drop  size  distributions  at  the  boundary  of  the  liquid  core 
were  extrapolated  back  from  the  position  nearest  to  the  injector 
where  they  were  measured,  assuming  an  average  drop  velocity 
and  the  linear  burning  rate  law,  see  Table  3  for  the  distribution 
used  for  d  =  0.6  mm.  This  translation  procedure  biases  out 
drops  that  bum  completely  before  reaching  the  first  meas¬ 
urement  location,  however,  measured  spray  properties  are 
dominated  by  large  drops  and  this  effect  is  expected  to  be 
small. 

Drop  size  estimates  obtained  at  the  boundary  of  the 
liquid  core  in  this  manner  appear  to  be  reasonable.  Using 
either  aerodynamic  breakup  theory,  due  to  Reitz  and  Bracco,12 
or  estimates  based  on  maximum  stable  drop  sizes  for  secon¬ 
dary  breakup  near  the  liquid  surface,  average  drop  sizes  near 
the  liquid  surface  for  pressure-atomized  nonevaporating 
sprays  are  correlated  by; 16 

dp*Vg  3  CfiO/(pu^)  (7) 

where  Os  is  roughly  10  and  40  for  slug  and  fully-developed 
flows  at  the  injector  exit.  Taking  the  latter  value  as 
representative  of  present  test  conditions  yields  dPavg  =  250 
(im,  which  is  comparable  to  present  estimates  of  SMD  along 
the  boundary  of  the  liquid  core,  see  Table  3. 

Computations 

Calculations  for  the  continuous  phase  were  performed 
with  a  modified  version  of  the  GENMJX  algorithm,^  with  a 
computational  grid  similar  to  past  work.17-22  Roughly  1000 
drop  groups  were  used  for  both  the  DSF  and  SSF  predictions, 
with  Eqs.  (3)-(5)  integrated  by  a  second-order  Runge-Kutta 
algorithm.  Numerical  accuracy  was  evaluated  by  doubling 
and  halving  the  grid  size  and  the  number  of  drop  groups, 
indicating  that  the  results  reported  here  are  numerically  closed 
within  3  percent. 

Results  and  Discussion 

Drop  Size  Properties 

Predicted  and  measured  SMD  over  the  spray 
crossection  are  plotted  as  a  function  of  x/d  in  Fig.  4  for 
injector  diameters  of  0.3  and  0.6  mm.  Predictions  include 
results  from  both  the  DSF  and  SSF  models.  The  first 
observation  is  that  liquid  extends  substantially  beyond  the 
region  where  liquid  was  seen  during  spray  shadowgraph 
observations  of  liquid  volume  fractions,3  e.g.,  x/d  =  500- 
1000  for  d  =»  0.3  and  0.6  mm,  as  opposed  to  x/d  ca.  350  for  a 
similar  range  of  pressures  and  injector  diameters.  It  is  clear 
that  earlier  spray  shadowgraphs  did  not  record  drops  within 
the  dilute  portions  of  the  spray  far  from  the  injector  exit  and 
really  don't  offer  a  quantitative  evaluation  of  predictions. 
Another  feature  of  the  measurements  is  that  SMD  does  not 
scale  with  x/d  as  d  is  changed,  which  would  be  indicative  of 
mixing-controlled  behavior  that  might  satisfy  the  LHF 
approximation.  Instead,  the  scaling  more  closely  approx¬ 
imates  completion  of  combustion  at  a  fixed  distance  from  the 
injector  with  x/d  at  the  end  of  the  liquid-containing  region 
increasing  from  ca.  580  to  ca,  980  as  d  is  reduced  from  0.6  to 
0.3  mm.  Since  initial  drop  sizes  and  velocities  of  the  two 
injectors  are  estimated  to  be  nearly  the  same,  this  behavior  is 
indicative  of  strong  effects  of  separated  flow. 

Predictions  using  the  DSF  and  SSF  models  are  not 
very  different  from  each  other  in  Fig.  4  and  are  in  reasonably 
good  agreement  with  the  measurements.  The  agreement 
between  predictions  and  measurements  is  best  toward  the 
downstream  end  of  the  spray  where  the  slide  impaction 


4 


75 


flow. 


technique  was  more  reliable  due  to  lower  liquid  fluxes. 
Discrepancies  between  predictions  and  measurements  are 
largest  for  near  injector  positions  and  d  =  0.3  mm.  where  drop 
sizes  are  felt  to  be  biased  downward  to  difficulties  of 
measuring  large  drops  near  the  axis  due  to  high  rates  of  drop 
impaction  since  these  positions  were  relatively  close  to  lhe 
injector  exit.  The  predictions  were  strongly  dominated  by 
effects  of  drop  banting  rates,  while  initial  drop  size 
distributions  were  calibrated;  therefore,  present  findings 
suggest  that  drop  burning  rates  observed  in  the  combusting 
sprays  are  consistent  with  the  single  drop  results  reported 
earlier.3 

Measured  and  predicted  (both  DSF  and  SSF  methods) 
drop  size  distributions  arc  illustrated  in  Figs.  4  and  5  for  the 
two  injector  diameters.  As  noted  earlier,  the  measurements  at 
positions  nearest  the  injector  were  extrapolated  to  the  injector 
exit  in  an  approximate  manner  to  provide  initial  conditions  for 
the  predictions.  Reversing  the  process  to  yield  the  predictions 
illustrated  in  Figs.  5  and  6  does  not  exactly  reproduce  the 
distributions  at  the  smallest  x/d  due  to  different  rates  of 
deceleration  experienced  by  various  drop  groups.  This  yields 
somewhat  whimsical  variations  of  predicted  drop  size  distri¬ 
butions  at  the  smallest  x/d,  due  to  limited  drop  size  statistics, 
but  the  predictions  are  still  in  reasonably  good  agreement  with 
the  measurements. 

Aside  from  limitations  in  resolving  drop  sizes  near  10 
Hm,  the  subsequent  evolution  of  the  drop  size  distributions 
with  increasing  x/d  is  also  predicted  reasonably  well  in  Figs.  5 
and  6.  In  particular,  the  distributions  become  skewed  toward 
smaller  drop  sizes  as  the  end  of  the  liquid-containing  region  is 
approached,  which  is  represented  reasonably  well  by  the 
predictions. 

Liquid  Flow  Properties 

As  noted  earlier,  measurements  of  liquid  flow  rates 
and  fluxes  are  less  accurate  than  measurements  of  SMD, 
however,  in  spite  of  these  limitations,  the  results  help  provide 
a  more  complete  picture  of  the  structure  of  the  sprays. 
Measured  and  predicted  normalized  liquid  flow  rates  over  the 
spray  crossection  are  plotted  as  a  function  of  x/d  in  Fig.  7  for 
the  two  injectors.  Measurements  at  one  position  for  each  of 
the  injectors  are  significantly  higher  than  the  rest  and  are  felt  to 
be  anomalous  due  to  problems  with  the  slide  impactor  in  the 
period  where  these  results  were  obtained.  The  remaining 
results  show  the  expected  trend  of  decreasing  liquid  flow  rates 
with  increasing  x/d.  Although  liquid  is  observed  for  x/d  > 
500,  the  flow-rate  is  generally  less  than  10  percent  of  the 
original  flow  rate  at  the  injector  exit  This  comes  about  due  to 
both  consumption  of  drop  liquid,  as  well  as  reduced  drop 
velocities  in  the  relatively  slow  gas  flow  far  from  the  injector. 

The  comparison  between  predictions  and  measure¬ 
ments  for  liquid  flow  rates  in  Fig.  7  is  not  as  good  as  for  SMD 
in  Fig.  4.  This  is  typical  of  evaluations  of  spray  models  using 
liquid  flow  rates  and  fluxes:  these  quantities  are  difficult  to 
measure  accurately  in  pressure  atomized  sprays  and  are 
particularly  sensitive  to  uncertainties  in  model  predic- 
nons.*3-17-22  as  before,  there  is  little  to  choose  between  the 
DSF  and  SSF  methods,  with  the  latter  yielding  a  somewhat 
shorter  spray  length.  This  comes  about  due  to  higher  mean 
drag  when  instantaneous  gas  velocities  are  used  to  estimate 
drop  trajectories  (since  drag  forces  are  nearly  quadratic 
functions  of  relative  velocities  for  present  conditions),  as  well 
as  effects  of  turbulent  dispersion  which  transports  drops 
laterally  into  slower  moving  gas  near  the  edge  of  the  flow. 
However,  effects  of  turbulent/drop  interactions  are  much 
reduced  for  the  present  monopropeUant  sprays  in  comparison 
to  evaporating  sprays,17  since  the  rate  of  drop  gasification  is 
not  influenced  by  the  radial  position  of  the  drops  within  the 


An  indication  of  the  width  of  the  sprays  can  be  seen 
from  the  plot  of  normalized  liquid  mass  flux  in  Fig.  8  for  d  = 

0.3  mm.  Results  for  d  =  0.6  mm  are  similar.18  These  results 
are  ploited  in  terms  of  the  radial  similarity  variable  for 
turbulent  jets,  r/x,  so  that  the  actual  width  of  the  flow  can  be 
seen.  However,  similarity  of  liquid  mass  fluxes  in  these 
coordinates  is  neither  expected  nor  observed.  Spray  widths 
extend  to  r/x  ca.  0.12  at  x/d  =  460,  but  the  liquid  becomes 
progressively  more  confined  near  the  axis  at  larger  distances 
from  the  injector. 

Liquid  fluxes  are  difficult  to  predict  accurately,  as 
noted  earlieT,13’17'22  and  there  is  only  qualitative  agreement 
between  predictions  and  measurements  in  Fig.  8.  In  fact, 
predictions  underestimate  the  spray  length  for  these  con¬ 
ditions,  see  Figs.  4  and  7,  so  that  no  predictions  are  available 
for  comparison  with  the  measurements  at  the  highest  position. 
Consideration  of  turbulent  dispersion  using  the  SSF  model 
does  not  yield  a  wider  liquid  flux  distribution.  Instead, 
separated  flow  effects  dominate  the  process,  causing  small 
drops  to  disappear  more  quickly  in  the  low  velocity  region 
near  the  edge  of  the  flow. 

Additional  Predictions 

The  comparison  between  predictions  and 
measurements  indicates  that  both  the  DSF  and  SSF  models 
provide  the  general  features  of  the  combusting  monopropeUant 
sprays,  particularly  for  drop  size  properties  where  predictions 
are  less  sensitive  to  model  parameters  and  measurements  are 
more  reliable.  Thus,  additional  computations  were  performed 
to  help  gain  insight  concerning  spray  properties  where  meas¬ 
urements  were  not  available.  This  included  predictions  of 
phase  velocities  and  the  sensitivity  of  predictions  to  variations 
of  model  parameters  and  spray  operating  conditions. 

Predicted  phase  velocities  for  present  test  conditions 
arc  plotted  in  Fig.  9.  NumbcT-averaged  drop  velocities  and 
time-averaged  gas  velocities  along  the  axis,  found  from  the 
DSF  model,  are  plotted  as  a  function  of  distance  from  the 
injector  ford  =  0.3  and  0.6  mm  Effects  of  separated  flow  are 
evident,  with  liquid  velocities  decaying  much  more  slowly 
than  gas  velocities.  The  relative  velocities  are  larger  for  d  =» 
0.3  mm,  since  the  smaller  diameter  injector  involves  higher 
rates  of  deceleration  of  the  gas,  which  cannot  be  followed  by 
the  drops  due  to  their  limited  response  properties.  In  fact, 
liquid  velocities  remain  nearly  constant  until  x/d  >  100,  since 
both  gas  velocities  and  drop  diameters  are  largest  in  this 
region.  Farther  from  the  injector,  drop  velocities  decrease 
rapidly  toward  gas  velocities,  due  to  increased  drop  response 
as  the  drops  become  smaller.  This  increases  drop  residence 
times  for  a  particular  penetration  length  and  makes  estimates  of 
properties  like  liquid  flow  rates  and  fluxes  at  large  x/d  very 
sensitive  to  model  parameters,  as  noted  earlier. 

The  sensitivity  of  predictions  to  variations  of  model 
parameters  and  spray  operating  conditions  can  be  seen  more 
quantitatively  from  the  results  summarized  in  Table  4.  The 
percent  increase  of  x/d  at  the  points  where  SMD  =  0  and  m/mo 
3  1/2  are  tabulated  for  various  changes  of  parameters.  These 
results  involve  use  of  the  DSF  model  (except  for  the  last  item) 
and  d  =  0.6  mm,  however,  they  are  typical. 

As  expected,  results  in  Table  4  show  that  the  drop 
burning  rate  has  a  strong  effect  on  the  axial  penetration  of 
spray  liquid,  with  a  nearly  one-to-one  correspondence  for 
changes  of  drop  burning  rates  by  a  factor  of  two.  Results  are 
less  sensitive  to  initial  SMD,  however,  since  reduction  of  drop 
sizes  increases  momentum  exchange  to  the  gas  causing  gas 
velocities  to  increase.  The  resulting  larger  drop  velocities  tend 
to  compensate  for  smaller  drop  lifetimes  when  evaluating 
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spray  penetration,  Raising  the  ambient  density  by  a  tactor  of 
10,  for  a  fixed  initial  SMD,  only  decreases  penetration  lengths 
by  30-40  percent,  largely  due  to  increased  entrainment  rates 
which  reduce  gas  phase  velocities  in  the  spray.  Varying 
ambient  density  is  actually  likely  to  have  a  greater  etfect, 
however,  by  reducing  drop  sizes  through  the  drop  breakup 
criterion  e«  Eq.  (7).  Based  on  effects  of  SMD  in  Table  4, 
this  would' considerably  shorten  spray  penetration  lengths 
from  results  measured  during  the  present  investigation. 
Effects  of  injector  diameter  variations  follow  from  the  nearly 
constant  spray  penetration  length  for  various  injector 
diameters,  i.e„  increasing  the  diameter  by  a  factor  of  ten 
reduces  the  penetration  length,  in  terms  of  x/d,  by  nearly  the 
same  ratio.  Finally,  effects  of  including  turbulent  dispersion 
are  relatively  small,  particularly  in  view  of  uncertainties  in 
estimating  initial  conditions  for  separated  flow  calculations. 


References 

'Kounalaids,  M.  E.  and  Faeth,  G.  M.,  "Combustion 
of  HAN-Based  Liquid  Monopropellants  near  the  Thermo¬ 
dynamic  Critical  Point,"  Comhust.1  Flame.  VoL  74,  November 
1988,  pp.  179-192. 

2Lee,  T,-W.,  Gore,  J.  P.,  Faeth,  G.  M.  and  Birk,  A., 
"Analysis  of  Combusting  High-Pressure  Monopropcllant 
Sprays,"  Combust.  Sci,  and  Tech..  Vol.  57,  1988,  pp.  OS- 
112. 


^Lee,  T.-W.,  Tseng,  L.-K.  and  Faeth,  G.  M., 
"Separated-:Flow  Considerations  for  Pressure-Atomized 
Corn-busting  Monopropellant  Sprays,"  J.  Pron.  Power,  in 
press. 


Conclusions 

Major  conclusions  of  the  study  are  as  follows: 

1 .  Present  slide  impaction  measurements  show  that  the 
liquid-containing  region  is  2-4  times  longer  than  found 
during  earlier  measurements  using  flash  shadowgraphs 
by  Lee  et  al.,3  due  to  the  inability  of  the  shadowgraph 
approach  to  detect  dilute  portions  of  the  spray. 

2.  Present  measurements  of  combusting  spray  properties 
show  that  tnonopropellant  drop  burning  rates  within 
the  sprays,  •  drp/dt  ca.  10  mm/s,  are  consistent  with 
earlier  single  drop  burning  rate  measurements  at  the 
same  pressure,  ca.  5  MPa.  Due  to  its  dominance  of 
spray  combustion  properties,  additional  study  of  the 
drop  combustion  properties  of  HAN-based  monopro- 
peliants  is  recommended. 

3.  The  combusting  monopropellant  sprays  exhibited 
strong  effects  of  separated  flow,  with  the  length  of  the 
liquid  containing  region  being  nearly  the  same,  ca  300 
mm,  for  injector  diameters  of  0.3  and  0.6  mm  when 
initial  SMD  and  injector  exit  velocities  were  nearly  the 
same. 

4.  Estimated  mean  drop  sizes  near  the  injector  exit  were 
consistent  with  recent  measurements  for  nonevap¬ 
orating  pres  sure -atomized  sprays.16  This  suggests 
that  initial  mean  drop  sizes  (SMD)  can  be  estimated 
from  either  the  aerodynamic  breakup  theory  of  Reitz 
and  Bracco,12  or  analogous  secondary  breakup  cri¬ 
teria,16  through  Eq.  (7).  However,  this  approach  is 
only  provisional,  pending  additional  study  of  effects  of 
liquid  turbulence  at  the  injector  exit  and  approach  of 
the  flow  to  the  LHF  limit 

5.  Predictions  using  the  DSF  and  SSF  models  were  in 
reasonably  good  agreement  with  each  other  and  with 
the  measurements,  particularly  in  view  of  the  relatively 
large  uncertainties  of  estimates  of  initial  drop  prop¬ 
erties  and  the  slide  impaction  measurements,  Thus, 
the  DSF  approach  is  favored,  since  it  requires  less 
extensive  computations.  The  relatively  small  effects  of 
drop/tuxbulence  interactions  follows  since  the  gas 
phase  has  constant  properties  and  drop  burning  rates 
are  relatively  independent  of  the  properties  of  their 
surroundings  for  combusting  tnonopropellant  sprays. 
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Table  1  Summary  of  test  conditions^ 


Injector 

Diameter1* 

(mm) 

.Ambient 

Pressure0 

Re 

(MPa) 

Pressure 

Drop 

Oh 

(MPa) 

Injector 

Velocity11 

Wef 

(m/s) 

Drop  Size  Measurements: 

0.58 

4.89 

3.13 

65.5 

7800 

0.030 

54200 

0,31 

4.69 

3.96 

73.8 

4700 

0.041 

36700 

Liquid  Flux  Measurements: 

0.58 

4.63 

3.56 

70  0 

8300 

0.030 

61700 

0.33 

4.75 

3.70 

71.4 

4800 

0.040 

36500 

Combustion  of  LGP  1845,  consisting  of  63.2%  HAN,  20% 
TEAN  and  16.8%  H^O  (by  mass).  Adiabatic  combustion 
yields  product  temperature  of  2150  K,  consisting  of  69.2% 
H20,  12.9%  CO2  and  17.4%  N2  (by  volume).  Liquid 
properties:  pf  =  1454  kg/m3,  |lf  =  0.0071  kg/ms,  c  =  0.067 
N/m. 

bInjector  length-to-diameter  ratio  of  42,  yielding  fully- 
developed  turbulent  pipe  flow  at  exit 
Combustion  product  environment  at  2790  K  and  consisting 
of  19.8%  H2O,  38.3%  N2  and  40.9%  A r  (by  volume). 
dFor  a  flow  coefficient  of  unity. 


Table  2  Source  terms  and  empirical  constants 
for  separated-flow  predictions 
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s* 

l 

n 

Z  ni  (miin-miouiyVj 
i=l 
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n 

2  ni((miUpi)in  -  (miUpi)out)j/Vj 

i=l 

k 

Pi(3u/3r)2  •  pe 

£ 

(Ceitii(3u/5r)2  -  Ce2pe)e/k 

Cp 

Cel  C*2  °k  Kb(mm/s) 

Cnpk2/e 

0.09 

1.44  1.87  1.0  1.3  10 
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Table  3  Initial  conditions  for  separaud-flow 


predictions® 


185 

190 

195 

200 

205 

210 

3 Phvin  oi-T 

1.20 

2.91 

5,32 

7.92 

9  73 

8.73 

is  /JiuhJIk.J _ r*  n  , 

215 

220 

225 

230 

235 

240 

7.72 

6,92 

6.42 

5.72 

5.22 

4.71 

245 

250 

255 

260 

265 

275 

4.21 

3.91 

3.51 

3.21 

2.81 

2.31 

285 

295 

305 

315 

325 

335 

2.01 

1.71 

1.40 

1.10 

0.80 

0.50 

-.-r  uuuiuuuon  tor  a  ue  mm  diameter 

Percentage  of  total  drop  number  flow  rate. 

injector. 

Table  4  Parameter  sensitivity  of  separated-flow 
predictions1 


Input  Factor  of 

Percent  increase  of  x/d 

Parameter  Variation 

SMD  =  0 

m  =  m<y2 

Drop  burning  rate 

1/2 

41 

56 

Drop  burning  rate 

2 

-  32 

•  40 

Initial  SMD 

1/2 

-28 

•38 

Ambient  density 

10 

-42 

-  30 

Injector  diameter 

10 

-80 

•  84 

Turbulent  dispersion** 

... 

-  9 

-  11 

Results  ford  =  0.6  mm 
noted 

using  the 

DSF  model 

except  as 

^e  effect  of  including 

turbulent 

dispersion  is  shown. 

Sketch  of  spray  combustion  apparatus. 
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Fig,  2  Sketch  of  slide  impactor  assembly. 


Fig.  3  Specification  of  initial  conditions  for  separated- 
flow  computations. 
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Fig.  5  Predicted  and  measured  drop  size  distributions  at 
various  stream  wise  distances,  d  =■  0.6  mm. 


Fig.  6  Predicted  and  measured  drop  size  distributions  at 
various  stream  wise  distances,  d  =  0,3  mm. 


Fig.  4  Predicted  and  measured  SMD  variation  with 
streamwise  distance. 
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